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Abstract — The  dynamic  characteristics  of  flexible  composite  beams  are  controlled  by  heating  sets 
of  shape  memory  alloy  (NITINOL)  fibers  embedded  along  the  neutral  axes  of  these  beams.  The 
activation  of  the  shape  memory  effect  of  the  fibers  increases  the  elastic  energy  and  enhances  the 
stiffness  of  the  composite  beams.  With  such  capabilities,  the  vibration  modes  of  the  beams  can  be 
tailored  and  shifted  away  from  the  excitation  frequency  band  in  order  to  avoid  undesirable 
vibrations. 

Emphasis  is  placed,  in  the  present  study,  on  the  effect  of  intentional  electrical  heating  of  a 
selected  subset  of  the  NITINOL  fibers  on  the  overall  dynamics  of  the  beams.  The  effect  of  the 
associated  thermal  energy  propagating  through  the  composite  on  the  unintentional  thermal 
aaivation  of  additional  subsets  of  the  NITINOL  fibers  is  accounted  for.  Such  an  effect  is  not  only 
significant  but  also  essential  to  the  thorough  undeistanding  of  the  operation  of  NITINOL- 
reinforced  composites. 

Finite  element  models  are  developed  to  describe  the  interaction  between  the  thermal  and 
dynamic  characteristics  of  the  NITINOL  composites  as  well  as  the  interaction  between  the 
intentional  and  unintentional  activation  of  the  NITINOL  fibers.  The  models  are  experimentally 
validated  and  close  agreement  is  obtained  between  the  theoretical  predictions  and  the  experimental 
results.  The  mathematical  models  and  procedures  described  in  this  paper  provide  an  invaluable 
means  of  predicting  realistic  performance  of  NITINOL-reinforced  composites. 


NOMENCLATURE 


[A]  interpolating  function  of  beam  defleaion 

[A,]  /th  element  of  [A] 

A„  cross-sectional  area  of  the  beam 

[B]  matrix  of  the  first  derivatives  of  the  nodal  interpolating  functions 

(CJ,  [D]  first  and  second  derivatives  of  the  interpolating  function  of  beam  deflection 
E„  Young’s  modulus  of  the  beam 

[F]  vector  of  external  loads  aaing  on  the  beam 
h  convective  heat  transfer  coefficient 

I„  area  moment  of  inertia  of  the  beam 

k  thermal  conductivity  of  the  beam 

f^f.:.3]  matrices  given  by  eqns  (16),  (17)  and  (18),  respectively 
[K,]  stiffness  matrix  of  the  beam  element 

ly ,  1.  direction  cosines  of  outward  normals  to  the  beam  boundaries 
L  length  of  the  beam  element  and  NITINOL  fiber 

M,  external  moment  acting  at  the  ith  node 

[M,]  mass  matrix  of  the  beam  element 
m,(i,j)  the  element  i,j  of  the  mass  matrix 
[yV'I  interpolating  function  of  the  beam  temperature 

N,  interpolating  function  of  the  /th  node 

p  number  of  vertices  of  the  element 

,  mechanical,  net  and  thermal  axial  loads  acting  on  the  beam 
[P'l  matrix  given  by  eqn  (19) 

q„  generalized  coordinate  of  the  nth  vibration  mode  of  the  NITINOL  fiber 

q„  generalized  acceleration  of  the  mh  vibration  mode  of  the  NITINOL  fiber 

Q  heat  flux  per  unit  area 

S,  ,  boundaries  of  the  NITINOL  fibers  and  the  beam,  respectively 

t  time 

Tq  initial  tension  in  a  NITINOL  fiber 

T,  total  tension  in  a  NITINOL  fiber 

V,  shear  force  acting  at  the  I'th  node 

w  transverse  deflection  of  the  beam  and  NITINOL  fibers 

X,  V,  c  Cartesian  coordinates  along  the  beam  neutral  axis  and  cross-section,  respectively 
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Creek  letters 

a  thermal  diffusivity  of  the  beam 

[J]  the  deflection  vector  of  the  beam  element 

dj  angular  deflection  of  the  /th  node 

0  temperature  at  any  location  (y,  z)  of  the  beam  cross-section 

©,  ambient  temperature 

[©']  vector  of  the  nodal  temperatures  of  the  element 

0,  mode  shape  of  the  /tth  mode 

aj„  natural  frequency  of  the  nth  mode. 

1.  INTRODUCTION 

Considerable  attention  has  been  devoted  recently  to  the  utilization  of  the  Shape  Memory 
Nickel-Titanium  alloy  (NITINOL)  in  developing  SMART  composites  that  are  capable  of 
adapting  intelligently  to  external  disturbances  (Ikegami  et  al.,  1990;  Rogers  et  al.,  1991; 
Baz  et  al.,  1990,  1991a, b).  Such  wide  acceptance  of  NITINOL  stems  from  its  unique 
behavior  when  it  is  subjected  to  particular  heating  and  cooling  strategies.  For  example, 
the  alloy  becomes  soft  when  it  is  cooled  below  its  martensite  transformation  temperature 
and  becomes  about  four  times  stiffer  when  it  is  heated  above  its  austenite  transformation 
temperature  (Funakubo,  1987).  Furthermore,  if  the  alloy  is  trained  to  have  a  particular 
shape  while  in  its  austenite  phase,  it  will  memorize  this  shape.  If  the  alloy  is  then  cooled  to 
its  martensite  phase  and  subjected  to  plastic  deformation,  it  will  return  to  its  memorized 
shape  when  it  is  heated  above  the  austenite  transformation  temperature.  The  phase  trans¬ 
formation  from  martensite  to  austenite  produces  significant  forces  as  the  alloy  recovers  its 
original  shape.  The  alloy  acts  as  an  actuator  transforming  thermal  energy  into  mechanical 
energy  (Perkins,  1975;  Duerig  et  al.,  1990).  Accordingly,  if  the  NITINOL  fibers  are 
embedded  inside  a  composite  matrix  at  optimal  locations,  they  can  be  used  to  control  the 
static  and  dynamic  characteristics  of  the  resulting  SMART  composite.  The  control  action 
is  generated  by  the  described  stiffening  of  the  NITINOL  fibers  and/or  the  shape  memory 
effect.  With  such  built-in  control  capabilities,  the  performance  of  the  SMART  composites 
can  be  optimized  and  tailored  to  match  changes  in  the  operating  conditions. 

It  is  therefore  the  purpose  of  this  study  to  develop  a  thorough  understanding  of  the 
fundamentals  governing  the  operation  of  NITINOL-reinforced  composite  beams.  The 
individual  contributions  of  the  composite  matrix,  the  NITINOL  fibers  and  the  shape 
memory  effect  to  the  overall  dynamic  performance  of  the  composite  beams  will  be 
determined.  Also,  the  influence  of  the  temperature  distribution  inside  the  composite, 
which  results  from  the  activation  of  a  small  subset  of  the  NITINOL  fibers,  on  the  overall 
performance  of  the  entire  beam  will  be  addressed  both  theoretically  and  experimentally. 
Such  an  important  interaction  between  the  thermal  and  dynamic  characteristics  of  the 
NITINOL-reinforced  composites  has  not  been  addressed  in  the  previous  analyses  of 
Rogers  et  al.  (1991)  and  Jia  and  Rogers  (1989).  In  these  studies,  the  NITINOL-reinforced 
composites  have  been  considered  to  operate  isothermally  even  though  the  activation  and 
deactivation  of  the  NITINOL  fibers  subject  these  composites  to  non-uniform  temperature 
fields.  Furthermore,  the  effects  of  intentional  electrical  heating  of  a  selected  subset  of  the 
NITINOL  fibers  and  the  associated  thermal  energy  propagating  through  the  composite  on 
the  unintentional  thermal  activation  of  additional  susbsets  of  the  fibers  have  not  been 
considered  by  Rogers  et  al.  (1991)  and  Jia  and  Rogers  (1989).  These  effects  significantly 
alter  the  dynamics  of  NITINOL-reinforced  composites,  particularly  those  made  of 
multi-lamina  where  the  intentional  electrical  activation  of  a  NITINOL  lamina,  by  an 
active  controller,  will  generate  enough  heat  to  activate  thermally  and  unintentionally  the 
adjacent  NITONOL  laminas.  The  phenomena  associated  with  intentional  electrical 
activation  and  the  associated  unintentional  thermal  activation  will  be  addressed,  in  detail, 
in  the  present  study. 

The  present  paper  is  organized  in  five  sections.  In  Section  1  a  brief  introduction  is 
given.  In  Sections  2  and  3  the  dynamic  and  thermal  models  of  the  NITINOL-reinforced 
composites  are  presented,  respectively.  The  experimental  behavior  of  a  single  and  two- 
layer  NITINOL-reinforced  composite  is  given  in  Section  4,  both  in  the  time  and  frequency 
domains.  Section  5  summarizes  the  results  and  the  conclusions  of  the  study. 


NITINOL-reinforced  composite  beams 
2.  THE  DYNAMICS  OF  NITINOL-REINFORCED  COMPOSITE  BEAMS 
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In  the  present  study,  the  NITINOL-reinforced  composites  are  made  by  embedding 
the  NITINOL  fibers  inside  vulcanized  rubber  sleeves  placed  along  the  neutral  axes  of 
these  composite  beams,  as  shown  in  Fig.  I.  In  this  arrangement,  the  NITINOL  fibers  are 
free  to  move  during  the  phase  transformation  process  in  order  to  avoid  degradation 
and/or  destruction  of  the  shape  memory  effect  which  may  result  when  the  fibers  are 
completely  bonded  inside  the  composite  matrix. 


Fig.  1.  A  schematic  drawing  of  the  cross-seaion  of  a  NITINOL-reinforced  composite  beam. 

The  dynamic  characteristics  of  this  class  of  NITINOL-reinforced  beams  are  obtained 
by  dividing  each  beam  into  finite  elements.  The  forces  typically  acting  on  any  of  these 
elements  are  displayed  in  Fig.  2  along  with  the  associated  nodal  displacements.  The 
stiffness  and  mass  matrices  of  each  beam  element  are  derived  in  Sections  2.1  and  2.2, 
respectively,  using  the  theory  of  Bemoulli-Euler  beams.  The  overall  stiffness  and  mass 
matrices  of  the  entire  beam  are  obtained  by  assembling  the  stiffness  and  mass  matrices  of 
the  individual  elements  (Paz,  1991). 


Fig.  2.  NITINOL-reinforced  beam  element  with  forces  and  resulting  displacements. 


2.1.  The  stiffness  matrix 

The  element  stiffness  [K,]  is  made  up  of  the  flexural  rigidity  of  the  beam,  the 
geometric  stiffness  that  accounts  for  the  axial  and  thermad  loading  as  well  as  the  siiffnc'i 
imparted  by  the  elasticity  of  the  NITINOL  fibers.  The  combined  stiffness  of  the  element 
is  obtained  using  the  principle  of  conservation  of  energy  and  equating  the  work  done  by 
external  loads  to  the  strain  energies  stored  in  the  element  as  follows: 

Sum  of  work  done  by  external  loads  =  Sum  of  stored  strain  energies 


or 


^[(5]^(F]  +  PJ2  (dw/cix)2  dx  P,/2 


(dw/dx)'  dx 


=  £,„/„/2  (d^w/dx  V  cLf  +  7;/2  (dw/dx)‘dx 


(1) 


where  the  terms  on  the  left-hand  side  denote  the  work  done  by  the  transverse  loads  and 
moments,  the  axial  mechanical  loads  and  axial  thermal  loads  P, ,  respectively.  The 
terms  on  the  right-hand  side  define  the  strain  energy  stored  by  vinue  of  the  flexural 
rigidity  E„I„  of  the  beam  and  the  energy  stored  in  the  NITINOL  fibers  due  to  their 

transverse  rfeflerrinn  w  while  under  fensinn  T  resnecTiveiv 
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In  eqn  (1),  the  thermal  loads  P,  is  generated  by  changes  in  the  temperature  A0  of 
the  element  caused  by  changes  in  the  ambient  temperature  or  by  the  activation  and 
de-activation  of  the  NITINOL  fibers.  It  is  given  by 


P,  =  Of  A0 


(2) 


where  a  is  the  thermal  expansion  coefficient  of  the  composite,  E„  is  its  modulus  of 
elasticity,  and  A„  is  the  cross-sectional  area  of  the  beam. 

Equation  (1)  reduces  to 


Ufim  =  £„/„  (d^w/dx^f  dx-P„  (dw/dxf  dx. 


(3) 


where  P„  is  the  net  axial  force  given  by 


P„  =  (P«  +  P,-  T,). 


(4) 


Defining  a  proper  displacement  function  for  the  composite  beam  element,  one  can 
write  the  deflection  w  as 


w  =  [A][Sl 


(5) 


where  the  elements  of  matrix  [A]  are  a  function  of  x  (Fenner,  1975). 

Accordingly,  dw/dx  and  d^w/dx^  can  be  obtained  by  differentiating  eqn  (5)  with 
respect  to  x  to  yield 

dw/dx  =[C][S]  and  d^w/dx^  =  [D][6].  (6) 

If  the  stiffness  matrix  [K^]  of  the  element  is  defined  by  the  following  relationship 

[F]  =  [K,][dl  (7) 

then  [Kf]  can  be  determined  by  combining  eqns  (3),  (6)  and  (7)  as  follows: 

[K,]  =  [DV{D]  dx-P,  [C]^[C]  dx.  (8) 

The  element  stiffness  matrix  [A,],  given  by  eqn  (8),  consists  of  two  components: 
the  conventional  transverse  stiffness  and  the  geometric  stiffness  that  combines  the  effect 
of  the  zocial  mechanical  loads,  axial  thermal  loads  and  the  tension  of  the  NITINOL 
reinforcing  fibers.  Equation  (8)  also  represents  the  basic  equation  for  understanding  the 
role  that  the  NITINOL  fibers  can  play  in  controlling  the  stiffness  of  the  composite  beam. 
For  example,  if  the  initial  fiber  tension  Tq,  resulting  from  the  pre-strain  alone,  is  high 
enough  to  counter-balance  the  mechanical  and  thermal  effects  (i.e.  P„  =  0),  then  the  beam 
stiffness  can  be  maintained  unchanged.  For  higher  pre-strain  levels,  the  beam  stiffness  can 
be  enhanced.  Further  enhancement  can  be  achieved  when  the  shape  memory  effect  of  the 
NITINOL  fibers  is  activated  by  heating  the  fibers  above  their  phase  transformation 
temperature.  The  additional  tension,  induced  into  the  fibers  by  the  phase  recovery  forces, 
makes  the  net  axial  load  P„  negative  and  accordingly  increases  the  overall  stiffness  of  the 
beam.  However,  it  is  essential  that  the  total  tension  in  the  NITINOL  fibers,  i.e.  the 
sum  of  the  tension  due  to  the  pre-strain  and  the  phase  recovery  force,  must  exceed  the 
mechanical  and  thermal  loads  and  compensate  for  the  softening  effect  in  the  matrix 
resulting  from  heating  the  NITINOL  fibers  inside  the  composite  matrix. 


% 


2.2.  The  mass  matrix 

The  element  nifQJ)  of  the  mass  matrix  [M,]  of  the  beam  element  is  obtained  using 
the  consistent  mass  formulation  (Zienkiewicz  and  Taylor,  1989)  as  follows: 

rrteiij)  *  p„A„  AiAjdx  (9) 

where  /I,  and  Aj  are  the  ith  and  yth  elements  of  the  vector  A  given  by  eqn  (5).  Also  in 
eqn  (9),  p„  denotes  the  density  of  the  composite  beam. 
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The  classical  finite  element  approach  is  used  to  form  the  equations  of  motion  of  the 
assembly  of  several  beam  elements  along  with  the  appropriate  boundary  conditions.  The 
solution  of  the  eigenvalues  of  the  resulting  homogeneous  equations  gives  the  natural 
frequencies  of  the  composite  beam  as  influenced  by  the  properties  of  the  matrix  and  the 
NITINOL  fibers.  It  is  important  to  note  that  these  properties  are  influenced  by  the 
temperature  distribution  inside  the  beam  which  is  developed  by  virtue  of  activating  and 
de-activating  the  NITINOL  fibers.  A  study  of  the  temperature  distribution  inside 
NlTINOL-reinforced  beams  is  presented  in  Section  3. 

The  analysis  presented  here  is  for  an  orthotropic  laminate  that  has  a  single  layer  of 
unidirectional  NITINOL  fibers.  Such  an  analysis  can  be  used  along  with  the  classical 
laminate  theory  to  assemble  the  stiffness  matrix  for  a  multi-laminate  composite  beam. 
A  similar  approach  is  carried  out  for  modeling  the  static  and  dynamic  characteristics  of 
NlTINOL-reinforced  composite  plates  (Ro,  1992). 

3.  THERMAL  CHARACTERISTICS  OF  NITINOL-REINFORCED  COMPOSITE  BEAMS 

The  thermal  characteristics  of  NlTINOL-reinforced  composite  beams  are  influenced 
primarily  by  the  temperature  distribution  inside  the  composite.  A  thermal  finite  element 
model  is  developed  to  determine  steady-state  and  transient  temperature  distributions 
resulting  from  different  activation  strategies  of  the  NITINOL  fibers.  The  theoretical 
predictions  are  compared  with  experimental  measurements  in  order  to  validate  the 
thermal  model. 

It  is  important  here  to  note  that  although  the  finite  element  model  used  in  predicting 
the  beam  dynamics  is  a  one-dimensional  model,  with  the  single  dimension  taken  along  the 
beam  neutral  axis,  the  thermal  model  is  considered  to  be  two-dimensional  model  to 
predict  the  temperature  distribution  over  the  beam  cross-section.  Such  a  distinction  is 
made  because  the  temperature  distribution,  along  the  beam  neutral  axis,  is  assumed  to  be 
uniform.  This  assumption  is  confirmed  experimentally  and  is  attributed  to  the  fact  that 
the  NITINOL  fibers  are  oriented  parallel  to  the  neutral  axis  (Baz  et  aL,  1992).  The  beam 
temperature,  however,  varies  only  over  the  cross-section  and  its  magnitude  depends  on  the 
number  and  location  of  the  activated  or  de-activated  NITINOL  fibers.  The  resulting 
temperature  distribution  can  be  used  to  compute  an  average  modulus  of  elasticity  of  the 
composite.  The  average  temperature  rise  above  ambient  can  also  be  used  to  compute  the 
axial  thermal  loading  on  the  beam  P,  which  results  from  fixing  the  two  ends  of  the  beam. 

The  two-dimensional  thermal  modeling  of  the  beam  is  favored  over  a  one-dimen¬ 
sional  lumped-parameter  approach  because  it  provides  a  more  accurate  simulation  of  the 
thermal  state  of  the  beam. 


THE  THERMAL  FINITE  ELEMENT  MODEL 


The  energy  balance  equation  that  governs  the  heat  transfer  across  the  beam  can  be 
written,  in  a  two-dimensional  Cartesian  coordinate  system,  as  follows 


d^Q  d^Q  Q  _ 

dy-  dz~  ^  k  a  dt 


(10) 


where  ©  is  the  beam  temperature  at  time  t  and  location  (y,  c)  as  defined  in  Fig.  3.  In 
eqn  (10),  Q  defines  the  rate  of  heat  generated  per  unit  area  during  the  activation  of 
the  NITINOL  fibers.  Also,  k  denotes  the  conductivity  of  the  beam  and  a  its  thermal 
diffusivity. 

The  above  equation  is  subject  to  the  following  boundary  and  initial  conditions 
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Q  =  0 


k 


f  ^ 

I 


K  - 


a© 

dz 


‘-A 

J 


on  boundary  5, . 
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Fig.  3.  Schematic  drawing  of  the  beam  cross-section  with  heat  transfer  boundaries. 


and 

©(y.  f  =  0)  =  ©o(y,  z)  on  the  beam  cross  section.  (13) 

Equation  (11)  defines  the  condition  at  the  NITINOL  fiber  circular  boundary  Si  on 
which  the  heat  flux  Q  is  specified,  and  eqn  (12)  specifies  the  conditions  at  the  beam  outer 
boundary  S2  where  the  interaction  with  the  ambient  temperature  ©»  is  through  convective 
heat  transfer  with  coefficient  h.  The  boundaries  Si  and  S2  are  defined  in  Fig.  3.  In 
eqns  (11)  and  (12),  ly  and  /j  denote  the  direction  cosines  of  the  outward  normals  to  the 
boundaries.  Equation  (13)  describes  the  initial  temperature  distribution  over  the  beam 
cross-section  at  time  t  =  0. 

Assuming  a  linear  interpolating  function  [N]  with  triangular  elements  that  have 
isotropic  thermal  properties,  then  the  temperature  ©  at  any  y,  2  and  t  can  be  expressed  in 
terms  of  the  nodal  temperatures  [©']  as  follows 


©  =  [N]m. 


(14) 


Using  the  Galerkin  method  along  with  the  assumed  interpolating  functions,  one  can 
write  the  following  finite  element  equation  (Rao,  1988) 


where 


and 


with 


[kUm  +  ([A:|]  +  [Arnxen  =  [P'l, 

[Arf]  ^  I  /:/a[iV]''[yV]  cL4, 

h[Nf[N]dS2, 

hi 

[kt]  =  J^jAr(5n5]cL4. 

[P'l  =  f  f  Q[Ni^dA  -  [  e[N]^d5,  +  he,[NVdS2 
U'  J  Jsf  Jsf 


[B]  = 


dNi/dy  dN^/dy  ...  dNp/dy 

\_dNi/dz  dN^/dz  ...  dNp/dz] 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


where  subscript  p  is  the  number  of  vertices  of  the  element  {p  =  3  for  a  triangular  element). 

The  Individual  element  equations  are  assembled  to  form  the  overall  equation  of  the 
NITINOL-composite  beam  which  can  be  solved  for  the  nodal  temperatures.  The  solution 
is  based  on  a  Crank-Nicolson  trapezoidal  scheme  (Hughes,  1977). 

The  thermal  finite  element  model  developed  in  Section  3  is  used  to  generate  the 
temperature  distribution  over  the  beam  cross-section  which  in  turn  is  utilized  to  compute 
the  average  properties  of  the  beam,  as  for  example  its  modulus  of  elasticity,  under 
different  operating  conditions.  The  theoretical  predictions  of  the  thermal  characteristics 
will  be  compared  with  experimental  measurements,  obtained  in  Section  4,  in  order  to 
check  the  validity  of  the  thermal  model. 
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4. 1 .  Experimental  set-up,  materials  and  methods 

The  characteristics  of  NITINOL-reinforced  beams  are  computed  using  the  developed 
dynamic  and  thermal  models.  The  theoretical  predictions  are  experimentally  validated 
using  a  composite  beam  made  of  randomly  oriented  glass  fibers  embedded  in  a  low  cure 
temperature  polyester  resin.  The  beam  is  30  cm  long,  2.5  cm  wide  and  0.156  cm  thick 
mounted  in  a  clamped-clamped  arrangement.  The  temperature  dependence  of  the 
modulus  of  elasticity  of  the  beam,  shown  in  Fig.  4,  is  obtained  experimentally  using  the 
Dynamic,  Mechanical,  and  Thermal  Analyzer  (DMTA)  of  Polymer  Laboratories  Ltd 
(1990). 

Four  NITINOL  55  fibers,  that  are  0.55  mm  in  diameter,  are  embedded  inside  the 
beam  through  vulcanized  rubber  sleeves  that  have  an  outer  diameter  of  0.95  mm.  The 
NITINOL  fibers  used  have  an  austenite  transformation  temperature  of  about  50° C.  The 
performance  of  the  composite  beam  is  monitored  when  different  subsets  of  the  NITINOL 
fibers  are  intentionally  activated  by  a  controlled  electrical  current.  Monitored  also  is 
the  unintentional  activation  of  the  remaining  NITINOL  fibers  which  results  from  the 
thermal  energy  propagating  from  the  electrically  activated  fibers.  The  effect  of  such  an 
unintentional  activation  on  the  overall  performance  of  the  NITINOL-reinforced  beam  is 
determined,  in  detail,  for  different  intentional  activation  strategies. 

The  experimental  set-up  used  in  measuring  the  interaction  between  the  intentional 
electrical  and  unintentional  thermal  activation  of  the  NITINOL  fibers  is  shown  in  Fig.  5. 
In  the  figure,  the  beam  is  clamped  in  a  fixed-fixed  arrangement  and  the  NITINOL  fibers 
are  divided  into  two  sets  which  are  separately  clamped.  The  first  set  is  electrically 
activated  whereas  the  second  set  is  activated  by  the  unintentional  thermal  energy  propa¬ 
gating  from  the  first  set  through  the  composite  beam.  In  the  arrangement  shown  in 
Fig.  5,  the  first  set  includes  fibers  number  1,  2  and  4  while  the  second  set  is  made  only  of 
fiber  number  3.  The  phase  recovery  forces  developed  by  each  set  of  fibers  are  separately 
measured  using  two  separate  load  cells. 

The  time  history  of  the  phase  recovery  forces  developed  by  each  set  of  fibers  is 
sampled  by  a  computer  when  the  first  set  is  subject  to  a  step  electrical  current.  The  corre¬ 
sponding  time  histories  of  the  temperatures  generated  in  the  NITINOL  fibers  and  the 
composite  beam  are  also  monitored  by  the  computer.  When  steady-state  conditions  are 
attained  the  modes  of  vibration  of  the  beam  are  measured  using  a  modally  tuned  impact 
hammer  (Ewins,  1984). 
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Fig.  Et'tect  of  operating  temperature  on  modulus  of  elasticity  of  the  test  fiberglas. 

composite  beam. 
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Fig.  5.  Experimental  set-up  used  in  monitoring  the  performance  of  NITINOL-reinforced  beams. 

4.2.  Experimental  results 

A.  Electrical  actr-ation  of  three  fibers 
A.l.  Phase  recovery  forces 

Figures  6a  and  6b  show  the  time  response  of  the  phase  recovery  forces  developed  by 
the  intentional  electrical  activation  of  fibers  1,  2  and  4,  and  the  unintentional  thermal 
activation  of  fiber  3,  respectively.  In  the  displayed  characteristics,  all  the  NITINOL  fibers 
are  initially  pre-tensioned  to  26.6  N  per  fiber  and  the  electrically  activated  fibers  are 
subjected  to  a  step  electric  heating  of  5.4  W  per  fiber.  Figure  6a  shows  that  the  phase 
recovery  forces  of  the  electrically  activated  fibers  quickly  rise  to  their  peak  value  of  66  N 
per  fiber.  On  the  other  hand,  Fig.  6b  indicates  that  fiber  number  3,  which  is  not  intention¬ 
ally  activated  by  the  electrical  current,  experiences  a  delay  of  about  100  s  before  its  starts 
developing  its  phase  recovery  forces.  This  time  delay  is  the  time  needed  for  the  thermal 
energy  propagating  from  the  electrically  activated  fibers  to  reach  and  heat  fiber  3  above 
its  phase  transformation  temperature.  Following  this  time  delay,  the  phase  recovery  force 
builds  up  relatively  slowly  to  a  peak  value  of  50.4  N  which  is  a  very  large  force  that  cannot 
be  neglected.  The  slow  rise  to  the  peak  force  is  attributed  primarily  to  the  thermal  capaci¬ 
tance  of  the  composite  beam  and  to  the  heat  lost  to  the  surroundings  across  the  beam 
surfaces.  Hence,  the  resulting  unintentional  thermal  activation  of  the  NITINOL  fiber 
number  3  changes  the  total  tensile  force  T,  of  the  reinforcing  fibers  and  the  total  elastic 
energy  of  the  beam.  These  changes  are  rather  complex  and  mainly  depend  on  the  tempera¬ 
ture  distribution  inside  the  composite  beam.  In  addition,  such  changes  are  significant  and 
must  be  accounted  for  as  they  alter  the  stiffness  and  the  dynamics  of  the  entire  beam,  as 
will  be  shown  in  what  follows. 

A. 2.  Modes  of  vibration 

The  effect  of  different  activation  strategies  on  the  first  threee  modes  of  vibration  of 
the  NITINOL-reinforced  composite  beam  is  shown  in  Figs  7a  and  7b  when  the  initial 
tension  is  26.6  N  per  fiber.  These  results  are  obtained  after  steady-state  temperatures 
and  forces  are  attained.  Comparisons  are  also  given  with  the  modes  of  the  beam  with 
unactivated  fibers.  For  instance,  when  fibers  1,  2  and  4  are  electrically  activated  while 
fiber  3  is  thermally  activated,  the  first  three  modes  of  vibration  are  found  to  be  94.9,  192.5 
and  318.7  Hz,  respectively.  This  is  in  comparison  to  87.9,  184.92  and  315.0  Hz  when  all 
the  fibers  are  unactivated.  The  increase  in  the  three  modes  of  vibration,  of  7.9%,  4.1% 
and  1.15%,  respectively,  is  attributed  to  the  increase  in  the  fiber’s  tension  developed 
by  both  the  intentional  electrical  and  the  unintentional  thermal  activation  of  these  fibers. 
In  order  to  isolate  the  two  effects,  fiber  3  is  replaced  by  another  NITINOL  fiber  which  has 
a  very  high  phase  transformation  temperature  of  100°C.  In  this  manner  the  thermal 
energy  propagating  through  the  beam  will  not  result  in  a  high  enough  temperature  to 
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Fig.  6.  Time  history  of  the  recovery  forces  with  three  activated  NITINOL  fibers. 

(a)  Electrically  activated,  (b)  thermally  activated. 

induce  the  thermal  activation  of  fiber  3.  Accordingly,  the  initial  tension  applied  to 
fiber  3  will  remain  unchanged  at  26.6  N  even  after  the  electrical  activation  of  fibers  1,  2 
and  4.  This  is  not  the  case  with  the  original  low  phase  transformation  temperature  fiber 
where  the  tension  increased  from  26.6  to  50.4  N  following  the  electrical  activation  of 
fibers  1,  2  and  4.  In  this  way  the  effect  of  unintentional  thermal  activation  on  the  beam 
dynamics  can  be  isolated  and  quantified. 

The  modes  of  vibration  measured  with  the  high  phase  transformation  temperature 
fiber  are  92.3,  182.5  and  291.3  Hz.  This  corresponds  to  an  increase  of  5.02'’7o,  -0.53% 
and  -7.5%  as  compared  to  the  modes  of  vibration  with  unactivated  fibers.  Accordingly, 
the  thermal  activation  accounts  for  about  36.4%,  113%  and  752%  of  the  total  increase  in 
the  normal  modes.  Such  percentages  are  very  significant  and  cannot  be  neglected.  For  the 
sake  of  completion,  the  effect  of  aaivating  all  the  four  fibers  electricaUy  on  the  natural 
frequencies  of  the  beam  is  shown  in  Fig.  7b  along  with  the  frequencies  obtained  with  all 
fibers  unactivated. 

A. 3.  Temperature  distribution 

In  order  to  develop  a  thorough  understanding  of  the  effect  of  unintentional  thermal 
activation,  it  is  necessary  to  consider  the  temperature  distribution  over  the  cross-seaion 
of  the  composite  beam  shown  in  Fig.  8a.  Figure  8b  shows  the  theoretical  temperature 
distribution  when  fibers  1,  2  and  4  are  electrically  activated,  and  Fig.  8c  shows  the 
corresponding  distribution  when  all  the  four  fibers  are  electrically  energized.  The 
distributions  displayed  represent  the  steady-state  distributions  obtained  by  the  thermal 
finite  element  model  after  700  s  from  the  initiation  of  the  step  heating  of  the  NITINOL 
fibers.  It  is  evident  from  Figs  8b  and  8d  that  the  electrical  activation  of  fibers  1,  2  and  4 
generates  temperatures  around  fiber  3  in  excess  of  its  phase  transformation  temperature. 
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Fig.  7.  Natural  frequencies  of  NITINOL-reinforced  beam,  (a)  Three  activated  fibers, 

(b)  four  activated  fibers. 

Under  such  conditions,  the  phase  transformation  is  induced  thermally  and  fiber  3  starts 
developing  its  phase  recovery  force.  The  extent  of  completion  of  the  phase  transforma¬ 
tion  process  and  the  magnitude  of  the  force  developed  depend  on  how  high  the  beam 
temperature  becomes  in  relation  to  the  phase  transformation  temperature. 

Figure  8d  shows  also  that  the  rise  of  the  beam  temperature,  near  fiber  3,  is  rather 
slow  because  of  the  thermal  capacitance  of  the  beam  and  also  because  of  the  heat  losses 
to  the  environment  across  its  surfaces.  Accordingly,  there  is  a  time  period  before  the 
beam  temperature  rises  to  the  phase  transformation  temperature  and  before  fiber  3 
starts  generating  its  phase  recovery  force.  During  that  period  of  time  and  until  steady- 
state  conditions  are  attained,  the  dynamics  of  the  beam  will  be  continuously  varying. 
Predictions  of  such  behavior  are  only  possible  through  the  interaction  between  the  thermal 
and  dynamic  finite  element  models. 

Comparison  between  the  theoretical  and  the  experimental  temperatures  of  the 
activated  NITINOL  fiber  and  the  beam  surface  are  shown  in  Figs  8d  and  8e  when  three 
and  four  fibers  are  electrically  activated,  respectively.  It  is  evident  that  theoretical 
predictions  are  in  close  agreement  with  the  experimental  results. 

B.  Electrical  activation  of  two  fibers 
B.l.  Phase  recovery  forces 

The  effect  of  electrically  activating  two  NITINOL  fibers  (1  and  2,  1  and  4,  or  2 
and  4)  on  the  phase  recovery  forces  developed  by  the  thermal  activation  of  the  remaining 
two  fibers  is  shown  in  Fig.  9.  Figure  9a  shows  the  forces  developed  by  the  electrically 
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Fig.  8.  Temperature  distribution  over  the  beam  cross-section  with  three  and  four  activated  fibers. 
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Fig.  9.  Time  history  of  the  recovery  forces  with  two  aaivated  NITINOL  fibers, 
(a)  Electrically  activated,  (b)  thermally  activated. 


activated  fibers  and  Fig.  9b  shows  the  corresponding  forces  developed  by  the  remaining 
thermally  activated  fibers.  The  magnitude  of  the  forces  developed  by  the  electrically 
activated  fibers  is  exactly  the  same,  irrespective  of  their  location.  However,  the  magnitude 
of  forces  generated  by  the  thermally  activated  fibers  depends  on  their  location,  inside 
the  beam,  relative  to  the  location  of  the  electrically  activated  fibers.  The  least  forces 
are  developed  when  the  electrically  activated  fibers  are  located  near  the  edges  of  the 
beam  (fibers  1  and  4)  where  most  of  their  generated  thermal  energy  will  be  lost  to  the 
surrounding.  The  remaining  portion  of  the  thermal  energy  will  propagate  inside  the  beam 
to  thermally  activate  the  remaining  two  fibers  (2  and  4).  Higher  forces  are  developed 
when  the  electrically  activated  fibers  (1  and  2,  or  2  and  4)  are  located  away  from  the 
beam  edges.  For  instance,  when  fibers  2  and  4  are  electrically  activated  they  produce  a 
total  force  of  about  132  N  while  the  thermally  activated  fibers  (1  and  3)  develop  a  total 
force  of  about  95  N.  Definitely,  the  magnitude  of  such  onintentionally  generated  forces 
cannot  be  neglected. 

It  is  again  evident  that  the  interaction  between  the  electrical  and  thermal  activations 
plays  a  very  important  role  in  determining  the  total  tension  developed  by  the  NITINOL 
reinforcement  which  in  turn  determines  the  dynamic  behavior  of  the  entire  composite 
beam. 


B.2.  Modes  of  vibration 

The  effect  of  diferent  electrical  activation  strategies  of  two  fibers  on  the  modes  of 
vibration  of  the  beam  is  shown  in  Fig.  10  and  the  results  are  summarized  in  Table  1. 
These  results  are  obtained  after  steady-state  temperatures  and  forces  are  attained. 
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Fig.  10.  Natural  frequencies  of  a  NITINOL-reinforced  beam  with  two  activated  fibers. 

(a)  Fibers  1  and  2,  (b)  flbers  1  and  4,  and  (c)  fibers  2  and  4. 

Table  I  and  Fig.  10  indicate  that  the  dynamic  performance  of  the  composite  beam, 
is  significantly  influenced  by  the  different  electrical  activation  strategies.  This  is  in  spite 
of  the  fact  that  two  fibers  are  always  activated  in  all  these  strategies.  The  studies  of  Rogers 
et  al.  (1991)  and  Jia  and  Rogers  (1989)  predict  that  the  beam  performance  remains  the 
same  for  all  the  cases  considered. 


Table  1.  Effect  of  activation  strategy  of  two  fibers  on  the 
modes  of  vibration  of  the  beam 


Number  of 
fibers 
aaivated 

Modes  of  vibration  (Hz) 

First 

Second 

Third 

None 

87.9 

184.9 

315.0 

All 

100.1 

202.5 

323.7 

1  and  4 

88.4 

182.5 

305.0 

1  and  2 

88.8 

183.8 

307.4 

2  and  4 

93.5 

190.1 

315.0 
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Furthermore,  it  is  interesting  to  note  that  when  the  effect  of  the  unintentional 
thermal  activation  is  negligible,  as  when  fibers  1  and  4  are  electrically  activated,  the 
beam  becomes  less  stiff  than  a  beam  with  unactivated  fibers.  This  is  clearly  manifested 
by  the  drop  in  the  frequencies  of  the  second  and  third  modes  of  vibration.  Such  a 
drop  is  attributed  to  the  fact  that  the  phase  recovery  forces  generated  by  fibers  1  and 
4  are  not  high  enough  to  compensate  for  the  softening  effect  of  the  matrix  modulus 
due  to  the  heating  of  the  fibers.  However,  when  the  unintentional  thermal  activation 
effect  becomes  significant,  as  when  fibers  2  and  4  are  electrically  activated,  the  beam 
performance  is  considerably  enhanced.  This  is  evident  from  the  significant  increase  in  the 
natural  frequencies  of  the  beam.  The  improved  performance  is  due  to  the  additional 
recovery  forces  developed  by  the  unintentional  thermal  activation  of  fibers  1  and  3. 

B .  3 .  Temperature  distribution 

The  interaction  between  the  thermal  and  dynamic  behavior  of  the  beam  can  best 
be  understood  by  considering  the  temperature  distributions  across  the  cross-section  of 
the  NITINOL-reinforced  beam  for  different  activation  strategies.  Figure  11  displays 
such  temperature  distributions  along  with  comparisons  between  the  theoretical  and 
experimental  temperatures  of  the  activated  NITINOL  fibers  and  the  beam  surface.  Close 
agreement  between  the  theoretical  and  experimental  results  is  evident. 

4.3.  Further  comparisons  between  theory  and  experiments 

Figure  12  shows  comparison  between  the  theoretical  and  the  experimental  natural 
frequencies  of  the  beam  for  different  activation  strategies  and  different  initial  fiber 
tensions.  The  effect  of  interaction  between  the  electrical  and  thermal  activation  of  the 
NITINOL  fibers  is  taken  into  consideration  in  these  results.  Close  agreement  is  also 
evident  between  the  theoretical  predictions  and  experimental  measurements. 


~=£C££n:v  ■  hc 

Fig.  12.  Comparisons  between  theoretical  and  experimental  natural  frequencies  of  the  beam 
(a)  Different  activation,  fbi  different  initial  tension 
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The  dynamic  and  thermal  characteristics  of  NITINOL-reinforced  composite  beams 
have  been  presented.  The  fundamental  issues  governing  the  behavior  of  this  new  class  of 
SMART  composites  have  been  introduced.  Particular  emphasis  has  been  placed  on  the 
interaction  between  the  intentional  electrical  activation  of  subsets  of  the  NITINOL  fibers, 
by  an  active  controller,  and  the  associated  unintentional  activation  of  neighboring 
NITINOL  fibers.  It  is  shown  that  such  an  interaction  plays  a  significant  role  in  deter¬ 
mining  the  dynamic  behavior  of  the  entire  composite  beam.  Furthermore,  such  interac¬ 
tion  cannot  be  neglected  as  in  the  previous  studies  of  Rogers  et  al.  (1991)  and  Jia  and 
Rogers  (1989). 

It  is  also  shown  that  the  dynamic  and  thermal  models  developed  in  this  study  enable 
the  accurate  prediction  of  the  behavior  of  NITINOL-reinforced  composite  beams.  With 
such  models  it  would  be  possible  to  design  NITINOL-reinforced  composites  that  have 
continuously  tunable  structural  characteristics  to  adapt  to  changes  in  the  operating 
conditions.  These  features  will  be  particularly  useful  in  many  critical  structures  that  are 
intended  to  operate  autonomoasly  for  long  durations  in  isolated  environments  such  as 
defense  vehicles,  space  structures  and  satellites. 
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1.  INTRODUCTION 


Considerable  attention  has  been  devoted  recently  to  the  utilization 
of  the  Shape  Memory  Nickel-Titanium  alloy  (NITINOL)  in  developing  SMART 
composites  that  are  capable  of  adapting  intelligently  to  external 
distvirbances  (Ikegami  et  al.  1990,  Rogers  et  al.l991,  and  Baz  et  al. 
1990  and  1991).  Such  wide  acceptance  of  NITINOL  stems  from  its  unique 
behavior  when  it  is  subjected  to  particular  heating  and  cooling 
strategies.  For  example,  the  alloy  becomes  soft  when  it  is  cooled 
below  its  martensite  transformation  temperature  and  becomes  about  four 
times  stiffer  when  it  is  heated  above  its  austenite  transformation 
temperature  (Funakubo  1987).  Furthermore,  if  the  alloy  is  trained  to 
have  a  particular  shape  while  in  its  austenite  phase,  it  will  memorize 
this  shape.  If  the  alioy  is  then  cooled  to  its  martensite  phase  and 
subject  to  plastic  deformation,  it  wiil  return  to  its  memorized  shape 
when  it  is  heated  above  the  austenite  transformation  temperature.  The 
phase  transformation  from  martensite  to  austenite  produces  significant 
forces  as  the  alloy  recovers  its  original  shape.  The  alloy  acts  as  an 
actuator  transforming  thermal  energy  into  mechanical  energy  (Perkins 
1975  and  Duerig  et  al.  1990).  Accordingly,  if  the  NITINOL  fibers  are 
embedded  inside  a  composite  matrix  at  optimal  locations,  they  can  be 
used  to  control  the  static  and  dynamic  characteristics  of  the  resulting 
SMART  composite.  The  control  action  is  generated  by  the  described 
stiffening  of  the  NITINOL  fibers  and/or  the  shape  memory  effect.  With 
such  built-in  control  capabilities,  the  performance  of  the  SMART 
composites  can  be  optimized  and  tailored  to  match  changes  in  operating 
conditions. 

EImphasis  is  placed,  in  the  present  work,  on  using  the  shape  memory 
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effect  of  the  NITINOL  fibers  to  control  the  performance  of  fiberglass 
composite  beams.  The  NITINOL  fibers  are  embedded  inside  vulcanized 
rubber  sleeves  placed  along  the  neutral  axes  of  these  composite  beams 
as  shown  in  Figure  (1).  In  this  arrangement,  the  fibers  are  free  to 
move  during  the  phase  transformation  process  in  order  to  avoid 
degradation  and/or  destruction  of  the  shape  memory  effect  which  may 
result  when  the  fibers  are  completely  bonded  inside  the  composite 
matrix. 


Figure  (1)  -  A  schematic  drawing  of  the  cross  section  of  a 

NITINOL-reinforced  composite  beam 


The  basic  phenomena  governing  the  thermo-dynamic  performance  of  the 
NITINOL  fibers  and  NITINOL-reinforced  composites  will  be  presented. 
The  NITINOL  fibers  will  be  utilized  to  control  the  buckling  and  the 
flow-induced  vibrations  of  NITINOL-reinforced  fiberglass  composite 
beams. 

2.  CHARACTERISTICS  OF  THE  NITINOL  FIBERS 

Knowledge  of  the  thermal  and  dynamic  behavior  of  the  shape  memory 
NITINOL  fibers  is  essential  to  the  understanding  of  their  role  in 
controlling  the  performance  of  NITINOL-reinforced  composites.  The 
thermo-dynamic  behavior  of  the  NITINOL  fibers  has  been  extensively 
studied  throughout  the  last  two  decades  (Funakubo  1987,  Perkins  1975, 
Jackson  et  al.  1972).  However,  we  will  present  a  different  outlook 
which  will  be  crucial  in  developing  the  basic  principles  governing  the 
performance  of  NITINOL-reinforced  composites. 

Emphasis  is  placed"  on  studying  the  effect  of  the  operating 
temperature  and  the  pre-strain  level  on  the  recovery  forces  and,  most 
importantly,  on  the  natural  frequencies  of  end-restrained  NITINOL 
fibers.  Such  end-restrained  fibers  constitute  the  basic  building 
block  of  NITINOL-reinforced  composites. 

Figure  (2)  shows  a  schematic  drawing  of  the  experimental  set-up 
used  to  determine  such  thermo-dynamic  characteristics.  In  the  set-up, 
the  NITINOL  fiber  is  clamped  in  a  holder  at  one  end  and  connected  to  a 
load  cell  at  the  other  end.  The  load  cell  monitors  the  pre-strain 
level  of  the  fiber  when  it  is  in  its  martensitic  phase.  It  also 
provides  continuous  measurements  of  the  recovery  force  when  the  fiber 
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undergoes  its  phase  transformation  due  to  external  electric  heating. 
The  fiber  temperature  is  monitored  by  a  thermocouple  bonded  to  the 
NITINOL  fiber.  The  fiber  assembly  is  mounted  on  an  sliding  table  which 
is  connected  to  a  mechanical  shaker.  Random  excitations  are  used  to 
drive  the  shaker  and  the  table,  thus  applying  a  transverse  displacement 
to  the  fixed  ends  of  the  NITINOL  fiber.  The  resulting  oscillations  of 
the  fiber  are  measured  by  a  non-contacting  magnetic  sensor  mounted  on 
the  table.  The  shape  memory  effect  of  the  NITINOL  fiber  is  energized 
electrically  and  the  resulting  fiber  temperature,  recovery  force  and 
amplitude  of  oscillation  are  continuously  sampled  by  a  digital 
computer. 


■  Power  supply 
To  Analyzer 


Figure  (2)  -  Experimental  set-up  for  monitoring  thermo-dynamic 
behavior  of  NITINOL  fibers. 

The  effect  of  the  pre-strain  level  on  the  recovery  force,  as  a 
function  of  time  during  a  heating  and  cooling  cycle,  and 
force-temperature  characteristics  are  shown  in  Figures  (3-a),  and  (3-b) 
respectively.  The  recovery  force  increases  almost  linearly  with 
increasing  pre-strain  levels.  Such  characteristics  conform  with 
published  results. 

The  new  outlook  on  the  thermo-dynamic  characteristics  of  the  NITINOL 
fibers  is  demonstrated  by  the  effect  of  the  pre-strain  on  the  natural 
frequencies  of  the  fiber  as  shown  in  Figure  (4).  The  figure  displays 
the  spectrum  analysis  of  the  amplitude  of  oscillation  of  the  fiber,  at 
different  pre-strain  levels,  with  and  without  the  activation  of  the 
shape  memory  effect. 

Figure  (4)  indicates  that  activating  the  shape  memory  effect 
results  in  a  significant  increase  of  the  natural  frequency  of  the  fiber 
which  becomes  more  pronounced  with  increasing  pre-strain  levels. 
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Figure  (3)  -  Effect  of  pre-strain  on  recovery  force,  and 

force-temperature  characteristics  of  a  30  cm  iong  22 
mil  NITINOL  fiber. 


Figure  (4)  -  Effect  of  pre-strain  on  the  natural  frequencies  of 
unactivated  (u)  and  activated  (a)  NITINOL  fibers. 

A  better  understanding  of  the  underlying  phenomena  can  be  obtained 
by  treating  the  NITINOL  fiber  using  the  classical  theory  of  vibrating 
strings.  The  wave  equation  for  transverse  vibrations  of  a  undamped 
freely  vibrating  string,  which  is  tightly  stretched  with  a  tension  T^, 
is  given  by  (James  et  al.  1989) 


a^w  j  ^w 


a  X 


c"  a  t^ 


(1) 


where  w  Is  the  transverse  deflection  of  the  string  at  a  distance  x 
along  the  string  at  time  t.  In  equation  (1),  the  constant  c  is  the 
wave  propagation  speed  given  by 

c  =  y  /  Pl  .  (2) 

where  Pj^  is  the  mass  per  unit  length  of  the  string. 


Using  conventional  separation  of  variables  approach  for  a  string  of 
length  L  which  is  fixed  at  both  ends,  the  transverse  deflection  w  is 
written  as 


w  =  J]  <In  .  (3) 

n=l 


where  and  q„  are  the  mode  shape  and  the  generalized  coordinate  for 
the  nth  mode  of  vibration  respectively.  Using  equation  (3),  equation 
( 1 )  reduces  to 


q„  +  (.)„  q^  =  0  ,  for  n=l,  2,  ,  .  (4) 

where  is  the  natural  frequency  of  a  fixed-fixed  string  given  by 


<j„  =  (  n  Tt  /  L  )  y  /  Pl  for  n=l,  2,  .  .  (5) 

Accordingly,  the  natural  frequency  of  the  string  is  proportional  to 
the  square  root  of  its  tension.  Using  this  relationship,  the  effect  of 
the  pre-strain  (or  the  initial  tension)  of  the  NITINOL  fibers  on  the 
first  natural  frequency  of  the  unactivated  and  activated  fiber  is  shown 
in  Figure  (5).  Two  distinct  linear  characteristics  are  observed  with  a 
significant  increase  in  the  natural  frequency  when  the  fiber  is 
activated.  A  unified  characteristic  caui  be  obtained  when  the  effect 
of  total  tension  T^,  which  is  the  sum  of  the  initial  tension  and  the 
phase  transformation  force  of  the  activated  fibers,  is  considered  as 
shown  in  Figure  (6).  The  natural  frequencies  of  the  unactivated  and 
activated  fiber  fall  on  a  single  straight  line  which  has  a  slope 
of(n/2L  v'pl  ) .  In  this  analysis  the  effect  of  thermal  expansion  on  the 
fiber  tension  is  negligible  as  compared  to  the  phase  transformation 
force. 

Therefore,  the  classical  theory  of  vibrating  strings  can  be  used  to 
predict  the  dynamics  of  unactivated  as  well  as  activated  NITINOL 
fibers.  Accordingly,  the  theory  of  vibrating  strings  can  be  utilized 
to  determine  the  strain  energy  stored  in  NITINOL  fibers  embedded  inside 
composite  beams  as  the  beams  deflect  from  their  equilibrium  position 
under  the  action  of  external  loads.  Using  this  approach  to  determine 
the  thermo -dynamic  behavior  of  NITINOL  fibers,  it  is  possible  to 
develop  a  thorough  xinder standing  of  the  static  and  dynamic  performance 
of  NITINOL-reinforced  composites. 
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Figure  (5)  -  Effect  of  pre-strain  (initial  tension)  on  the  first 
natural  frequency  of  unactivated  and  activated  NITINOL 
fiber. 


Total  Tension^/2 


Figure  (6)  -  Effect  of  total  tension  on  the  first  natural  frequency  of 
unactivated  and  activated  NITINOL  fiber. 
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3.  CHARACTERISTICS  OF  NITINOL-REINFORCED  COMPOSITE  BEAMS 

3.1.  STATIC  CHARACTERISTICS 


The  static  characteristics  of  NITINOL-reinforced  composite  besuns 
are  primarily  governed  by  their  stiffness.  The  beam  stiffness  is  made 
up  of  different  components  which  Include:  the  flexural  rigidity  of  the 
beam,  the  geometric  stiffness  that  accounts  for  the  axial  and  thermal 
loading  as  well  as  the  stiffness  imparted  by  the  elasticity  of  the 
NITINOL  fibers.  The  individual  components  of  the  beam  stiffness  can  be 
determined  by  considering  the  NITINOL-reinforced  beam  element  shown  in 
Figure  (7)  with  the  forces  acting  on  it  amd  the  associated 
displacements.  The  combined  stiffness  of  the  element  caui  be  obtained 
using  the  principle  of  conservation  of  energy  and  equating  the  work 
done  by  external  loads  to  the  strain  energies  stored  in  the  element.  In 
the  present  analysis,  the  theory  of  Bernoulli -Euler  beams  is  used  with 
the  assumption  of  small  deflections. 


Figure  (7)  -  NITINOL-reinforced  beam  element  with  forces  and 

resulting  displacements 

3.1,1.  EXTERNAL  WORK 

The  work  done  by  the  external  loads  includes: 
a.  work  done  by  transverse  loads  and  moments  (V^) 

This  work  is  given  by 

Wi  =  1/2  [dp  [FI,  (6) 

where  [5]  and  [F]  are  the  displacement  and  transverse  loads  vectors, 
respectively,  given  by 


(7) 


[3]  =  [Wj  Wj 

and  [F]  =  [V^  M,  Vj  Mj]^  (8) 

with  Wj  and  dj  are  the  linear  and  angular  deflections  of  node  i. 
respectively  and  and  M|  are  the  shear  and  moment  acting  at  node  i, 
respectively. 

b.  work  done  by  the  axial  nechanical  loads  (V^) 


Wo-  is  given  by 


W2,  =  P,  /  2  J  (dw/dx)2  dx. 


where  P,  is  the  external  axial  compressive  load  acting  along  the 
neutral  axis  of  the  beam  element. 

c.  work  done  by  the  axial  thermal  loads  (V2t) 

W2t  represents  the  work  done  by  the  thermal  loads  P^  on  the  beam 
element  due  to  changes  in  the  temperature  A0  of  the  element  caused  by 
changes  in  the  ambient  temperature  or  during  the  activation  and 
de-activation  of  the  NITINOL  fibers.  It  is  given  by 


Wat  “  ^  J  (dw/dx)^  dx. 


where  is  given  by 

Pt  =  a  i8  E„  A„,  (11) 

where  a  is  the  thermal  expansion  coefficient  of  the  composite,  EI^  is 
its  modulus  of  elasticity  and  A„  is  the  beam  cross  sectional  area. 

3.1.2.  STORED  STRAIN  ENERGY 

The  stored  strain  energy  consists  of  two  components: 
a.  strain  energy  of  beam  (V3} 

The  energy  stored  in  the  beam  element  due  to  its  bending  is  given 


W3  =  E„  /  2  J  (d^#/dx2)2  dx 


where  E„  is  the  flexural  rigidity  of  the  beam. 


b.  strain  energy  of  NITINOL  fibers 

Considering  the  NITINOL  fiber  as  a  string  with  a  tension  T  which  is 
displaced  laterally  a  distance  w  from  the  neutral  axis  of  the  bean. 
Then  its  stored  strain  energy  is  given  by 

L 

W4  =  T  /  2  I  (dw/dx)2  dx.  (13) 

0 

Equating  the  sum  of  the  work  done  by  the  external  forces  F,  and 
Pt  to  the  sum  of  the  strain  energies  stored  in  beam  and  NITINOL  fibers 
gives 


*  ^2^  *  Wat  =  W3  +  W4.  (14) 

Substituting  equations  (6).  (9),  (10),  (12)  and  (13)  into  equation 
(14)  yields 


L  L 

=  E.  I,  I  (d2w/dx2)2  dx  -  P„  I  (dw/dx)2  dx,  (15) 

0  0 


where  P„  is  the  net  axial  force  give  by 

Pn  =  (  P-  +  Pt  -  T).  (16) 

Defining  a  cubic  displacement  function  for  the  composite  beam 
element,  of  the  following  form 

w=a+bx+cx^+dx^  (17) 

where  a,  b,  c  and  d  are  constants  that  can  be  calculated  in  terms  of 
the  deflections  of  the  nodes  i  and  j  bounding  the  beam  element.  Then 
equation  (17)  can  be  rewritten  as 

w  =  (AI  (5],  (18) 

where  the  elements  of  matrix  [Al  are  function  of  x  (Fenner  1975). 

Accordingly,  dw/dx  and  d^/dx^  can  be  obtained  by  differentiating 
equation  (18)  with  respect  to  x  to  yield 

dw/dx  =  (Cl  (31  and  (i^/dx^=  [D]  [5],  (19) 

where  the  matrices  (Cl  and  (Dj  are  given  by 


and 


[D]  = 


(20) 


[C]  = 


dx 


((AJ) 


d2 

55^2 


([A]). 


The  following  relationships  can  also  be  obtained  from  ecjuation  (19) 

(dw/dx)2  =[5]‘r[C]'^rc)(5]  and  (d^w/dx^)^  »[5]T[d]T[d]  [5] .  (21) 

If  the  stiffness  matrix  (K^l  of  the  element  is  defined  by  the 
following  relationship 

[F]  =  [Kel  [51.  (22) 

then,  [Kg]  can  be  determined  by  combining  equations  (15),  (21)  and  (22) 
as  follows 


L  L 

[Kel  =  I«  J  [D]T[r]  dx  -  P„  I  [C]’’[C]  dx  (23) 
0  0 

It  can  be  seen  from  equation  (23),  that  the  element  stiffness 
matrix  [Kg]  consists  of  two  components:  the  conventional  transverse 
stiffness  and  the  geometric  stiffness  that  combines  the  effect  of  the 
axial  mechanical  loads,  axial  thermal  loads  and  the  tension  of  the 
NITINOL  reinforcing  fibers.  Equation  (23)  also  represents  the  basic 
equation  for  understanding  the  role  that  the  NITINOL  fibers  can  play  in 
controlling  the  static  characteristics  of  the  composite  beam.  For 
example,  if  the  beam  is  not  reinforced  by  NITINOL  fibers  (i.e.  T  =  0) 
and  the  mechanical  and  thermal  loads  induce  compressive  stresses  in  the 
beam,  then  the  geometric  stiffness  will  increase  and  the  total  element 
stiffness  will  decrease.  When  the  combined  effect  of  the  mechanical 
and  thermal  loads  reaches  a  critical  magnitude  such  that  the  geometric 
stiffness  becomes  equal  to  the  flexural  stiffness  of  the  beam,  the 
beam  stiffness  vanishes  and  the  beam  becomes  elastically  unstable. 
Subjecting  the  beam  to  any  additional  external  disturbance  will  cause 
the  beam  to  buckle. 

It  should  be  pointed  out  that  the  thermal  loading,  as  it  increases 
the  geometric  stiffness,  also  decreases  the  flexural  stiffness  of  the 
beam  because  it  reduces  its  effective  Young’s  modulus  Such  a  dual 
effect  makes  the  beam  buckle  under  smaller  thermal  loads  than  under 
pure  mechanical  loading. 

However,  the  critical  load  of  the  un-reinforced  beam  can  be 
increased  by  embedding  pre-strained  NITINOL  fibers  into  the  beam.  If 
the  tension  T,  resulting  from  the  pre-strain  alone,  is  high  enough  to 
counter-balance  the  mechanical  and  thermal  effects  then  the  beam 
stiffness  can  be  maintained  unchanged.  For  higher  pre-strain  levels, 
the  beam  stiffness  can  be  enhanced.  Further  enhancement  can  be 
achieved  when  the  shape  memory  effect  of  the  NITINOL  fibers  is 
activated  by  heating  the  fibers  above  their  phase  transformation 
temperature.  The  additional  tension,  induced  into  the  fibers  by  the 
phase  recovery  forces,  makes  the  net  axial  load  P„  negative  and 


increases  accordingly  the  overall  stiffness  of  the  beam  element. 
However,  it  is  essential  that  the  total  tension  in  the  NITINOL  fibers, 
i.e. ,  the  sum  of  the  tension  due  to  the  pre-strain  and  the  phase 
recovery  force,  must  exceed  the  mechauiical  and  thermal  loads  and 
compensate  for  the  softening  effect  in  the  matrix  resulting  from 
heating  the  NITINOL  fibers  inside  the  composite  matrix. 

Therefore,  effective  control  of  the  stiffness  of  NITINOL-reinforced 
composites  can  be  achieved  by  proper  selection  of  the  initial 
pre-strain  level  of  the  NITINOL  fibers.  This  selection  is  particularly 
crucial,  in  view  of  the  results  of  Figure  (3),  as  the  pre-strain  level 
determines  the  generated  levels  of  recovery  forces. 

The  finite  element  model  of  the  NITINOL-reinforced  beams  describes 
the  interaction  between  the  external  loads,  operating  conditions  and 
the  geometrical  and  physical  parameters  of  the  composite  beam  and  the 
NITINOL  fibers.  It  defines  how  the  NITINOL  fibers  can  be  utilized  to 
tailor  the  stiffness  of  the  composite  to  compensate  for  environmental 
and  operating  conditions  and  disturbances.  The  stiffness  obtained  for 
the  individual  elements  of  the  beam  can  be  assembled  using  the 
classical  finite  element  approach  (Fenner  1975).  The  assembled  model 
can  then  be  subjected  to  the  appropriate  boundary  conditions  in  order 
to  compute  the  deflections  corresponding  to  particular  external  loading 
conditions.  The  analysis  presented  is  for  an  orthotropic  laminate  that 
has  a  single  layer  of  unidirectional  NITINOL  fibers.  Such  an  analysis 
can  be  used  along  with  the  classical  laminate  theory  to  assemble  the 
stiffness  matrix  for  a  multi-laminate  composite  beam.  Similar  approach 
can  be  carried  out  for  modeling  the  static  and  dynamic  characteristics 
of  NITINOL-reinforced  composite  plates. 

The  finite  element  model  developed  will  be  validated  with 
experimental  results  obtained  with  fiberglass  composite  beams. 


3.2.  DYNAMIC  CHARACTERISTICS 


The  dynamic  characteristics  of  NITINOL-reinforced  beams  are 
obtained  by  combining  the  stiffness  matrix  [K^]  with  the  mass  matrix 
[Mg]  of  the  beam  to  form  the  following  element  equation  of  motion 

[Mg]  [5]  +  [Kg]  [5]  =  [F],  (24) 

where  [5]  is  the  nodal  acceleration  vector.  The  elements  mg(i,j)  of  the 
element  mass  matrix  [Mg]  are  obtained  using  the  consistent  mass 
formulation  (Zienkiewicz  and  Taylor  1989)  as  follows 


where  [Aj]  and  [Aj]  are  the  1th  and  Jth  elements  of  the  vector  A  given 
by  equation  (18). 

The  classical  finite  element  approach  is  used  to  form  the  equations 
of  motion  of  the  assembly  of  several  beam  elements  and  the  appropriate 
boundary  conditions  are  then  applied.  The  solution  of  the  eigenvalues 
of  the  resulting  homogeneous  equations  give  the  natural  frequencies  of 
the  composite  beam  as  influenced  by  the  properties  of  the  matrix  and 
the  NITINOL  fibers.  It  is  important  to  note  that  these  properties  are 
influenced  by  the  temperature  distribution  inside  the  beam  which  is 
developed  by  virtue  of  activating  and  de-activating  the  NITINOL  fibers. 
A  study  of  the  temperature  distribution  Inside  NITINOL-reinforced  beams 
will  be  presented  in  the  following  section. 


3.3.  THERMAL  CHARACTERISTICS 


The  thermal  characteristics  of  NITINOL-reinforced  composite  beams 
are  influenced  primarily  by  the  temperature  distribution  inside  the 
composite.  A  thermal  finite  element  model  is  developed  to  determine 
steady-state  and  transient  temperature  distribution  resulting  from 
different  activation  strategies  of  the  NITINOL  fibers.  The  theoretical 
predictions  are  compared  with  experimental  measurements  in  order  to 
validate  the  thermal  model. 

It  is  important  here  to  note  that  although  the  finite  element  model 
used  in  predicting  the  beam  dynamics  is  a  one-dimensional  model,  with 
the  single  dimension  taken  along  the  beam  neutral  axis,  the  thermal 
model  is  considered  to  be  a  two-dimensional  model  to  predict  the 
temperature  distribution  over  the  beam  cross  section.  Such  a 
distinction  is  made  because  the  temperature  distribution,  along  the 
beam  neutral  axis,  is  assumed  to  be  uniform.  This  assumption  is 
confirmed  experimentally  and  is  attributed  to  the  fact  that  the  NITINOL 
fibers  are  oriented  parallel  to  the  neutral  axis.  The  beam 
temperature,  however,  varies  only  over  the  cross  section  and  its 
magnitude  depends  on  the  number  and  location  of  the  activated  or 
de-activated  NITINOL  fibers.  The  resulting  temperature  distribution 
can  be  used  to  compute  an  average  Young’s  modulus  of  the  composite. 
The  average  temperature  rise  above  the  ambient  can  also  be  used  to 
compute  the  axial  thermal  loading  on  the  beam  which  results  from 
fixing  the  two  ends  of  the  beam. 

The  two-dimensional  thermal  modeling  of  the  beam  is  favored  over  a 
one-dimensional  lumped-parameter  approach  because  it  provides  more 
accurate  simulation  of  the  thermal  state  of  the  beam. 

3.3.1  THERMAL  FINITE  ELEMENT  MODEL 

The  energy  balance  equation  that  governs  the  heat  transfer  across 
the  beam  can  be  written,  in  two-dimensional  cartesian  coordinate  system 
,  as  follows 
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a^e  32q  q  13  0 

-  +  -  +  -  =  -  -  (26) 

8y2  az^  k  aat 

where  6  Is  the  beam  temperature  at  time  t  and  location  (y.  z)  as  defined 
in  Figure  (8).  In  equation  (26),  Q  defines  the  rate  of  heat  generated 
per  unit  area  during  the  activation  of  the  NITINOL  fibers.  Also,  k 
denotes  the  conductivity  of  the  beam  and  a  its  thermal  dlffuslvity. 

The  above  equation  is  subject  to  the  following  boundary  and 
initial  conditions 

P  a  e  a  e  . 

k  -  ly  +  -  £*  +  Q  =  0  on  boundary  Sj,  (27) 

Lay  3  z  J 

r  a  0  a  0 

k  -  ty  +  -  +  h  (0  -  0^)  =0  on  boundary  S2,  (28) 

Lay  a  z  J 


and  0  (y,  z,  t=0)  =  0^  (y,  z)  on  beam  cross  section  (29) 

Equation  (27)  defines  the  condition  at  the  NITINOL  fiber  circular 
boundary  on  which  the  heat  flux  Q  is  specified,  and  equation  (28) 
specifies  the  conditions  at  the  beam  outer  boundary  Sg  where  the 
interaction  with  the  ambient  temperature  0,  is  through  convective  heat 
transfer  with  coefficient  h.  The  boundaries  Sj  and  Sg  are  defined  in 
Figure  (8),  In  equations  (27)  and  (28),  ty  and  t^  denote  the  direction 
cosines  of  the  outward  normals  to  the  boundaries.  Equation  (29) 
describes  the  initial  temperature  distribution  over  the  beam  cross 
section  at  time  t  =  0. 

y 


Figure  (8)  -  Schematic  drawing  of  beam  cross  section  with  heat  transfer 
boundaries. 
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Assxuning  a  linear  interpolating  function  [Nl  with  triangular 
elements  that  have  Isotropic  thermal  properties,  then  the  temperature  6 
at  any  y,  z  and  t  can  be  expressed.  In  terms  of  the  nodal  temperatures 
[6*]  as  follows 


0  =  (Nl  (0*],  (30) 

Using  Galerkin  method  along  with  assumed  Interpolating  functions, 
one  can  write  the  following  finite  element  equation  (Rao  1988) 


(ktl 

[0*]  +  |lk|l  ♦  [kg]  j  [0«]  =  IP*], 

(31) 

where 

[k*] 

=  1  1  k/a  IN]T  (Nl  dA  , 

A* 

(32) 

[k|] 

=  1  h  [N]^  IN]  dSz  . 

sl 

(33) 

(kal 

=  1  1  k  IB]"^  IB]  dA  , 

A* 

(34) 

and  (P®1  = 

A* 

Q  IN)’’  dA  -  r  0  IN]'^  dSj  +  J  h  0,  [N]^  dS^. 
sl  sl 

(35) 

with 

IB] 

»  r  aNj/ay  aNg/ay  . . .  aN/ay 

(36) 

-  aNj/az  aN2/az  . . .  aNp/az  -I 

where  subscript  p  is  the  number  of  vertices  of  the  element  (p=3  for 
triangular  element). 

The  individual  element  equations  are  assembled  to  form  the  overall 
equation  of  the  NITINOL-composite  beam  equation  which  cam  be  solved  for 
the  nodal  temperatures.  The  solution  is  based  on  a  Cramk-Nlcolson 
trapezoidal  scheme  (Hughes  1977). 


3.4.  NUMERICAL  AND  EXPERIMENTAL  RESULTS 
3.4.1.  BASIC  CHARACTERISTICS  OF  BEAM 


The  characteristics  of  MITINOL-reinforced  beams  are  computed  using 
the  developed  static,  dynamic  and  thermal  models.  The  theoretical 
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predictions  are  compared  with  experimental  results  obtained  with  a 
composite  beam  made  of  randomly  oriented  glass  fibers  embedded  in  a  low 
cure  temperature  polyester  resin.  The  beam  is  30  cm  long,  2. 5  cm  wide 
and  0. 156  cm  thick  mounted  in  clamped-c  lamped  arrangement.  The 
temperature  dependence  of  Young's  modulus  of  the  beam,  shown  in  Figure 
(9),  is  obtained  experimentally  using  the  Dynamic,  Mechanical,  and 
Thermal  Analyzer  (DMTA)  of  Polymer  Laboratories.  Ltd  (1990). 


TEMPERATURE  (  'C  ) 

Figure  (9)  -  Effect  of  operating  temperature  on  Young’s  modulus  of  the 
test  fiberglass  composite  beam 


Four  NITINOL  55  fibers,  that  are  0.55  mm  in  diameter,  are  embedded 
inside  the  beam  through  vulcanized  rubber  sleeves  that  have  outer 
diameter  of  0.95  mm.  Two  sets  of  NITINOL  fibers  were  used.  The  first 
set  consisted  of  trained  fibers  that  have  austenite  transformation 
temperature  of  50  C.  The  second  set  however  is  untrained  and 
accordingly,  the  shape  memory  effect  has  not  been  imparted  to  it.  The 
two  sets  are  inserted,  one  at  a  time,  inside  the  sleeves  and  the  effect 
of  the  shape  memory  and  the  associated  phase  recovery  forces  on  the 
performance  of  the  composite  beam  are  monitored  when  the  beam  is 
exposed  to  different  ambient  temperatures.  The  experimental  set-up, 
shown  in  Figure  (10),  is  placed  inside  a  temperature-controlled  chamber 
to  determine  the  natural  frequencies  of  the  fixed-fixed  beam  as  a 
function  of  the  ambient  temperature.  The  set-up  is  very  similar  to 
that  used  in  studying  the  thermo-dynaunic  characteristics  of  the  NITINOL 
fibers.  However,  instead  of  activating  the  NITINOL  fibers 
electrically,  the  fibers  are  activated  thermally  by  controlling  the 
temperature  of  the  environmental  chamber.  The  measurements  are  carried 
out  after  steady-state  and  thermal  equilibrium  conditions  are  attained. 
Under  these  conditions,  the  composite  matrix  and  the  NITINOL  fibers  are 
all  at  the  same  equilibrium  temperature.  At  each  equilibrium 
temperature,  the  composite  beam  is  subjected  to  random  vibrations  and 
the  resulting  response  is  monitored  by  an  micro-accelerometer  bonded  to 
the  beam.  The  response  is  analyzed  in  the  frequency  domain  to 
determine  the  modes  of  vibration  of  the  composite  beam. 
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Figure  (10)  -  Experimental  set-up  used  in  monitoring  the  performance  of 
NITINOL-reinforced  beams 


3.4.2.  THE  STATIC  AND  DYNAMIC  CHARACTERISTICS 


Figure  (11-a)  shows  the  measured  changes  in  the  first  natural 
frequency  of  the  beam  when  it  is  reinforced  with  untrained  fibers  which 
are  pre-strained  at  different  levels.  The  changes  are  normalized  with 
respect  to  the  natural  frequency  w,,  of  the  un-reinforced  beam  measured 
at  2S°C.  The  normalized  characteristics  of  the  un-reinforced  beam  are 
also  plotted  to  serve  as  a  datum  for  defining  the  effect  of  reinforcing 
the  beam  with  NITINOL  fibers  and  also  the  effect  of  the  pre-strain 
level.  It  can  be  seen  that  the  frequency  of  the  un-reinforced  beam 
drops  as  the  ambient  temperature  increases  and  when  the  temperature 
exceeds  40  C  the  beam  losses  its  elastic  stability  and  start  to  buckle. 
The  drop  in  the  natural  frequency  of  the  un-reinforced  beam  is 
attributed  to  the  softening  of  the  matrix  which  is  clearly  demonstrated 
by  the  loss  in  the  Young’s  modulus  of  the  beam  as  shown  in  Figure  (9). 

Reinforcing  the  beam  with  pre-strained  untrained  NITINOL  fibers 
considerably  increases  the  natural  frequency  of  the  beam.  The  extent 
of  the  upward  shift  in  natural  frequency  increases  with  increased 
pre-strain  level.  An  increase  of  about  40/i  is  obtained  at  room 
temperature  when  the  pre-strain  level  is  only  0.Z6V,.  However,  as  the 
ambient  temperature  increases  the  frequency  shift  drops  in  a  manner 
similar  to  the  characteristics  of  the  plain  un-reinforced  beam.  Such  a 
drop  is  again  attributed  to  the  softening  effect  of  the  matrix  and  the 
fact  that  the  untrained  NITINOL  fibers  act  as  a  static  pre-tensioning 
device  that  produce  constant  tension  which  is  independent  of 
temperature.  Therefore,  the  frequency  enhancement  is  only  generated  by 
the  reinforcement  and  the  pre-strain  effects,  and  not  by  the  shape 
memory  effect. 
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Figure  (11)  -  Effect  of  the  ambient  temperature  and  pre-strain  ievel  on 
enhancing  the  first  mode  of  vibration  of  clamped- clamped 
composite  beam  reinforced  with  NITINOL  fiber  without  (a) 
and  with  (b)  shape  memory. 


ic  IS  importemc  to  note  tnat  consiaeraoiy  nigner  increases  in 
natural  frequency  can  be  obtained  by  further  Increases  of  the 
pre<r£traln  level  to  Its  maximum  permissible  level  of  6%. 

However,  a  greater  frequency  shift  can  be  achieved  by  imparting  the 
shape  memory  effect  to  the  NITINOL  fibers.  The  fibers  are  trained  over 
250  cycles  using  the  procedure  outlined  by  Johnson  (1984).  The  trained 
fibers  are  inserted  into  the  composite  beeim  to  replace  the  untrained 
set.  the  frequency  shifts  become  significant,  particularly  at  high 
ambient  temperatures.  This  is  clearly  demonstrated  in  Figure  (11-b). 
For  temperatures  between  room  temperature  and  40  C,  the  frequency 
shifts  obtained  are  similar  to  those  with  the  untrained  fibers  within 
exgerimental  accuracy.  But,  once  the  ambient  temperature  exceeds  the 
50°C  temperature,  which  is  the  austenite  phase  transformation 
temperature,  the  frequency  shift  characteristics  chzuiges  from  a 
gradually  decaying  trend  and  develops  a  gradually  increasing  profile. 
Such  a  sudden  change  is  a  reflection  of  the  contribution  of  the  phase 
recovery  forces  developed  by  the  shape  memory  effect  which  is 
illustrated  in  Figure  (12). 


Temperature  ('’C) 


Figure  (12)  -  Effect  of  pre-strain  level  and  ambient  temperature  on  the 
phase  recovery  forces  of  trained  NITINOL  fibers. 


The  shape  memory  effect  generates  strain  energy  in  the  NITINOL 
fibers  to  counterbalance  the  softening  effect  of  the  composite  matrix 
with  increasing  temperature.  The  amount  of  strain  energy  developed, 
depends  on  the  initial  pre-strain  level,  it  can  merely  compensate  for 
the  softening  effect  to  maintain  the  beam  frequency  at  nearly  a 
constant  value  which  is  independent  of  ambient  temperature  as  shown  for 
pre-strain  level  of  0.078/i.  It  can  also  increase  the  beam  frequency 
as  the  ambient  temperature  increases  as  indicated  for  pre-strain  levels 
of  0.22  and  0.26^.  With  pre-strain  level  of  0.26%,  the  frec^ency 
increase  reaches  about  70%  of  that  at  ambient  temperature  of  90  C  as 
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compared  to  18%  increase  when  xmtrained  fibers  are  used.  In  this  way, 
the  individual  contributions  of  the  pre-strain,  matrix  softening  and 
shape  memory  effect  on  the  frequency  shift  are  isolated.  This 
facilitates  checking  the  validity  of  the  mathematical  models  against 
the  experimental  results. 

Comparisons  between  the  theoretical  predictions  and  the 
measurements  are  shown  in  Figures  (13-a)  and  {13-b)  for  NITINOL  fibers 
without  and  with  shape  memory  effect,  respectively.  The  figures 
include  comparisons  for  the  first  and  second  modes  of  vibrations. 
Close  agreement  between  theory  auid  experiments  is  evident. 


Figure  (13)  -  Comparison  between  the  theoretical  and  experimental 
frequencies  of  composite  beam  reinforced  with  NITINOL 
fibers  without  (a)  and  with  shape  memory  (b). 


3.4.3.  THERMAL  CHARACTERISTICS  OF  BEAM 


The  thermal  finite  element  model  developed  in  section  3.3.1  is  used 
to  generate  the  temperature  distribution  over  the  beam  cross  section 
which  in  turn  is  utilized  to  compute  the  average  Young’s  modulus  of  the 
beam  under  different  operating  conditions.  Such  thermal 
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characteristics  will  be  presented  in  this  section  and  compared  with 
experimental  results  in  order  to  check  the  validity  of  the  theoretical 
predictions. 

Figure  (14)  shows  a  sample  of  the  finite  element  mesh  used  in  the 
analysis  of  the  t:?mperature  field  in  the  NITINOL-reinforced  beam. 


Figure  (14)  -  The  finite  element  mesh  of  the  thermal  model  of  the 
NITINOL-reinforced  beam. 


Figures  (15),  (16)  and  (17)  show  the  steady-state  temperature  and 
Young’s  modulus  distributions  over  the  cross  section  of  the  beam  when 
activating  all  the  four  fibers,  the  two  extreme  fibers,  and  the  middle 
two  fibers,  respectively.  It  is  evident  that  the  activation  strategy, 
as  well  as  the  number  and  location  of  the  activated  fibers,  result  in 
dramatic  variations  of  the  temperature  and  Young’s  modulus 
distribution.  These  variations  influence  the  static  and  dynamic 
characteristics  of  the  composite  beam.  Hence,  integration  of  the 
thermal  and  mechanical  models  is  essential  to  the  understanding  and  the 
prediction  of  the  behavior  of  NITINOL-reinforced  composites. 

Comparisons  are  shown  in  Figures  (15)  through  (17)  with  the 
experimental  results  when  the  fibers  are  activated  electrically  with 
8.3  watts/fiber  and  when  steady-state  conditions  are  attained  after  720 
seconds.  The  temperature  distribution  is  also  monitored  during  the 
720  second  period  required  to  reach  steady-state  at  six  different 
locations.  Figure  (18)  displays  the  spatial  distribution  of  these 
measurement  stations. 
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ram  raoa  raEn  ram 


Figure  17  -  Temperature  and  Young’s  modulus  distribution  over  beam 
cross  section  when  the  two  middie  NITINOL  fibers  are 

activated 


TEMPEIiATURE  {  “C  ) 


y 


\ 

\ 

\ 

v 

_  Nitinol  i 

Wire 

i  _  Fiberglass 

y'  Beam 

6 

5 

^ - 

1  2*  3*  V 

Figure  (18)  -  Spatial  distribution  of  temperature  measurement  stations 


Figures  (19-a)  and  (19-b)  show  comparisons  between  the  experimental 
and  theoretical  temperature  distributions  over  the  beam  cross  section 
at  three  stations  located  on  the  outer  surface  of  the  beam.  Close 
agreement  between  theory  and  experiment  is  evident. 


THEORY  :  HEATTXG  4  FIBERS 


Figure  (19)  -  Comparison  between  theoretical  amd  experimental 

surface  temperature  distributions  over  beam  cross 
sections  with  all  four  NITINOL  fibers  activated. 


Figure  (20)  shows  the  experimental  surface  temperature  of  the  beam 
along  its  longitudinal  axis.  It  is  clear  that  the  temperature  is 
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practically  constant  along  the  beam  longitudinal  axis.  Such 
measurements  validate  the  assumption  used  in  deriving  the  thermal  model 
and  Justify  the  use  of  the  two-dimensional  simulation  of  the  beam. 


EXPERIMENT  :  HEATING  4  HBERS 


Figure  (20)  -  Experimental  surface  temperature  of  the  beam  along 
longitudinal  axis. 


its 


4.  APPLICATIONS  OF  NITINOL-REINFORCED  COMPOSITE  BEAMS 


The  feasibility  of  utilizing  NITINOL-reinforcing  fibers  to  actively 
control  the  buckling  and  flow— induced  vibrations  of  composite  beams  are 
demonstrated  in  this  section. 

4.1.  ACTIVE  BUCKLING  CONTROL 


The  new  trend  of  designing  light  weight  and  large  structures  render 
these  structures  to  be  more  susceptible  to  failure  due  to  buckling. 
Baz  and  Tampe  (1989)  were  successful  in  to  enhancing  the  buckling 
characteristics  of  long  slender  beams  by  using  external  helical  shape 
memory  actuators.  In  the  present  study,  actuators  in  the  form  of 
NITINOL  fibers  are  embedded  inside  the  long  slender  beams.  With  such 
a  configuration,  beams  can  be  manufactured  from  light  weight  sections 
that  have  built-in  capabilities  for  withstanding  failure  due  structural 
instabilities.  It  was  shown  theoretically  in  section  2  that 
NITINOL-reinforced  composite  beams  can  have  enhanced  buckling 
characteristics  depending  on  the  pre-strain  level  of  the  NITINOL  fibers 
in  comparison  with  the  external  mechanical  and  thermal  loads  acting  on 
these  beams.  The  validation  of  such  theoretical  model  is 
experimentally  demonstrated  in  this  section. 
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4.1.1.  THE  EXPERIMENTAL  SET-UP 


Figure  (21)  shows  a  schematic  drawing  of  the  experimental  set-up 
used  in  actively  controlling  the  buckling  of  a  NITINOL-reinforced  beam. 
The  beam  dimensions  are  63.75  cm  by  2.5  cm  by  0.44  cm.  The  beam  is 
reinforced  by  eight  o.55  mm  NITINOL  fibers  which  are  embedded 
symmetrically  along  the  neutral  axis  of  the  beam. 


Figure  (21)  -  Schematic  drawing  of  the  active  buckling  control 

experiment. 


One  end  of  the  beam  is  clamped  to  a  fixed  base  and  the  other  end  is 
connected  to  the  piston  of  a  loading  cylinder.  The  cylinder  is 
pressurized  by  compressed  air  from  the  storage  tank  of  an  air 
compressor.  The  increasing  compressive  load  applied  by  the  load 
cylinder  to  the  beam  will  eventually  cause  the  beam  to  buckle.  The 
resulting  deflection  of  the  beam  is  monitored  continuously  by  two 
non-contacting  sensors  which  are  placed  on  both  sides  of  the  beam.  The 
sensors  also  serve  as  physical  stops  to  prevent  excessive  deflection 
once  buckling  has  occured.  The  output  signals  of  the  sensors  are  sent 
to  a  micro-computer  via  a  set  of  analog-to-digital  converters.  The 
processing  these  signals  is  shown  in  the  controller  block  diagram  shown 
in  Figure  (22).  When  the  beam  deflection  exceeds  a  pre-set  value  of  a 
dead-band,  the  controller  is  turned  on  using  a  proportional  controller 
with  a  saturation  limit.  The  control  action  is  sent  via  a 
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Volt 


Figure  (22)  -  Block  diagram  of  the  active  buckling  control  system 


digital-to-analog  converter  to  a  power  amplifier  to  activate  the 
NITINOL  fibers  embedded  inside  the  compressively  loaded  composite  beam. 
The  activation  of  the  NITINOL  fibers  will  compensate  for  the  monitored 
deflection  and  the  phase  recovery  forces  developed  in  the  fibers  will 
attempt  to  bring  the  beam  back  to  its  undeflected  position. 

It  is  important  here  to  note  that  the  controller  dead-band  is 
essential  to  prevent  chattering  of  the  controller  as  observed  by  Baz, 
Iman  and  McCoy  (1990).  Also,  the  saturation  limiting  of  the  maximum 
voltage  applied  to  the  NITINOL  fibers  is  necessary  to  avoid  destruction 
of  the  shape  memory  effect  of  the  NITINOL  fibers  due  to  excessive 
heating. 

In  the  active  buckling  control  system  described,  the  NITINOL 
reinforcing  fibers  are  clamped  at  one  end  to  the  fixed  base  and  at  the 
other  end  to  pre-tensioning  cylinder  via  a  load  cell.  The  load  cell 
monitors  the  initial  value  of  the  pre-tension  applied  to  the  NITINOL 
fibers  by  the  pre-tensioning  cylinder.  The  load  cell  also  continuously 
measures  the  phase  recovery  forces  developed  in  the  NITINOL  fibers  as 
they  undergo  their  phase  transformation. 


4.1.2.  EXPERIMENTAL  RESULTS 


Figure  (23)  shows  a  comparison  of  the  performance  with  and  without 
the  active  buckling  control.  The  results  displayed  are  for  a 
NITINOL-reinforced  beam  with  each  of  the  eight  fibers  has  an  initial 
tension  of  33.7  N  which  corresponds  to  an  initial  pre-strain  of  0.3SX. 
For  the  controlled  cases,  the  controller  dead  baind  corresponds  to 
deflection  error  of  0.0176  mm  and  the  controller  gain  is  2727  volt/mm. 
The  saturation  limit  is  6  volts/fiber  and  the  maximum  current  is  1.6A. 


In  Figure  (23-a),  the  deflection  of  the  beam  resulting  from  the 
application  of  a  gradually  increasing  compressive  load  is  shovm.  The 
rate  of  load  increase  is  500  N/min  as  shown  in  Figure  (23-b).  It  can 
be  observed  that  the  uncontrolled  beam  buckles  when  the  compressive 
axial  loads  starts  exceeding  330  N.  However,  when  all  eight  NITINOL 
fibers  are  activated,  the  beam  can  withstand  axial  loads  up  to  950  N 
before  it  starts  to  buckle.  Therefore,  with  the  NITINOL  reinforcement, 
it  is  possible  to  triple  the  critical  buckling  load  of  the  beam.  With 
four  activated  fibers,  the  critical  buckling  load  is  about  700  N  which 
corresponds  to  about  double  the  critical  load  of  the  uncontrolled  beam. 
Figure  (23-c)  shows  the  corresponding  time  history  of  the  temperature 
of  the  NITINOL  fibers  due  to  the  activation  and  de-activation  of  the 
controller.  In  Figure  (23-d),  the  time  history  of  the  phase  recovery 
forces  developed  in  the  NITINOL  fibers  is  shown.  The  tension  in  the 
fibers  remains  equal  to  the  initial  pre-tension,  i.e.  33. 7N  x  8  = 
269. 6N,  for  the  uncontrolled  case.  However,  the  tension  increases  to 
approximately  1000  N  when  eight  NITINOL  fibers  are  activated. 

With  such  capability,  it  is  possible  to  energize  different  sets  of 
the  NITINOL  fibers  to  counterbalance  the  external  loading  condition  and 
avoid  buckling  of  the  composite  beam.  Therefore,  for  small  external 
loads,  it  is  only  necessary  to  energize  a  few  fibers,  but  as  the  load 
increases,  the  controller  can  energize  a  larger  number  of  fibers  to 
maintain  the  beam  in  its  undefelected  form. 

The  effect  of  varying  the  controller  parameters  on  the  performance 
of  the  active  control  system  is  shown  in  Figures  (24),  (25)  and  (26). 
In  Figure  (24),  the  effect  of  varying  the  control  dead-band  on  the 
system  performance  is  shown.  In  this  case,  the  controller  will  be  off 
until  the  beam  deflection  exceeds  the  dead  band.  Once  the  deflection 
exceeds  the  dead  band  ,  the  controller  is  energized.  This  is 
accomplished  by  the  heating  of  the  NITINOL  fibers  and  the  development 
of  the  phase  recovery  forces  as  shown  in  Figures  (24-c)  and  (24-d), 
respectively.  For  the  range  of  dead  bands  considered,  between  0.0176 
mm  and  0.528  mm.  the  effect  on  the  critical  buckling  loads  is 
insignif leant. 

The  effect  of  varying  the  controller  gain  from  2727  volt/mm  to 
136.35  volt/mm  on  the  system  performance  is  shown  in  Figure  (25).  This 
effect  varies  the  slope  of  temperature  rise  of  the  NITINOL  fibers  and, 
in  turn,  the  rate  at  which  the  corresponding  phase  transformation 
forces  are  recovered.  Changing  the  controller  gain  is  found  to 
influence  to  some  extent  the  critical  buckling  load.  For  example,  when 
the  controller  gain  is  2727  voLt/mm  the  critical  buckling  load  is  950 
N  and  when  the  gain  is  reduced  to  136.35  volt/mm  the  critical  buckling 
load  become  about  850.  Therefore,  reducing  the  controller  gain  by  a 
factor  of  1/20  only  results  in  a  10.5%  reduction  in  the  critical 
buckling  load. 

The  effect  that  the  pre-tension  has  on  the  system  performance  is 
shown  in  Figure  (26)  for  dead  band  of  0.0176  nun  and  controller  gain  of 
2727  volts/mm.  It  is  clear  from  the  results  obtained  that  ,  the 
pre-tension  plays  the  most  crucial  role  in  controlling  the  buckling  of 
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the  beam.  Increasing  the  tension  from  0.0  N/fiber  to  33.7  N/fiber 
increases  the  critical  buckling  load  from  about  100  N  to  950  N, 
respectively. 

4.1.3.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENTS 

The  mechanism  of  actively  concrolling  the  buckling  of  the 
NITINOL-reinforced  beam  can  best  be  understood  by  considering  Figure 
(27-a).  The  figure  represents  the  theoretical  prediction  of  the 
buckling  characteristics  of  actively  controlled  NITINOL-reinforced 
beams.  In  the  figure,  the  applied  axial  load  is  increased  gradually  at 
a  linear  rate  of  500  N/min.  For  the  uncontrolled  beam,  the  critical 
buckling  load  is  fixed  at  320  N  and  remain  unchanged  with  time.  This 
load  corresponds  to  a  pre-tension  of  33.5  N/fiber.  When  the  applied 
load  becomes  equal  to  the  critical  buckling  load,  the  beam  is  on  the 
verge  of  elastic  instability.  The  beam  will  buckle  once  the  applied 
load  exceeds  the  fixed  critical  buckling  load.  For  the  case  of 
controlled  beam,  the  beam  starts  at  time  t  =  0  with  the  same  buckling 
load  as  the  uncontrolled  beam.  When  the  controller  senses  any 
deflection  greater  than  the  dead  band  due  to  the  application  of  the 
external  load,  the  buckling  characteristics  of  the  beam  is  enhanced  as 
represented  by  the  dashed  characteristics.  The  activation  of  the 
NITINOL  fibers  makes  the  beam  less  susceptible  to  buckling  as  the 
critical  buckling  load  is  increased  to  become  850  N  instead  of  the 
original  uncontrolled  load  of  320  N.  Accordingly,  the  controlled  beam 
will  not  buckle  until  the  applied  load  exceeds  the  theoretically 
predicted  limit  of  850  N. 

The  effect  of  varying  the  pre-tension  levels  on  the  theoretical 
prediction  of  the  critical  buckling  load  is  shown  in  Figure  (27-b) 
along  with  the  corresponding  experimental  results.  It  is  evident  that 
there  is  a  close  agreement  between  theory  and  experiment. 

4.2.  ACTIVE  CONTROL  OF  FLOW-INDUCED  VIBRATIONS 

The  phenomenon  of  vibrations  induced  by  the  flow  of  fluids  past 
flexible  structures  has  been  of  concern  for  many  years.  This  concern 
is  attributed  to  the  detrimental  effects  that  such  vibrations  can  have 
on  the  integrity  of  these  structures.  Several  attempts  have  been  made 
to  passively  and  actively  control  the  flow-induced  vibrations  of 
various  structural  members.  For  example,  Baz  and  Ro  (1991)  utilized  a 
direct  velocity  feedback  controller  to  control  the  vortex-induced 
vibrations  of  a  flexible  cylinder.  The  control  system  relied  in  its 
operation  on  an  electromagnetic  actuator  to  provide  the  control  action 
necessary  to  resist  the  flow-induced  vibration.  Baz  and  Kim  (1992) 
developed  a  modal  space  control  method  to  control  the  vortex-induced 
vibrations  of  a  flexible  cylinder  using  piezo-electric  actuators. 


CRITICAL  BUCKLING  LOAD  (  N  ) 


A 


Figure  (27-a)  -  Theoretical  prediction  of  critical  buckling  load  of 
controlled  and  uncontrolled  beams. 


INITIAL  FIBER  PRE -TENSION  {  N  ) 

Figure  (27-b)  -  Comparison  between  theoretical  and  experimental 

critical  buckling  loads  as  function  of  the  initial 
fiber  pre-tension. 
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In  this  study,  NITINOL  fibers  are  used  as  embedded  actuators  to 
control  the  flow-induced  vibrations  of  NITINOL-reinforced  composite 
beams.  The  control  action  used  is  of  the  ON-OFF  type.  When 
undesirable  vibrations  are  detected,  the  NITINOL  fibers  are  activated. 
The  strain  energy  generated  in  the  fibers  by  the  phase  transformation 
makes  the  beam  stiffer  and  shifts  its  natural  frequencies  away  from  the 
excitation  frequency.  This  modal  tuning  of  the  beam  dynamics  in 
response  to  the  external  disturbances  can  result  in  reducing  the 
amplitude  of  vibration  of  the  beam  to  acceptable  limits.  This  will 
only  be  true  if  the  flow-induced  vibrations  have  a  narrow  frequency 
band.  For  broad  band  excitations,  the  modal  tuning  mechanism  will  not 
be  effective  in  attenuating  the  induced  vibrations  as  it  merely  shifts 
the  vibration  energy  to  higher  frequency  bands.  However,  in 
NITINOL-reinforced  composites  an  additional  mechanism  can  play  a 
dominant  role  in  the  suppression  of  broad  band  vibrations.  This 
mechanism  is  generated  by  the  temperature-dependent  damping 
characteristics  of  the  composite  matrix  as  shown  in  Figure  (28).  These 
characteristics  are  obtained  experimentally  using  the  Dynamic, 
Mechanical  ,  and  Thermal  Analyzer  (DMTA). 


Temperature  {°C) 

Figure  (28)  -  Loss  coefficient  of  fiberglass  composite  beam 


When  the  NITINOL  fibers  are  activated,  the  temperature  of  the 
matrix  increases  as  Indicated  in  section  3.4.3.  Such  a  temperature 
increase  is  accompanied  with  an  increase  in  the  loss  coefficient  of  the 
matrix  as  displayed  in  Figure  (28).  Operation  at  a  temperature 
corresponding  to  the  maximum  loss  coefficient  is  essential  to  achieve 
maximum  structural  damping.  At  that  temperature,  the  dissipation  of 
the  vibration  energy  will  be  maximum  and  the  attenuation  of  broad  band 
vibrations  will  also  be  maximum. 
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therefore.  In  NITINOL-reinforced  composites,  the  interaction 
between  the  modal  timing  and  the  enhanced  damping  characteristics  is 
crucial  in  controlling  both  narrow  and  broad  band  vibrations. 


4.2.1.  TEST  BEAM  AND  FACILITY 


Figure  (29)  shows  a  schematic  drawing  of  the  test  beam  mounted  in  a 
clamped-clamped  configuration  inside  a  low  speed  wind  tunnel.  The  beam 
is  30  cm  long,  2.5  cm  wide  and  0.156  cm  thick.  It  is  made  of 
fiberglass/polyester  resin  composite  with  four  embedded  NITINOL  fibers. 
The  fibers  are  0. 55  mm  in  diameters  and  are  inserted  inside  0. 95  mm 
vulcanized  rubber  sleeves.  The  elastic  and  damping  characteristics  of 
the  beam  are  shown  in  Figures  (9)  and  (28)  respectively. 

The  beam  is  mounted  inside  a  30  cm  x  30  cm  test  section  of  a  low 
speed  wind  tunnel  and  is  subjected  to  flow  speeds  between  4. 8  and  8. 3 
m/s. 
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Figure  (29)  -  A  schematic  drawing  of  a  NITINOL-reinforced  beam  mounted 
inside  a  low-speed  wind  tunnel. 


The  NITINOL  fibers  have  an  initial  pre-tension  of  17.6  N  each  and 
are  electrically  activated  by  applying  a  voltage  of  4. 5  V  across  each 
fiber.  This  generates  a  current  of  1.85  A  and  a  total  of  8.325  watts 
are  dissipated  in  the  composite  beam.  The  resulting  shift  of  the  first 
three  natural  frequencies  of  the  beam  are  shown  in  Figure  (30).  The 
figure  shows  the  time  history  of  the  frequency  shift  when  all  the  four 
fibers  are  activated  for  two  minutes  and  then  de-activated  for  the 
remaining  six  minutes. 

The  effect  of  varying  the  number  of  activated  fibers  on  the 
frequency  shift  of  the  first  three  modes  is  shown  in  Figure  (31-a) 


Acj/wq  (x'100%) 
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Figure  (30)  -  The  shift  of  the  first  three  modes  of  the 

NITINOL-reinforced  beaim  during  the  activation  and 
de-activation  cycles  of  all  the  four  NITINOL  fibers 


Number  of  activated  fibers  Number  of  activated  fibers 


Figure  (31)  -  Effect  of  number  of  activated  NITINOL  fibers  on  the 
frequency  shift  and  damping  ratio  of  the  first  three 
modes  of  the  NITINOL-reinforced  beam. 
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after  steady-state  conditions  are  attained.  The  results  displayed  in 
the  figure  are  normalized  with  respect  to  the  natural  frequencies 
measured  at  25°C.  The  corresponding  variation  of  the  modal  damping  is 
shovm  in  Figure  (31-b).  The  figure  clearly  indicates  that  increasing 
the  number  of  activated  fibers  results  in  enhancing  the  damping 
characteristics  of  the  beam.  This  is  attributed  to  the  increase  in  the 
beam  temperature  when  the  number  of  activated  fibers  is  increased. 
This  in  turn  makes  the  composite  matrix  of  the  beam  operate  near  the 
point  of  maximum  loss  coefficient  shown  in  Figure  (28). 


4.2.2.  EXPERIMENTAL  RESULTS 


The  performance  of  the  NITINOL-reinforced  composite  beam  subjected 
to  flow-induced  vibrations  is  monitored  at  different  flow  speeds  with 
and  without  the  activation  of  the  NITINOL  fibers.  The  experiments  aim 
at  demonstrating  the  feasibility  of  NITINOL  fibers  in  attenuating  the 
flow-induced  vibrations. 

Figure  (32)  shows  the  spectra  of  the  amplitude  of  vibration  of  the 
beam  as  measured  at  different  flow  speeds  ranging  between  4.8  m/s  and 
8.32  m/s.  The  figure  also  shows  a  comparison  between  the  amplitudes  of 
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Figure  (32)  -  Spectra  of  the  amplitudes  of  vibration  of  the 

NITINOL-reinforced  beam  at  different  flow  speeds  with 
and  without  the  activation  of  the  four  NITINOL  fibers. 
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Figure  (34)  -  Effect  of  activating  different  numbers  of  NITINOL  fibers 
on  the  spectra  of  the  amplitude  of  flow-induced  vibration 
at  flow  speed  =  8.32  m/s 


Number  of  activated  fibers 


Figure  (35)  -  Effect  of  number  of  activated  NITINOL  fibers  on  the 
amplitude  of  flow-induced  vibration  at  flow  speed  of  8.32 
m/s 


Emphasis  has  been  placed  in  the  presentation  on  the  actuation 
capabilities  of  the  NITINOL  fibers.  Extensive  efforts  are,  however,  in 
progress  to  use  the  NITINOL  fibers  to  extract  modal  and  physical 
displacements  of  structures  with  multi-modes  of  vibration  (Baz,  Poh  and 
Gilheany  1991). 

With  such  built-in  sensing  and  controlling  capabilities, 
NITINOL-reinforced  composites  can  provide  a  means  for  continuously 
tuning  the  structural  characteristics  to  adapt  to  changes  in  the 
operating  conditions.  These  features  will  be  particularly  useful  in 
mamy  critical  structures  that  are  intended  to  operate  autonomously  for 
long  durations  in  isolated  environment  such  as  defense  vehicles,  space 
structures  and  satellites. 
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NOMENCLATURE 


[A]  interpolating  function  of  beam  deflection 

(All  ith  element  of  [A] 

A,,  cross  sectional  area  of  beam 

[B]  matrix  of  first  derivatives  of  the  nodal  interpolating 
functions 

c  wave  propagation  speed 

IC],  Id]  first  and  second  derivatives  of  interpolating  function 
of  beam  deflection 
Young' s  modulus  of  beam 

[F]  vector  of  external  loads  acting  qn  beam 

h  convective  heat  transfer  coefficient 

I,  area  moment  of  inertia  of  beam 

k  thermal  conductivity  of  beam 

[k"^2.3l  matrices  given  by  equations  (32),  (33)  and  (34) 
respectiveiy 

(K^l  stiffness  matrix  of  beam  element 

ty,  direction  cosines  of  outward  normals  to  beam  boundaries 

L  length  of  beam  element  aind  NITINOL  fiber 

Mj  external  moment  acting  at  ith  node 

[M^I  mass  matrix  of  beam  element 

m^(i,j)  the  element  i,j  of  the  mass  matrix 

[N]  interpoiating  function  of  beam  temperature 

Nj  interpoiating  fxonction  of  ith  node 

p  number  of  vertices  of  element 

mechanical,  net  and  thermal  axial  loads  acting  on  beam 
IP®]’  matrix  given  by  equation  (35) 
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ABSTRACT 

The  buckling  characteristics  of  flexible  fiberglass  composite  beams 
are  actively  controlled  by  activating  optimal  sets  of  a  shape  memory  alloy 
(NITINOL)  wires  which  are  embedded  along  the  neutral  axes  of  the  beams. 
With  such  active  control  capabilities,  the  beams  can  be  manufactured  from 
light  weight  sections  without  compromising  their  elastic  stability.  This 
feature  will  be  invaluable  in  building  light  weight  structures  that  have 
high  resistance  to  failure  due  to  buckling. 

A  finite  element  model  is  developed  to  analyze  the  individual 
contributions  of  the  fiberglass-resin  laminate,  the  NITINOL  wires,  and  the 
shape  memory  effect  to  the  overall  performance  of  the  composite  beams.  A 
closed-loop  computer-controlled  system  is  built  to  validate  the  finite 
element  model.  The  system  is  used  to  control  the  buckling  of  a  fiberglass 
polyester  resin  beam  which  is  63.75  cm  long,  0.45  cm  thick  and  2.54  cm  wide 
reinforced  with  8  NITINOL  -  55  wires  that  are  0.55  mm  in  diameter.  The 
results  obtained  confirm  the  developed  theoretical  model  and  indicate  that 
the  critical  buckling  load  can  be  increased  three  times  when  compared  to 
the  uncontrolled  beam. 


1,  INTRODUCTION 

Considerable  attention  has  been  devoted  recently  to  the  utilization  of 
the  Shape  Memory  Nickel-Titanium  alloy  (NITINOL)  in  developing  SMART 
composites  that  are  capable  of  adapting  intelligently  to  external 
disturbances  (Ikegami  et  al.  1990,  Rogers  et  al.l99l,  and  Baz  et  al.  1990 
and  1991).  Such  wide  acceptance  stems  from  the  fact  that  NITINOL  acts  as  an 
actuator  converting  thermal  energy  into  mechanical  energy  (Perkins  1975  and 
Duerig  et  al.  1990)  as  it  undergoes  its  unique  phase  transformation  from 
low  temperature  martensite  to  high  temperature  austenite.  During  this 
phase  transformation  process  large  phase  recovery  forces  are  generated  and 
thereby  alter  the  strain  energy  of  the  composite  inside  which  the  NITINOL 
fibers  are  embedded.  With  such  capabilities,  the  static  and  dynamic 
performance  of  the  SMART  composites  can  be  optimized  and  tailored  to  match 
changes  in  the  operating  conditions. 

Emphasis  is  placed,  in  the  present  work,  on  using  the  shape  memory 
effect  of  the  NITINOL  fibers  In  controlling  the  buckling  of  fiberglass 
composite  beams.  The  NITINOL  fibers  are  embedded  Inside  vulcanized  rubber 
sleeves  placed  along  the  neutral  axes  of  these  composite  beams  as  shown  In 
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Figure  (1).  In  this  arrangement,  the  fibers  are  free  to  move  during  the 
phase  transformation  process  in  order  to  avoid  degradation  and/or 
destruction  of  the  shape  memory  effect  which  may  result  when  the  fibers  are 
completely  bonded  inside  the  composite  matrix.  The  NITINOL  fibers  are 
trained  to  memorize  the  shape  of  the  unbuckled  beam  and  when  the  beam  is 
deflected  under  the  action  of  external  compressive  loads,  the  controller 
activates  the  NITINOL  fibers  by  heating  them  above  their  transformation 
temperature.  The  generated  phase  recovery  forces  bring  the  beam  back  to 
its  memorized  undeflected  position.  The  present  study  is  motivated  by  the 
work  of  Baz  and  Tampe  (1989)  where  external  helical  shape  memory  actuators 
are  used  to  enhance  the  buckling  characteristics  of  long  slender  beams. 


urtecrums  »tnxcTn 

^  LWeONTKOUfD  — »  COffTHOLLtD 


Figure  (1)  -  Principle  of  buckling  control  of  NITINOL-reinforced  composites 

\ 

2.  STATIC  CHARACTERISTICS  OF  NITINOL-REINFORCED  BEAMS 

The  static  characteristics  of  NITINOL-reinforced  composite  beams  are 
primarily  governed  by  their  stiffness.  The  overall  beam  stiffness  is  made 
up  of  the  following  components:  the  flexural  rigidity  of  the  beam,  the 

geometric  stiffness  that  accounts  for  the  axial  and  thermal  loading  as  well 

^  as  the  stiffness  imparted  by  the  elasticity  of  the  NITINOL  fibers.  These 

individual  components  of  the  beam  stiffness  can  be  determined  by 

considering  the  NITINOL-reinforced  beam  element  shown  in  Figure  (2)  with 
the  forces  acting  on  it  and  the  associated  displacements.  The  combined 
stiffness  of  the  element  can  be  obtained  using  the  principle  of 

conservation  of  energy  and  equating  the  work  done  by  external  loads  to  the 
^  strain  energies  stored  in  the  element.  In  the  present  analysis,  the  theory 

of  Bernoulli-Euler  beams  is  used  with  the  assumption  of  small  deflections. 


(a)  (b) 

Figure  (2)  -  NITINOL-reinforced  beam  element  with  forces  and  displacements 
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2. 1  External  Work 

The  work  done  by  the  external  loads  includes: 

a.  work  done  by  transverse  loads  and  moments  (W^):  This  work  is  given  by 

=  1/2  [6]T  (F],  (1) 

where  [5]  is  deflection  vector  and  [F]  is  transverse  load  vector. 

b.  work  done  by  the  axial  mechanical  loads  (V^a):  given  by 

L 

Wja  »  P,  /  2  J  (dw/dx)2  dx.  (2) 

0 

where  is  the  external  axial  compressive  load  acting  on  the  beam. 

c.  work  done  by  the  axial  thermal  loads  (U2t):  ^2t  represents  the  work 

done  by  the  thermal  loads  P^  on  the  beam  element  due  to  changes  in  the 
temperature  A9  of  the  element  caused  by  changes  in  the  ambient  temperature 
or  during  the  activation  and  de-activation  of  the  NITINOL  fibers.  It  is 
given  by 

L 

Wjt  =  Pt  /  2  I  (dw/dx)^^  dx.  (3) 

0 

where  is 

P^  =  a  A0  E„  A,.  (4) 

where  a  is  the  thermal  expansion  coefficient  of  the  composite,  £„  is  its 
modulus  of  elasticity  and  A„  is  the  beam  cross  sectional  area. 

2.2  Stored  strain  energy 

The  stored  strain  energy  consists  of  two  components: 

a.  strain  energy  of  beam  (W3);  The  energy  stored  in  the 

to  its  bending  is  given  by  ^ 

Wj  =  E,  I„  /  2  J  (d2w/dx2)2  dx, 

0 

where  I„  is  the  flexural  rigidity  of  the  beam. 

b.  strain  energy  of  NITINOL  fibers  (W^);  Considering  the  NITINOL  fiber  as 
a  string  with  a  tension  T  which  is  displaced  laterally  a  distance  w  from 
the  neutral  axis  of  the  beam.  Then  its  stored  strain  energy  W^  is  given  by 

L 

W4  =  T  /  2  J  (dw/dx)2  dx.  (6) 

0 

Equating  the  sum  of  the  work  done  by  the  external  forces  F,  P,  and  Pj, 
to  the  sura  of  the  strain  energies  stored  in  the  beam  and  the  NITINOL  fibers 
gives 


beam  element  due 

(5) 


U, 


'at  -  W3  *  W,. 


(7) 
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Substituting  equations  (1),  (2),  (3).  (5)  and  (6)  into  equation  (7) 

yields 


L  L 

[5]^[F]  =  E.  I.  J  (d^w/dx^)^  dx  -  P„  J  (dw/dx)^  dx,  (8) 

0  0 
where  P„  is  the  net  axial  force  give  by 


Pn  =  (  P-  Pt  -  T). 


(9) 


Defining  a  proper  displacement  function  for  the  composite  beam  element, 
one  can  write  the  beam  deflection  w  as 


w  =  (A)  [5],  (10) 

where  the  elements  of  matrix  [A]  are  function  of  x  (Fenner  1975). 

Accordingly,  dw/dx  and  d^/dx^  can  be  obtained  by  differentiating 
equation  (10)  with  respect  to  x  to  yield 

dw/dx  =  [C]  (51  and  d^/dx^  [D]  [5].  (11) 

If  the  stiffness  matrix  [K^l  of  the  element  is  defined  by  the  following 
relationship 


(FI  =  (KJ  (51,  (12) 

then,  (Kg)  can  be  determined  by  combining  equations  (8),  (11)  and  (12)  as 
follows 

L  L 

(K,)  =  E„  I„  J  (D]T(Dl  dx  -  P„  I  (C]T(cj  dx.  (13) 

0  0 

The  element  stiffness  matrix  (Kg)  of  equation  (13)  consists  of  two 
components;  the  conventional  transverse  stiffness  and  the  geometric 
stiffness  that  combines  the  effect  of  the  axial  mechanical  loads,  axial 
thermal  loads  and  the  tension  of  the  NITINOL  reinforcing  fibers.  Equation 
(13)  also  represents  the  basic  equation  for  understanding  the  role  that  the 
NITINOL  fibers  can  play  in  controlling  the  static  characteristics  of  the 
composite  beam.  For  example,  if  the  beam  is  not  reinforced  by  NITINOL 
fibers  (i.e.  T  =  0)  and  the  mechanical  and  thermal  loads  induce  compressive 
stresses  in  the  beam,  then  the  geometric  stiffness  will  increase  and  the 
total  element  stiffness  will  decrease.  When  the  combined  effect  of  the 
mechanical  and  thermal  loads  reaches  a  critical  magnitude  such  that  the 
geometric  stiffness  becomes  equal  to  the  flexural  stiffness  of  the  beam, 
the  beam  stiffness  vanishes  and  the  beam  becomes  elastically  unstable. 
Subjecting  the  beam  to  any  additional  external  disturbance  will  cause  the 
beam  to  buckle. 

It  should  be  pointed  out  that  the  thermal  loading,  as  it  Increases  the 
geometric  stiffness,  also  decreases  the  flexural  stiffness  of  the  beam 
because  it  reduces  its  effective  modulus  of  elasticity  E^.  Such  a  dual 
effect  makes  the  beam  buckle  under  smaller  thermal  loads  than  under  pure 
mechanical  loading. 
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However,  the  critical  load  of  the  un-relnforced  beam  can  be  Increased 
by  embedding  pre-strained  NITINOL  fibers  into  the  beam.  If  the  tension  T, 
resulting  from  the  pre-strain  alone,  is  high  enough  to  counter-balance  the 
mechanical  and  thermal  effects  then  the  beam  stiffness  can  be  maintained 
unchanged.  For  higher  pre-strain  levels,  the  beam  stiffness  can  be 
enhanced.  Further  enhancement  can  be  achieved  when  the  shape  memory  effect 
of  the  NITINOL  fibers  is  activated  by  heating  the  fibers  above  their  phase 
transformation  temperature.  The  adoitional  tension,  induced  into  the 
fibers  by  the  phase  recovery  forces,  makes  the  net  axial  load  negative 
and  accordingly  increases  the  overall  stiffness  of  the  beam  element. 
However,  it  is  essential  that  the  total  tension  in  the  NITINOL  fibers, 
i.e.  .  the  sum  of  the  tension  due  to  the  pre-strain  and  the  phase  recovery 
force,  must  exceed  the  mechanical  and  thermal  loads  and  compensate  for  the 
softening  effect  in  the  matrix  resulting  from  heating  the  NITINOL  fibers 
inside  the  composite  matrix. 

Therefore,  effective  control  of  the  stiffness  of  NITINOL-reinforced 
composites  can  be  achieved  by  proper  selection  of  the  initial  pre-strain 
level  of  the  NITINOL  fibers.  This  selection  is  particularly  crucial  as  the 
pre-strain  level  determines  the  generated  levels  of  recovery  forces. 

3.  THE  EXPERIMENTAL  SET-UP  AND  RESULTS 

3.1.  Experimental  set-up 

Figure  (3)  shows  a  schematic  drawing  of  the  experimental  set-up  used  in 
actively  controlling  the  buckling  of  a  NITINOL-reinforced  beam.  The  beam 
dimensions  are  63.75  cm  by  2.5  cm  by  0.44  cm.  The  beam  is  reinforced  by 
eight  0.55  mm  NITINOL  fibers  which  are  embedded  symmetrically  along  the 
neutral  axis  of  the  beam. 
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Figure  (3)  -  Schematic  drawing  of  the  active  buckling  control  experiment. 


The  right  end  of  the  beam  is  clamped  to  a  fixed  base  and  the  left  end 
is  connected  to  the  piston  of  a  loading  cylinder.  The  cylinder  is 
pressurized  by  compressed  air  from  the  storage  tank  of  an  air  compressor. 
The  increasing  compressive  load  applied  by  the  load  cylinder  to  the  beam 
will  eventually  cause  the  beam  to  buckle.  The  resulting  deflection  of  the 
beam  is  monitored  continuously  by  a  non-contacting  sensor  which  is  placed 
at  the  mid-span  of  the  beam.  The  sensor  also  serves  as  physical  stop  to 
prevent  excessive  deflection  once  buckling  has  occured.  The  output  signals 
of  the  sensor  is  sent  to  a  micro-computer  via  a  set  of  analog-to-dlgital 
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converters.  The  processing  of  the  sensor’s  signal  is  shovm  in  the 
controller  block  diagram  shown  in  Figure  (4).  When  the  beam  deflection 
exceeds  a  pre-set  value  of  a  dead-band,  the  controller  is  turned  on  using  a 

V 


DB  Yt-'f 


Figure  (4)  -  Block  diagram  of  the  active  buckling  control  system 

proportional  controller  with  a  saturation  limit.  The  control  action  is 
sent  via  a  digital-to-analog  converter  to  a  power  amplifier  to  activate  the 
NITINOL  fibers  embedded  Inside  the  compressively  loaded  composite  beam. 
The  activation  of  the  NITINOL  fibers  will  compensate  for  the  monitored 
deflection  and  the  phase  recovery  forces  developed  in  the  fibers  will 
attempt  to  bring  the  beam  back  to  its  undeflected  position. 

In  the  active  buckling  control  system  described,  the  NITINOL 
reinforcing  fibers  are  clamped  at  one  end  to  the  fixed  base  and  at  the 
other  end  to  pre-tensioning  cylinder  via  a  load  cell.  The  'load  cell 
monitors  the  initial  value  of  the  pre-tension  applied  to  the  NITINOL  fibers 
by  the  pre-tensioning  cylinder.  The  load  cell  also  continuously  measures 
the  phase  recovery  forces  developed  in  the  NITINOL  fibers  as  they  undergo 
their  phase  transformation. 

3.2.  Experimental  Results 

Figure  (5)  shows  a  comparison  of  the  performance  with  and  without  the 
active  buckling  control.  The  results  displayed  are  for  a 
NITINOL-reinforced  beam  with  each  of  the  eight  fibers  has  an  initial 
tension  of  33.7  N  which  corresponds  to  an  Initial  pre-strain  of  0.  35/i.  For 
the  controlled  cases,  the  controller  dead  band  corresponds  to  deflection 
error  of  0.0176  mm  and  the  controller  gain  is  2727  volt/mm.  The  saturation 
limit  is  6  volts/fiber  and  the  maximum  current  is  1.6A. 

-  '  ;  4  flWr*  ^  ;  •  FtWe* 


Figure  (5)  -  Performance  of  buckling  controller  with  0,  4  and  8  NITINOL 
fibers  activated  (dead  band  =  0.0176  ram,  gain  *  2727  volts/mm). 
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Figure  (6)  -  Effect  of  dead  band  (gain  »  2727  volts/mm,  initial 

pre- tension/fiber  =  33. 7N  and  8  fibers  are  activated). 


Figure  (7)  -  Effect  of  controller  gain  (dead  band  =  0.0176  mm,  initial 
pre-tension/fiber  *  33. 7N  and  8  fibers  are  activated). 
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In  Figure  (5-a),  the  deflection  of  the  beam  resulting  from  the 
application  of  a  gradually  increasing  compressive  load  is  shown.  The  rate 
of  load  increase  is  500  N/min  as  shown  in  Figure  (5-b).  It  can  be  observed 
that  the  uncontrolled  beam  buckles  when  the  increasing  compressive  axial 
load  exceeds  330  N.  However,  when  all  eight  NITINOL  fibers  are  activated, 
the  beam  can  withstand  axial  loads  up  to  950  N  before  it  begins  to  buckle. 
Therefore,  with  the  NITINOL  reinforcement  it  is  possible  to  almost  triple 
the  critical  buckling  load  of  the  beam.  With  four  activated  fibers,  the 
critical  buckling  load  is  about  700  N  which  corresponds  to  about  double  the 
critical  load  of  the  uncontrolled  beam.  Figure  (5-c)  shows  the 
corresponding  time  history  of  the  temperature  of  the  NITINOL  fibers  due  to 
the  activation  and  de-activation  of  the  controller.  In  Figure  (5-d),  the 
time  history  of  the  phase  recovery  forces  developed  in  the  NITINOL  fibers 
is  shown.  The  tension  in  the  fibers  remains  equal  to  the  initial 
pre-tension,  i.e.  33. 7N  x  8  =  269. 6N,  for  the  uncontrolled  case.  However, 
the  tension  increases  to  approximately  1000  N  when  eight  NITINOL  fibers 
are  activated. 

It  is  also  possible  to  energize  different  sets  of  the  NITINOL  fibers  to 
counterbalance  the  external  loading  condition  in  order  to  prevent  buckling 
of  the  composite  beam.  For  small  external  loads  it  is  only  necessary  to 
energize  a  few  fibers,  but  if  the  load  increases  the  controller  can 
energize  a  larger  number  of  fibers  to  maintain  the  beam  in  its  undefelected 
form. 

The  effect  of  varying  the  controller  parameters  on  the  performance  of 
the  active  control  system  is  shown  in  Figures  (6).  (7)  and  (8).  In  Figure 
(6),  the  effect  of  varying  the  control  dead-band  on  the  system  performance 
is  shown.  In  this  case,  the  control  action  is  only  generated  when  the  beam 
deflection  exceeds  the  dead  band.  This  is  accomplished  by  the  heating  of 
the  NITINOL  fibers  and  the  development  of  the  phase  recovery  forces  as 
shown  in  Figures  (6-c)  and  (6-d),  respectively.  For  the  range  of  dead 
bands  considered,  between  0.0176  mm  and  0.528  mm,  the  effect  on  the 
critical  buckling  loads  is  insignificant. 

The  effect  of  decreasing  the  controller  gain  from  2727  volt/mm  to 
136.35  volt/mm  on  the  system  performance  is  shown  in  Figure  (7).  This 
effect  decreases  the  slope  of  temperature  rise  of  the  NITINOL  fibers  and, 
in  turn,  the  rate  at  which  the  corresponding  phase  transformation  forces 
are  recovered.  Changing  the  controller  gain  is  found  to  influence  to  some 
extent  the  critical  buckling  load.  For  example,  when  the  controller  gain 
is  2727  volt/mm  the  critical  buckling  load  is  950  N  and  when  the  gain  is 
reduced  to  136.35  volt/mm  the  critical  buckling  load  uecome  about  850N. 
Therefore,  reducing  the  controller  gain  by  a  factor  of  1/20  only  results  In 
a  10.5%  reduction  in  the  critical  buckling  load. 

The  effect  that  the  pre-tension  has  on  the  system  performance  is  shown 
in  Figure  (8)  for  a  dead  band  of  0.0176  mm  and  controller  gain  of  2727 
volts/mm.  It  is  clear  from  the  results  obtained  that  ,  the  pre-tension 
plays  the  most  crucial  role  in  controlling  the  buckling  of  the  beam. 
Increasing  the  tension  from  0.0  N/fiber  to  33.7  N/fiber  increases  the 
critical  buckling  load  from  about  100  N  to  950  N,  respectively. 

4.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENTS 

The  mechanism  of  actively  controlling  the  buckling  of  the 
NITINOL-reinforced  beam  can  best  be  understood  by  considering  Figure 
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Figure  (9-a)  *  Theoretical  prediction  of  critical  buckling  load  of 

controlled  and  uncontrolled  beams. 


Figure  (9-b)  -  Comparison  between  theoretical  and  experimental 

critical  buckling  loads  as  function  of  the  Initial 
fiber  pre-tension. 


ana  Auupiive  structures 


5.  CONCLUSIONS 

The  buckling  charact^eristics  of  NITINOL-reinforced  composite  beams  have 
been  presented.  The  fundatr«ntal  issues  governing  the  behavior  of  this  new 
class  of  SMART  composites  have  been  introduced.  Applications  of  NITINOL 
reinforcing  fibers  in  the  control  of  buckling  have  been  successfully 
demonstrated. 

Emphasis  has  been  placed  in  the  presentation  on  the  actuation 
capabilities  of  the  NITINOL  fibers.  Extensive  efforts  are,  however,  in 
progress  to  use  the  NITINOL  fibers  to  extract  modal  and  physical 
displacements  of  structures  with  multi-modes  of  vibration  (Baz,  Poh  and 
Gilheany  1991). 

With  such  built-in  sensing  and  controlling  capabilities. 
NITINOL-reinforced  composites  can  provide  a  means  for  continuously  tuning 
their  structural  characteristics  to  adapt  to  changes  in  the  operating 
conditions.  These  features  will  be  particularly  useful  in  many  critical 
structures  that  are  intended  to  operate  autonomously  for  long  durations  in 
isolated  environment  such  as  defense  vehicles,  space  structures  and 
satellites. 
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ABSTRACT 

A  new  class  of  distributed  sensors  is  presented  which  can  measure 
both  the  modal  and  physical  displacements  of  vibrating  composite  beams. 
The  sensor  relies  in  its  operation  on  a  set  of  super-elastic  Shape  Memory 
Alloy  (SMA)  wires  which  are  embedded  off  the  neutral  axes  of  the  vibrating 
beams.  The  wires  are  arranged  in  a  special  manner  which  allows  continuous 
monitoring  of  the  deflection  curve  of  the  beam.  The  output  signals  of  the 
SMA  wires  are  processed  to  determine  the  modal  displacements  of  the  beam 
and  the  physical  displacements  at  critical  discrete  points  along  the  beam 
axis. 

The  theoretical  and  experimental  performance  of  the  sensor  are 
presented  in  both  the  time  and  frequency  domains.  Comparisons  are  given 
between  the  experimental  performance  of  the  SMA  distributed  sensor  and 
that  of  conventional  laser  sensors  in  order  to  demonstrate  the  accuracy 
and  merits  of  the  distributed  SMA  sensor.  The  results  obtained  suggest 
the  potential  of  this  new  class  of  sensors  as  a  viable  means  for 
controlling  the  vibrations  of  flexible  composite  beams  and  plates 
particularly  with  modal  control  algorithms.  [  Work  supported  by  grant  from 


ARO]. 


1.  INTOODUCTION 


Considerable  interest  has  been  directed  recently  towards  the 
development  of  a  wide  variety  of  distributed  sensors  for  monitoring  the 
vibration  of  composite  flexible  structures.  Distinct  among  these  sensors 
are  those  relying  in  their  operation  on  optical  fibers  [1-2], 
piezo-electric  films  [3-4]  and  Shape  Memory  Alloy  (SMA)  wires  [5-6]  which 
are  embedded  in  the  composite  matrix  of  the  flexible  structures.  By  far 
the  most  commonly  used  and  researched  class  of  embedded  sensors  is  the 
optical  fibers  sensors.  These  sensors  have  been  developed  to  monitor  not 
only  structural  vibrations  [1]  but  also  temperature  distribution  [7], 
cracks  and  defects  [8]  as  well  as  curing  of  composite  matrices  [9].  In 
spite  of  such  wide  acceptance,  the  complexity  of  the  instrumentation  and 
the  signal  processing  algorithms  associated  with  optical  fibers  sensors 
still  poses  many  serious  challenges  that  remain  to  be  addressed  [10]. 

Such  problems  are  avoided  with  the  embedded  piezo-electric  film 
sensors  developed  by  Lee  [3].  These  sensors  utilize  specially  shaped 
piezo-electric  films  to  monitor  the  modal  coordinates  of  vibrating  beams 
and  plates.  The  shaping  of  the  films  aims  at  isolating  the  individual 
modal  signals  in  a  manner  similar  to  the  modal  filtering  technique  of 
Meirovitch  [11].  Such  specially  configured  piezo-electric  films  have  also 
been  used  by  Lee  et  al  [4]  to  monitor  uniaxial  and  pure  shear  strain  rates 
in  vibrating  beams.  In  spite  of  their  on-line  and  real  time  measurement 
capability  without  the  need  for  any  signal  processing,  the  configured 
piezo-electric  films  are  restricted  to  the  singular  structural 
configuration  that  they  are  designed  for  and  are  limited  to  simple  and 
time  invarying  structures. 

The  use  of  shape  memory  wires  as  embedded  distributed  sensors  provides 


a  viable  alternative  to  optical  fibers  and  piezo-electric  film  sensors. 
The  SMA  sensors  require  very  simple  instrumentation  and  when  properly 
designed,  they  require  very  simple  signal  processing  algorithms  to  extract 
modal  and  physical  coordinates  of  the  vibrating  structures.  The  sensors 
rely  in  their  operation  on  the  changes  in  the  electrical  resistance  of  the 
SMA  wires  when  subjected  to  a'dal  loading  in  a  manner  similar  to 
conventional  discrete  strain  gage  sensors.  This  physical  property  has 
been  successfully  utilized  in  monitoring  the  position  of  robotic  joints 
[12].  The  SMA  sensor  have  also  been  used  at  Virginia  Polytechnic 
Institute  [5]  and  Boeing  Aerospace  [6]  to  measure  a  single  mode  of 
vibration  of  composite  flexible  beams  [5-6].  No  attempt  has  been  made 
however  to  use  the  SMA  sensor  to  monitor  several  modes  of  vibration  or  to 
extract  the  modal  and  physical  parameters  of  vibrating  structures. 

It  is  therefore  the  purpose  of  this  study  to  develop  such  a  SMA  sensor 
which  is  capable  of  monitoring  multi-modes  of  vibration  of  flexible 
composite  beams  in  order  to  extract  individual  modal  and  physical 
parameters  of  these  beams.  Furthermore,  this  study  documents  the 
performance  of  such  SMA  sensors,  both  in  the  time  and  frequency  domains, 
and  validates  their  performance  by  comparison  with  conventional  sensors. 
It  is  important  to  note  that  the  SMA  sensor  relies  in  its  operation  on  the 
use  of  the  classical  modal  decomposition  approach  [13].  This  approach  has 
been  recently  utilized,  by  many  investigators,  to  extract  the  modal 
displacements  of  vibrating  beams  [14-16]  and  plates  [17]  from  the 
measurements  of  a  set  of  discrete  sensors.  The  placement  strategy  of  these 
discrete  sensors  is,  however,  very  critical  in  order  to  avoid  nodes  of 
vibration.  This  is  in  contrast  to  the  distributed  SMA  sensor  which  can  be 
placed  at  any  place  along  the  beam  without  any  concern  about  the  location 
of  the  nodes  of  vibration.  Such  versatility  makes  the  SMA  sensor  suitable 


for  monitoring  the  vibration  of  more  complex  structures  where  the  location 
of  the  nodes  of  vibration  can  not  be  accurately  determined  particularly 
for  higher  order  modes. 

Accordingly,  this  paper  is  organized  in  four  sections  and  an  appendix. 
In  section  1,  a  brief  introduction  is  given  and  section  2  presents  the 
theoretical  principles  of  the  sensor.  The  experimental  performance  of  the 
sensor  is  given  in  section  3  and  the  conclusions  are  summarized  in 
section  4.  The  calibration  of  the  sensor  is  described  in  the  appendix. 

2.  THEORETICAL  MODELING  OF  THE  SENSOR 

2.1.  Monitoring  the  modal  coordinates 

The  multi-mode  distributed  SMA  sensor  relies  in  its  operation  on  a  SMA 
wire  embedded  inside  a  composite  beam  as  shown  in  Figure  (1).  The  wire  is 
placed,  at  a  distance  a,  off  the  neutral  axis  of  the  beam.  The  wire  is 

tapped  along  its  length  at  different  locations  1,  2,  ..  ,  N  where  N  is 

equal  to  both  the  number  of  modal  coordinates  to  be  extracted  and  the 
number  of  physical  coordinates  to  be  estimated.  The  changes  AL^,  AL2, 
..and  AL^  in  the  lengths  L^,  L2,  ..  and  Lj,  of  the  different  segments  of 

the  wire,  due  to  the  vibration  of  the  beam,  are  monitored  continuously  by 

measuring  the  changes  AR^,  AR2,  ..  and  AR^  in  the  electrical  resistances 
of  these  segments  such  that 

ALj  =  (ARi  A)  /  for  i=l, . . ,  N  (1) 

where  A  is  the  wire  cross  sectionai  area  and  is  its  electrical 
resistivity. 

The  changes  in  the  lengths  ALj  are  related  to  the  curvature  of  the 
deflection  curve  of  the  beam  as  it  vibrates  under  the  action  of  the 
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Figure  (1)  -  Schematic  drawing  of  the  SMA  distributed  sensor 


external  loads.  Assuming  small  deflections  and  that  the  SMA  wire 
experiences  the  same  linear  strains  as  the  host  beam  in  which  the  wire  is 
embedded,  then  the  changes  ALj  can  be  calculated  from 


AL,  = 


0 


e  dx 


dx 


(2) 


where  e  is  the  axial  strain  in  the  wire  and  y  is  the  beam  deflection 
in  the  transverse  direction  which  obeys  the  following  partial  differential 
equation  describing  the  dynamics  of  a  Bernoulli -Euler  beam  [13],  is 


3*y  3^y 

E  I  -  +  p  Ah  -  =  y  5(x-Xi)  +  y  M,  5’(x-x,)  (3) 

X  d  t 

where  El  is  the  flexural  rigidity  of  the  beam,  p  is  the  mass  density, 
Afj  is  the  beam  cross  sectional  area,  t  is  the  time,  Fj  is  a  point  load 
acting  on  the  beam  at  location  Xj,  Mj  is  a  moment  acting  on  the  beam  at 
location  Xj,  8  is  the  kronecker  delta  function  and  S'  is  the  unit  doublet 
function.  Applying  the  separation  principle  [13]  to  equation  (3),  the 
transverse  deflection  y(x, t),  at  any  location  x  and  time  t,  can  be  written 
as  a  linear  combination  of  the  mode  shapes  0j(x)  of  the  beam  as  follows 

N 

y(x,t)  =  <^i(x)  UjCt)  (4) 

1  =  1 


where  Uj(t)  and  0i(x)  are  the  generalized  modal  coordinate  and  the 
corresponding  mode  shape  respectively  for  the  ith  mode.  The  mode  shape 
0j(x)  for  a  cantilevered  beam  is  given  by 


0i(x)  =  [coshCkjx]  -  cos(kjx)]-ai  [sinh(kjx)-sin(kjx) ]  (5) 


with  OLy  given  by 
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ttj  =  [cosh(kjLi)  +  cosCkjLj)]  /  [sinh(kjLi )+sin(kiLi ) ]  (6) 

and  kj  is  defined  as  the  wave  number  of  the  ith  mode  which  satisfies  the 
following  characteristic  equation 

cos  (kjLj)  cosh  (kjLj)  +1=0  (7) 

Substituting  equations  (5),  (6)  and  (7)  into  equation  (4)  the 

transverse  deflection  y  at  any  x  and  t  can  be  determined.  The  resulting 

2  2 

equation  can  be  used  to  calculate  the  beam  curvature  3  y/3x  which  is  in 
turn  can  be  substituted  into  equation  (2)  to  calculate  the  change  in 
■lengths  ALj. 

After  some  manipulations,  the  vector  of  length  changes  ALj  can  be 
calculated  from 


ALj  “  ^  ^  kj  ^|sinh(kjLt  )+sin(kjL| )  -  aj  |cosh(kjLj  )-cos(kjLj )  j  j Uj  (8) 

J=i 

Letting  AL  be  the  vector  of  length  changes  given  by 


AL  =  [ALj  AL2  AL^]^, 


(9) 


and 


C  be  a  (NxN)  matrix  whose  constant  and  known  entries  Cjj  are 
given  by 

Cjj  =a  kj ^|sinh(kjLj )+sin(kjLj ) j  -  aj |cosh(kjLj )-cos (kjLj ) j j ,  (10) 

U  be  the  vector  of  the  beam  modal  coordinates  given  by 

U  =  [  Uj  U2  .  .  u^ )  ’*'  (11) 


Then  substituting  equations  (8),  (10)  and  (11)  into  equation  (9),  it 

yields  the  following  matrix  equation  for  the  vector  of  length  changes 


AL  =  C  U 


(12) 
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indicating  that  the  changes  in  the  lengths  AL  of  the  different  segments  of 
the  SMA  wire  are  linearly  related  to  the  modal  coordinates  U  of  the 
composite  beam. 

Defining  AR  as  the  vector  of  resistance  changes  of  the  different 
wire  segments  ,  i.e. 

AR  =  EARi  AR2  ...  AR^]”^  (13) 

and  substituting  for  ARj  from  equation  (1),  the  above  equation  reduces  to 

AR  =  A/Pe  AL  (14) 

Combining  equations  (12)  and  (14),  gives  the  modal  coordinate  vector  U 
in  terms  of  the  resistance  changes  vector  AR  as  follows 

U  =  Pg/A  AR  (15) 

i.e.  measuring  the  changes  AR  of  the  N  segments  of  the  SMA  wire  can  be 
used  to  directly  calculate  the  N  modal  coordinates  of  the  vibrating 
composite  beam  as  the  matrix  equation  (15)  represents  a  set  of  N  linear 
equations  in  the  N  unknowns  u^,  U2,  ..  and  u^. 

It  should  be  noted  that  the  matrix  C  is  always  non-singular  i.e.  a 
solution  for  equation  (15)  always  exists,  unless  two  tapping  points  of  the 
SMA  wire  coincide.  This  means  that  the  length  of  the  corresponding  two 
wire  segments  is  equal  resulting  in  a  matrix  C  with  two  equal  rows. 
Accordingly,  the  rank  of  matrix  C  is  reduced  to  N-1.  Of  course  such 
condition  will  not  occur  on  purpose.  However,  if  it  does  occur  the  sensor 
can  still  extract  (N-1)  modal  and  (N-1)  physical  displacements.  In  other 
words,  the  tapping  points  of  this  distributed  sensor  can  be  placed  any 
where  along  the  beam  even  coinciding  with  any  node  of  vibration.  This  is 
not  the  case  with  discrete  sensors  such  as  conventional  strain  gages  and 
accelerometers  where  these  sensors  should  not  be  placed  at  nodes  of 
vibration.  Such  critical  placement  requirement  of  the  discrete  sensors 
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limits  their  use  to  simple  structures  and  favors  the  use  of  the 


distributed  sensor  described  in  this  study. 

2.2.  Monitoring  the  physical  coordinates 

In  a  similar  fashion,  the  measurements  of  AR  can  also  be  used  to 
compute  the  N  physical  coordinates  of  the  vibrating  beam.  This  can  be 
achieved  by  defining  Y  as  a  vector  of  physical  coordinates  given  by 

Y  =  [  y(Li)  ylLg)  ....  y(L(,)]^  (16) 


Therefore,  from  equation  (4)  the  following  matrix  equation  can  be 
written 


Y 

i  U 

(17) 

where  $ 

= 

r  ^i(Li)  . 

•  1 

(18) 

(pxiLz) 

0j(Ljj) 

- 

From  equations  (15) 

and 

(17),  the 

physical  transverse 

displacement 

vector  Y  of  the  N  taping  points  can  be  calculated  from  the  measurements  of 
the  resistance  changes  AR  of  the  N  segments  of  the  SMA  wire  as  follows 

Y  =  pg  /  A  41  C'^  AR  (19) 

Equations  (15)  and  (19)  constitute  the  basic  equations  that  govern 
describing  the  capability  of  the  distributed  SMA  sensor  in  extracting  the 
modal  and  physical  coordinates  of  a  vibrating  composite  beam  from  the 
measurements  of  resistance  changes  of  N  segments  of  the  wire.  In  these 
equations,  the  calculation  of  both  the  modal  and  physical  coordinates  of  a 
beam  is  shown  to  be  very  simple  as  it  is  reduced  to  merely  the  solution  of 


-  9  - 


two  sets  of  simultaneous  linear  equations. 

Equations  (15)  and  (19)  will  be  utilized  in  the  following  section  to 
extract  the  coordinates  of  an  experimental  beam  and  the  results  obtained 
will  be  compared  to  the  results  measured  by  conventional  laser  sensors. 

The  above  analysis  of  the  SHA  sensor  indicates  that  the  sensor  is  in 
effect  a  multi-mode  distributed  strain  gage.  Such  distributed  nature  of 
the  sensor  has  many  inherent  advantages.  First,  it  makes  the  placement  of 
the  sensor  and  its  tapping  points  insensitive  to  the  location  of  the  nodes 
of  vibration.  This  of  course  is  not  the  case  for  conventional  discrete 
strain  gages  where  the  sensors  can  not  be  placed  at  these  nodes.  Second, 
because  the  SMA  sensor  relies  on  integrating  the  strain  along  the  wire 
segments,  as  described  by  equation  (2),  its  output  signal  will  be  less 
noisy  than  that  of  conventional  strain  gages  and  high  signal-to-noise 
ratios  can  be  obtained.  Third  and  because  of  its  embedded  nature  the  SMA 
sensor  can  also  detect  structural  failures  by  monitoring  the  failure  of 
any  wire  segment.  These  advantages  make  the  SMA  sensor  suitable  for 
accurately  monitoring  the  vibration  and  integrity  of  complex  SMART 
composites. 

Worth  noting  also  is  that  the  theory  developed  for  the  SMA  sensor 
is  based  on  the  knowledge  of  the  mode  shapes  of  the  beam  or  the  structure 
inside  which  it  is  embedded.  However,  the  cantilevered  beam  example  given 
in  this  study  can  be  easily  extended  to  beams  with  other  boundary 
conditions. 


3.  EXPERIMENTAL  PERFORMANCE  OF  THE  DISTRIBUTED  SENSOR 
3.1.  Experimental  set-up 


Figure  (2)  shews  a  schematic  layout  of  the  experimental  set  up  used 


K 


Figure  (2)  -  Schematic  drawing  of  the  experimental  set-up  used 

in  testing  the  SMA  sensor. 


in  evaluating  the  performance  of  the  SMA  distributed  sensor  in  both  the 
time  and  frequency  domain.  The  sensor  is  made  of  a  0. 125  mm  diameter  SMA 
wire  embedded  inside  a  composite  beam  whose  physical,  geometrical  and 
dynamical  properties  are  given  in  Table  1.  The  SMA  wire  itself  is 
manufactured  from  a  super-elastic  Nickel-Titanium  alloy  called  NITINOL 
which  has  a  resistivity  of  0.644x10“^  ohm/m.  Such  an  alloy  is  selected 
because  its  super-elasticity  allows  the  wire  to  experience  large 
deformations,  about  ten  times  more  than  conventional  materials,  and  still 
completely  spring  back  to  its  original  undeformed  shape  [18].  This 
feature  makes  the  sensor  suitable  for  monitoring  the  vibration  of  very 
flexible  structures  without  exhibiting  any  plastic  deformations. 


Table  1  -  Physical,  geometrical  and  dynamical  properties  of 

test  composite  beam 


Length  Width 

thickness 

density 

Young’ s  Mod. 

1st  Mode 

2nd  Mode 

3rd  Mode 

(cm) 

(cm) 

(cm) 

(gm/cc) 

(GN/m  ) 

(Hz) 

(Hz) 

(Hz) 

48.0 

5.0 

0.25 

0.73 

4.35 

1.88 

11. 75 

33.  15 

The  composite  beam,  with  the  embedded  SMA  sensor,  is  mounted  in  a 
cantilevered  manner  on  an  oscillating  platform  which  can  oscillate  freely 
on  two  guide  rails.  The  rails  are  set  parallel  to  the  transverse 
vibration  direction  of  the  beam.  A  shaker  (Wilkoxon  Research  model  F3/F9, 
Bethesda,  MD)  is  driven  with  a  sinusoidal  function  generator  to  provide 
controlled  vibration  of  the  platform  and  in  turn  the  composite  beam.  The 
resulting  beam  vibration  is  monitored  by  the  SMA  sensor  which  is  provided 
with  only  two  segments.  The  first  segment  extends  between  the  beam  tip  1 
and  its  fixed  root  whereas  .e  second  segment  extends  between  the  beam 
mid-span  point  2  and  the  fixed  base.  The  resistances  of  the  two  segments 
are  measured  using  amplifiers  model  P-3500  from  Measurements  Group, 
Raleigh,  North  Carolina.  The  amplifiers  operate  with  9  Volt  DC  power 


source.  The  conditioned  outputs  of  the  amplifiers  are  sampled  by  a 
386-microprocessor  via  an  input/output  board  (model  DASH-16  from 
METRABYTE,  Taunton,  MA).  The  sampling  rate  used  is  190  Hz  with  an 
accuracy  of  12  bits  which  corresponds  to  4.88x10“^  mm.  The  displacements 
of  the  first  and  second  modes  are  computed  from  the  measurements  of  the 
resistance  of  the  two  segments  using  equation  (15).  Also,  the  physical 
displacements  of  the  beam  tip  and  mid-span  are  computed,  in  the 
microprocessor,  using  equation  (19).  The  physical  displacements  obtained 
from  the  SMA  sensor  are  compared  with  those  measured  by  two  laser  sensors 
(model  IIIB-LA40HR  ,  Aromat  Corp. ,  New  Providence  ,  NJ).  The  laser 
sensors  are  placed  on  the  oscillating  platform  facing  the  beam  tip  and 
mid-span  points.  The  sensors  have  accuracy  of  20  pm  over  a  frequency 
band  between  0-1000  Hz. 

The  calibration  of  the  SMA  sensor  is  carried  out  by  comparing  the 
sensor  output  with  the  output  of  the  laser  sensors  according  to  the 
procedure  outlined  in  the  appendix.  Such  a  procedure  is  important  to 
account  for  errors  in  the  sensor’s  parameters  as  the  fiber  position  a  and 
segments  length  Lj®.  The  laser  sensors  are  calibrated  by  reflecting  their 
light  beam  off  a  micrometer  head.  The  spacing  between  the  head  and  the 
laser  emitter  is  varied  and  the  output  of  the  laser  sensor  is  recorded. 

3.2.  Performance  of  the  SMA  sensor 
3.2.1.  in  the  time  domain 

The  composite  beam  with  the  embedded  SMA  sensor  is 
subjected  to  three  different  initial  sinusoidal  excitations  at  1.9,  7  and 
11.8  Hz  respectively.  The  first  and  third  excitation  frequencies  are 
selected  in  order  to  resonant  the  beam  at  its  first  and  second  modes  of 


vibration.  Following  the  initial  excitation  period,  the  shaker  driving  the 


platform  is  turned  off  and  the  beam  is  left  to  vibrate  freely. 


The 


signals  of  the  SMA  sensor  and  the  laser  sensors  are  then  sampled  by  the 
computer  over  a  period  of  1.5  seconds.  All  the  sampling  and  the  results, 
reported  in  this  study,  are  carried  out  during  the  transient  period  which 
follows  the  steady-state  excitation  by  the  shaker. 

Figure  (3)  shows  the  corresponding  length  changes  of  the  two  SMA 
wire  segments  at  the  three  excitation  conditions.  These  length  changes  AL 
are  used  along  with  equation  (19)  to  compute  the  physical  displacements  at 
the  tip  and  mid-span  of  the  beam.  Figure  (4)  shows  the  computed  physical 
displacements  as  obtained  by  the  SMA  sensor  along  with  the  corresponding 
displacements  of  the  laser  sensors  when  the  first  vibration  mode  of  the 
beam  is  excited.  Such  a  comparison  indicate  excellent  agreement  between 
the  measurements  of  the  SMA  and  the  laser  sensors. 

Figures  (5)  and  (6)  display  similar  comparisons  between  the 
displacements  of  the  SMA  and  the  laser  sensors,  at  the  beam  tip  and 
mid-span,  when  the  excitation  frequencies  are  7  and  11.8  Hz  respectively. 
The  two  figures  clearly  demonstrate  the  accuracy  of  the  SMA  in  monitoring 
the  vibration  of  the  beam  in  the  time  domain. 

It  is  essential  to  note  that  all  the  measurements  and  the 
calculations  of  the  modal  and  physical  displacements  are  carried  out,  in 
real  time,  with  a  sample  interval  of  5.2  ms.  This  interval  is 
corresponding  to  a  sampling  rate  of  190  Hz  which  is  at  least  8  times 
greater  than  the  Nyquist  frequency  limit  necessary  to  avoid  aliasing. 

3.2.2.  in  the  frequency  domain 

The  performance  of  the  SMA  in  the  frequency  domain  is 
evaluated  by  impacting  the  composite  beam  with  a  modally  tuned  impact 
hammer  (model  PCB-  GK291B02,  PCB  Piezotronics,  Depew,  NY).  The  response 
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Figure  (5)  -  Comparison  between  displacement  measurements  of  SMA 
and  laser  sensors  at  tip  and  mid-span  of  beam  when 
subjected  to  sinusoidal  exclatation  of  7.0  Hz. 


Mid-Span  Displacement  -  yg  (mm)  Tip  Displacement  -y  ,  (mm) 


Figure  (6)  -  Comparison  between  displacement  measurements  of  SMA 
^  and  laser  sensors  at  tip  and  mid-span  of  beam  when 

subjected  to  sinusoidal  exciatation  of  11.8  Hz. 


of  the  SMA  sensor  as  well  as  that  of  the  laser  sensors  to  such  impacts  are 
analyzed  using  a  dual  channel  spectrum  analyzer.  These  responses  are  used 
to  predict  the  frequency  response  of  the  different  vibration  modes. 
Figures  (7)  and  (8)  show  the  normalized  power  spectra  of  the  SMA  and  the 
laser  sensors  outputs  at  the  beam  tip  and  its  mid-span  respectively.  It 
is  evident  that  the  response  of  the  SMA  sensor  matches  that  of  the  laser 
sensors  in  the  frequency  domain.  Furthermore,  the  experimental 
frequencies  obtained  by  the  analysis  of  the  SMA  sensor  output  are  in  close 
agreement  with  the  theoretical  predictions  Wj  which  are  given  by 

1/2 

=  (ki  Li)2/  (Li^)  (EI/pA)  (20) 

where  kj  is  the  wave  number  of  the  ith  mode  given  by  equation  (7).  In 
particular,  the  first  three  modes  are  found  experimentally  to  be  1.875, 
11.76  and  33.15  Hz  whereas  the  corresponding  theoretical  predictions  are 
1.89,  11.85  and  33.18  Hz  respectively. 

However,  it  is  essential  to  note  that  the  bandwidth  of  the 
prototype  SMA  sensor  is  limited  by  the  fact  that  it  has  two  sampling 
points.  This  implies  that  it  can  only  measure  excitations  up  to  the 
second  mode  of  vibration  which  is  11.8  Hz. 


4.  CONCLUSIONS 

The  feasibility  of  a  new  class  of  distributed  SMA  sensors  in 
measuring  modal  and  physical  displacements  of  composite  beams  is 
successfully  demonstrated.  The  sensors  utilize  SMA  wires  which  are 
embedded  inside  the  matrix  of  the  composite  beams.  Such  an  arrangement 
can  be  augmented  with  sets  of  embedded  actuators,  as  SMA  or 
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Figure  (7)  -  Comparison  between  frequency  response  of  SMA  and  the 
laser  sensors  at  tip  of  beam  when  impacted  with  an 
impact  hammer. 
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Figure  (8)  -  Comparison  between  frequency  response  of  SMA  and  the 
laser  sensors  at  mid-span  of  beam  when  impacted  with 
an  impact  hammer. 


piezo-electric,  in  order  to  develop  a  class  of  SMART  composite  beams  which 
have  built-in  vibration  and  shape  control  capabilities.  These  capabilities 
can  be  utilized  along  with  any  of  the  modal  control  algorithms,  such  as 
those  developed  by  Baz  et  al  [15-16],  to  control  the  structural  vibrations 
and  shape  of  these  SMART  composites. 

The  general  theory  of  operation  of  the  distributed  sensor  is 
presented  for  monitoring  N  modal  and  physical  coordinates.  The 
experimental  validation  of  the  sensor  performance  is  demonstrated  for  a 
composite  beam  where  the  first  two  modes  of  vibrations  dominate  its 
dynamic  response.  The  physical  displacements  of  the  beam  tip  and  mid-span 
as  obtained  by  the  two-mode  distributed  sensor  are  found  to  be  in 
excellent  agreement  with  the  measurements  of  conventional  laser  sensors. 
Emphasis  however  should  be  placed  on  the  fact  that  the  distributed  sensor 
is  tested  with  excitations  within  its  frequency  band  capabilities.  If  the 
sensor  is  designed  to  monitor  low  frequency  excitations  and  the  external 
frequency  band  exceeds  its  maximum  frequency  limit,  then  it  is  essential 
to  use  low  pass  filters  to  filter  out  the  high  frequency  contents  in  the 
sensor  signals  in  order  to  avoid  conventional  observation  spillover 
problems  [19]  and  spatial  aliasing  [17].  This  typical  instrumentation 
problem  can  be  avoided  by  proper  selection  of  the  sensor’s  bandwidth 
to  be  compatible  with  the  expected  excitation  frequency  band.  However, 
for  more  accurate  monitoring  of  the  displacements  it  is  recommended  to 
provide  the  SMA  sensor  with  more  bandwidth  than  the  excitation  frequency 
band,  i.e  number  of  sampling  points  N  should  be  greater  than  the  number  of 
displacements  P  to  be  determined.  In  this  manner,  equation  (15)  becomes: 

U  =  [uj  U2  ..  Up]^  =  p,/A  [C^  C]-^  AR  (21) 

where  U  is  the  vector  of  P  modal  displacements,  C  is  a  matrix  (PxN)  whose 
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entries  are  given  by  equation  (10)  and  C]"^  is  the  pseudo-inverse 
of  C.  Such  a  statistical  approach,  as  N>P,  of  the  over  determined  equation 
C  U  =  Pg/A  AR  ensures  that  the  modal  displacements  U  extracted  are  the 
least-squares  estimates.  This  will  be  particularly  important  in  the 
case  of  noisy  measurements.  The  modal  displacements  U  extracted  can  then 
be  used  to  compute  the  physical  displacements  Y  of  the  beam  using  the 
following  transformation  equation 

Y  =  i  U  (22) 

where  $  is  the  mode  shape  matrix  (NxP)  whose  entries  correspond  to  the 
shapes  of  the  jth  mode  at  location  i. 

It  should  be  pointed  out  also  that  the  theory  developed  in  this 
study  can  be  easily  extended  to  utilize  the  distributed  SMA  sensor  to 
monitor  the  vibration  of  two  or  three  dimensional  motions  as  for  example 
of  plates  or  more  complex  structures. 
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APPENDIX 


CALIBRATION  OF  THE  SNA  DISTRIBUTED  SENSOR 

The  linear  nature  of  the  SMA  distributed  sensor,  as  defined  by 
equations  (15)  and  (19),  makes  its  calibration  a  rather  simple  process. 
Such  a  calibration  is  necessary  to  account  for  errors  in  the  parameters  of 
the  sensor  such  as  the  position  and  length  of  the  SMA  wire  and  its 
different  segments.  The  calibration  procedure  is  based  on  the  use  of  a 
conventional  laser  sensor  placed  at  location  i  along  the  beam  to  monitor 
the  beam  physical  displacement  at  time  sample  k.  The  beam  is 

subjected  to  a  general  excitation  and  the  corresponding  resistance  changes 

AR^  of  the  different  wire  segments  are  determined  at  the  same  sample 

interval  k.  The  collection  of  a  set  of  z  measurements  of  the  physical 
displacement  yjj^  and  the  resistance  changes  AR^  may  be  concatenated,  based 
on  equation  (19),  to  form  the  following  equation; 

Y,  =  A??  Cl  (A-1) 

where 

Yi  =  vector  of  concatenated  physical  displacement  measurements 

given  by 

T 

“  j^YKk)  Ynk+i)  ••••  Yuk+z-i)  j  ’  (A-2) 

A3?  =  matrix  of  concatenated  resistance  changes,  given  by 


■^^l(k) 

^f^2(k) 

.  .  .  .  AR^dc) 

- 

^^l(k+l) 

^2(k+l) 

.  ...  AR„(,j+l, 

(A-3) 

^^^Kk+z-l)  ^f^2(k+z-l)  ••••  ^^M(k+z-l) 

and  Cj  =  vector  of  sensor  parameters  at  location  i,  given  by 
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'12 


(A-4) 


The  elements  of  the  vector  Cj  can  then  be  computed  using  the  least 
mean  squared  estimates  as  follows 

Cl  =  [A3?^  Yj  (A-5) 


The  above  procedure  is  repeated  by  moving  the  laser  sensor  to 
location  j  and  estimating  the  vector  Cj  of  sensor  parameters  in  a  similar 
manner.  When  all  the  rows  Cj®  are  estimated,  the  matrix  C  is  formed  and 
the  calibration  process  of  the  distributed  sensor  is  completed  provided 
that  the  modal  shape  matrix  4  is  already  known  or  is  experimentally 
determined  using  the  classical  modal  analysis  methods  [20]. 
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NOMENCLATURE 


Latin  letters 


A 

C 

C 

Cl 

E 

F 

I 

k 

Li 

M 

N 

P 

Ri 

t 

T 

U 

U 

X 

y 

Y 

Yi 

z 

Greek  letters 

“i 

ALi 

AL 

ARj 

AR 

A3? 

e 

P 

Pe 

<Pi 

0 


distance  between  SMA  wire  and  neutral  axis 

cross  sectional  area  of  wire 

cross  sectional  area  of  beam 

matrix  given  by  equation  (10)  (NxN) 

matrix  given  in  equation  (21)  (PxN) 

ith  row  of  C  matrix 

Young’s  modulus  of  composite  beam 

force  acting  on  beam 

area  moment  of  inertia  of  beam 

sample  order 

the  wave  number  for  the  ith  mode  given  by  equation  (7) 
length  of  ith  segment 
moment  acting  on  beam 

number  of  wire  segments,  number  of  measured  modes, 
and  number  of  measured  physical  coordinates 
number  of  accurate  modal  displacements 
resistance  of  ith  wire  segment 
time 

beam  thickness 
modal  coordinate  of  ith  mode 
modal  coordinate  vector  (Nxl) 
accurate  modal  coordinate  vector  (Pxl) 
distance  along  the  neutral  axis  of  beam 
deflection  in  transverse  direction 
vector  of  physical  deflection  of  tapping  points 
vector  of  concatenated  physical  displacements  at  ith 
location 

number  of  calibration  samples 


constant  for  ith  mode  given  by  equation  (6) 
change  in  length  of  ith  segment 
vector  of  length  changes  (Nxl) 
changes  in  resistance  of  ith  segment 
vector  of  resistance  changes  (Nxl ) 
matrix  of  concatenated  resistance  changes 
strain  in  SMA  wire 
mass  density  of  composite  beam 
electrical  resistivity  of  SMA  wire 
mode  shape  for  ith  mode 

modal  shape  matrix  (NxN)  given  in  equation  (17) 
modal  shape  matrix  (NxP)  given  in  equation  (22) 
the  ith  mode 
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ABSTRACT 

The  static  and  dynamic  characteristics  of  flexible  fiberglass 
composite  beams  are  controlled  by  activating  optimal  sets  of  shape  memory 
alloy  (NITINOL)  wires  which  are  embedded  along  the  neutral  axes  of  these 
beams.  The  underling  phenomena  influencing  the  behavior  of  this  class  of 
composite  structural  members  are  presented.  The  individual  contributions 
of  the  fiberglass-resin  laminate,  the  NITINOL  wires  and  the  shape  memory 
effect  to  the  overall  performance  of  the  composite  beam  are  determined  at 
different  operating  temperatures  and  initial  preloads  of  the  wires.  The 
modes  of  vibration  of  the  fiberglass  beams  are  measured  with  and  without 
the  NITINOL  reinforcement  at  various  operating  conditions.  With  properly 
designed  NITINOL  reinforcement,  it  is  shown  that  the  beams  can  become 
stiffen  and  less  susceptible  to  buckling.  The  modes  of  vibrations  of  the 
activated  Nitinol-reinforced  composite  beams  can  also  be  shifted  to  higher 
frequency  bands  relative  to  those  of  the  unactivated  or  un-reinforced 
beams.  Finite  element  models  are  developed  to  describe  the  static,  dynamic 
and  thermal  interaction  between  the  NITINOL  wires  and  the  fiberglass-resin 
laminates.  Close  agreement  is  obtained  between  theoretical  predictions  and 
experimental  results.  With  such  tunable  characteristics,  the 
NITINOL-reinforced  composite  beams  can  be  effective  in  attenuating  the 
vibrations  induced  by  various  external  disturbances. 


1.  INTRODUCTION 


Considerable  attention  has  recently  been  devoted  to  the  utilization  of 
the  Shape  Memory  Nickel-Titanium  alloy  (NITINOL)  in  developing  SMART 
composites  that  are  capable  of  adapting  intelligently  to  external 
disturbances  (Ikegami  et  al.  [1],  Rogers  et  al.  [2],  and  Baz  et  al.  [3-4]). 
Such  wide  acceptance  of  NITINOL  stems  from  its  unique  behavior  when  it  is 
subjected  to  particular  heating  and  cooling  strategies.  For  example,  the 
alloy  becomes  soft  when  it  is  cooled  below  its  martensite  transformation 
temperature  and  becomes  about  four  times  stiffer  when  it  is  heated  above 
its  austenite  transformation  temperature  (Funakubo  [5]).  Furthermore,  when 
the  alloy  is  trained  to  have  a  particular  shape  while  in  its  austenite 
phase,  it  will  memorize  this  shape.  If  the  alloy  is  then  cooled  to  its 
martensite  phase  and  subject  to  plastic  deformation,  it  will  return  to  its 
memorized  shape  when  it  is  heated  above  the  austenite  transformation 
temperature.  The  phase  transformation  from  martensite  to  austenite  produces 
significant  forces  as  the  alloy  recovers  its  original  shape.  The  alloy 
acts  as  an  actuator  transforming  thermal  energy  into  mechanical  energy 
(Perkins  [6]  and  Duerig  et  al.  [7]).  Accordingly,  if  the  NITINOL  fibers 
are  embedded  inside  a  composite  matrix  at  optimal  locations,  they  can  be 
used  to  control  the  static  and  dynamic  characteristics  of  the  resulting 
SMART  composite.  The  control  action  is  generated  by  the  described 
stiffening  of  the  NITINOL  fibers  and/or  the  shape  memory  effect.  With  such 
built-in  control  capabilities,  the  performance  of  the  SMART  composites  can 
be  optimized  and  tailored  to  match  changes  in  operating  conditions. 

Emphasis  is  placed  in  the  present  work  on  using  the  shape  memory  effect 
of  the  NITINOL  fibers  to  control  the  performance  of  fiberglass  composite 
beams.  The  NITINOL  fibers  are  embedded  inside  vulcanized  rubber  sleeves 
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placed  along  the  neutral  axes  of  these  composite  beams  as  shown  in  Figure 
(1).  In  this  arrangement,  the  fibers  are  free  to  move  during  the  phase 
transformation  process  in  order  to  avoid  degradation  and/or  destruction  of 
the  shape  memory  effect  which  may  result  when  the  fibers  are  completely 
bonded  to  the  composite  matrix. 

The  basic  phenomena  governing  the  static  and  dynamic  performance  of  the 
NITINOL-reinf orced  composites  will  be  presented,  both  theoretically  and 
experimentally,  in  this  paper. 

The  paper  is  organized  in  five  sections.  In  section  1,  a  brief 
introduction  is  given.  In  sections  2  and  3,  the  finite  element  models 
describing  the  static  and  dynamic  characteristics  of  the  NITINOL-reinforced 
beams  are  introduced  respectively.  The  experimental  validation  of  the 
finite  element  models  is  presented  in  section  4  and  the  conclusions  are 
summarized  in  section  5. 

2.  STATIC  CHARACTERISTICS  OF  NITINOL-REINFORCED  BEAMS 

The  static  characteristics  of  NITINOL-reinforced  composite  beams  are 
primarily  governed  by  their  stiffness.  The  beam  stiffness  is  made  up  of 
different  components  which  include:  the  flexural  rigidity  of  the  beam,  the 
geometric  stiffness  that  accounts  for  the  axial  and  thermal  loading,  as 
well  as  the  stiffness  imparted  by  the  elasticity  of  the  NITINOL  fibers. 
The  individual  components  of  the  beam  stiffness  can  be  determined  by 
considering  the  NITINOL-reinforced  beam  element  shown  in  Figure  (2)  with 
the  forces  acting  on  it  and  the  associated  displacements.  The  combined 
stiffness  of  the  element  can  be  obtained  using  the  principle  of 
conservation  of  energy  and  equating  the  work  done  by  external  loads  to  the 
strain  energies  stored  in  the  element.  In  the  present  analysis,  the  theory 
of  Bernoulli-Euler  beams  is  used  with  the  assumption  of  small  deflections. 
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FIBERGLASS  COMPOSITE  BEAM 


Figure  (1)  -  A  schematic  drawing  of  the  cross  section  of 

NITINOL-reinforced  composite  beam. 


X 

A 


Figure  (2)  -  NITINOL-reinforced  beam  element  with  forces  and 

resulting  displacements. 


2.1.  EXTERNAL  WORK 


» 


» 


i 
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The  work  done  by  the  external  loads  includes; 

a.  work  done  by  transverse  loads  and  moments  (Wj ) 

This  work  is  given  by 

Wi  =  1/2  [5]"^  [F],  (1) 

where  [5]  and  [F]  are  the  displacement  and  transverse  loads  vectors, 
respectively,  given  by 


II 

dj  Wj 

(2) 

and 

[F]  =  [Vj 

Ml  Vj 

Mj]"^. 

(3) 

where  Wj  and 

are  the 

linear 

and  angular 

deflections 

of 

node 

respectively  and 

and  Ml 

are  the  shear  and 

moment  acting 

at 

node 

respectively. 

b.  work  done  by  the  axial  mechanical  loads  (Wg^) 

W2n,  is  given  by  [8]  as  ^ 

W2m  =  Pn.  /  2  I  (dw/dx)2  dx,  (4) 

0 

where  is  the  external  axial  compressive  load  acting  along  the  neutral 
axis  of  the  beam  element. 

c.  work  done  by  the  axial  thermal  loads  (Wj^,) 

^2t  represents  the  work  done  by  the  thermal  loads  on  the  beam 
element  due  to  changes  in  the  temperature  LB  of  the  element  caused  by 
changes  in  the  ambient  temperature  or  during  the  activation  and 
de-activation  of  the  NITINOL  fibers.  It  is  given  by 

L 

Wgt  =  Pt  /  2  J  (dw/dx)^  dx,  (5) 

0 

where  is  given  by 

P,  =  a  A0  E„  A„,  (6) 


i 
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where  a  is  the  thermal  expansion  coefficient  of  the  composite,  is  its 
modulus  of  elasticity  and  A^,  is  the  beam  cross  sectional  area. 

2.2.  STORED  STRAIN  ENERGY 

The  stored  strain  energy  consists  of  two  components: 

a.  strain  energy  of  beam  (W3) 

The  energy  stored  in  the  beam  element  due  to  its  bending  [9]  is  given 
by 

L 

W3  =  Em  Ira  /  2  J  (d^w/dx^)^  dx  (7) 

0 

where  is  the  flexural  rigidity  of  the  beam. 

b.  strain  energy  of  NITINOL  fibers  (W4) 

Considering  the  NITINOL  fiber  as  a  string,  as  shown  by  Baz  et  al  [10], 
with  a  tension  T  which  is  displaced  laterally  a  di  '•ance  w  from  the  neutral 
axis  of  the  beam,  then  its  stored  strain  energy  W4  is  given  by  [9] 

L 

W4  =  T  /  2  I  (dw/dx)2  dx.  (8) 

0 

Equating  the  sum  of  the  work  done  by  the  external  forces  F,  and 
to  the  sum  of  the  strain  energies  stored  in  beam  and  NITINOL  fibers  gives 

Wi  +  W2„  +  U2,  =  W3  +  W4.  (9) 

Substituting  equations  (1),  (4),  (5),  (7)  and  (8)  into  equation  (9)  yields 

L  L 

[5]^[F]  =  E^  I^  J  (d^/dx^)^  dx  -  Pn  J  (dw/dx)^  dx,  (10) 

0  0 


where  P„  is  the  net  axial  force  give  by 


I 


I 


► 


i 


Pn  =  {  Pm  ^  Pt  -  T).  (11) 

Defining  a  cubic  displacement  function  for  the  composite  beam  element, 
of  the  following  form 

w  =  a  +  b  X  +  c  +  d  x^.  (12) 

where  a,  b,  c  and  d  are  constants  that  can  be  calculated  in  terms  of  the 
deflections  of  the  nodes  i  and  j  bounding  the  beam  element,  equation  (12) 
can  be  rewritten  as 

w  =  [A]  [8],  (13) 

where  the  elements  of  matrix  [A]  are  function  of  x  (Fenner  [11]). 

Accordingly,  dw/dx  and  d^/dx^  can  be  obtained  by  differentiating 
equation  (13)  with  respect  to  x  to  yield 

dw/dx  =  [C]  [5]  and  d^/dx^=  [D]  [5],  (14) 

where  the  matrices  [C]  and  [D]  are  given  by 

[C]=^([A])  and  (DI  *  ^2  ([A]).  (15) 

The  following  relationships  can  also  be  obtained  from  equation  (14) 

(dw/dx)^  ==[6]^[C]^[C]  [5]  and  (d^/dx^)^  =[6]^[D]^[D]  [5] .  (16) 

If  the  stiffness  matrix  [K^]  of  the  element  is  defined  by  the  following 
relationship 

[FI  =  [K,]  [5],  (17) 

then,  [Kg]  can  be  determined  by  combining  equations  (10),  (16)  and  (17)  as 

follows 

L  L 

[Kg]  =  I„  I  [D]'^[D]  dx  -  P„  J  [C]^[C]  dx.  (18) 

0  0 

It  can  be  seen  from  equation  (18),  that  the  element  stiffness  matrix 
[Kg]  consists  of  two  components:  the  conventional  transverse  stiffness  and 
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the  geometric  stiffness  that  combines  the  effect  of  the  axial  mechanical 
loads,  axial  thermal  loads  and  the  tension  of  the  NITINOL  reinforcing 
fibers.  Equation  (18)  also  represents  the  basic  equation  for  understanding 
the  role  that  the  NITINOL  fibers  can  play  in  controlling  the  static 
characteristics  of  the  composite  beam.  For  example,  if  the  beam  is  not 
reinforced  by  NITINOL  fibers  (i.e.  T  =  0)  and  the  mechanical  and  thermal 
loads  induce  compressive  stresses  in  the  beam,  then  the  geometric  stiffness 
will  increase  and  the  total  element  stiffness  will  decrease.  When  the 
combined  effect  of  the  mechanical  and  thermal  loads  reaches  a  critical 
magnitude  such  that  the  geometric  stiffness  becomes  equal  to  the  flexural 
stiffness  of  the  beam,  the  beam  stiffness  vanishes  and  the  beam  becomes 
elastically  unstable.  Subjecting  the  beam  under  this  condition  to  any 
additional  external  disturbance  will  cause  the  beam  to  buckle. 

It  should  be  pointed  out  that  the  thermal  loading,  as  it  increases  the 
geometric  stiffness,  also  decreases  the  flexural  stiffness  of  the  beam 
because  it  reduces  its  effective  modulus  of  elasticity  E^.  Such  a  dual 
effect  makes  the  beam  buckle  under  smaller  thermal  loads  than  under  pure 
mechanical  loads. 

However,  the  critical  load  of  the  un-reinforced  beam  can  be  increased 
by  embedding  pre-strained  NITINOL  fibers  into  the  beam.  If  the  tension  T, 
resulting  from  the  pre-strain  alone,  is  high  enough  to  counter-balance  the 
mechanical  and  thermal  effects  then  the  beam  stiffness  can  be  maintained 
unchanged.  For  higher  pre-strain  levels,  the  beam  stiffness  can  be 
enhanced.  Further  enhancement  can  be  achieved  when  the  shape  memory  effect 
of  the  NITINOL  fibers  is  activated  by  heating  the  fibers  above  their 
austenite  phase  transformation  temperature.  The  additional  tension, 
induced  into  the  fibers  by  the  phase  recovery  forces  makes  the  net  axial 
load  Pn  negative,  and  accordingly  increases  the  overall  stiffness  of  the 
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beam  element.  However,  it  is  essential  that  the  total  tension  in  the 
NITINOL  fibers,  i.e. ,  the  sum  of  the  tension  due  to  the  pre-strain  and  the 
phase  recovery  force,  must  exceed  the  mechanical  and  thermal  loads  and 
compensate  for  the  softening  effect  produced  in  the  matrix  by  the  heating 
the  NITINOL  fibers.  Therefore,  effective  control  of  the  stiffness  of 
NITINOL-reinforced  composites  can  be  achieved  by  proper  selection  of  the 
initial  pre-strain  level  of  the  NITINOL  fibers.  This  selection  is 
particularly  crucial  as  the  pre-strain  level  determines  the  generated 
levels  of  recovery  forces  [10]. 

The  finite  eiement  model  of  the  NITINOL-reinforced  beams  describes  the 
interaction  between  the  external  loads,  operating  conditions  and  the 
geometrical  and  physical  parameters  of  the  composite  beam  and  the  NITINOL 
fibers.  It  defines  how  the  NITINOL  fibers  can  be  utilized  to  tailor  the 
stiffness  of  the  composite  to  compensate  for  environmental  and  operating 
conditions  and  disturbances.  The  stiffness  obtained  for  the  individual 
elements  of  the  beam  can  be  assembled  using  the  classical  finite  element 
approach  (Fenner  [11]).  The  assembled  model  can  then  be  subjected  to  the 
appropriate  boundary  conditions  in  order  to  compute  the  deflections 
corresponding  to  particular  external  loading  conditions.  The  analysis 
presented  is  for  an  orthotropic  laminate  that  has  a  single  layer  of 
unidirectional  NITINOL  fibers.  Such  an  analysis  can  be  used  along  with  the 
classical  laminate  theory  (Vinson  and  Sierakowski  [12])  to  assemble  the 
stiffness  matrix  for  a  multi-laminate  composite  beam.  A  similar  approach 
has  been  carried  out  for  modeling  the  static  and  dynamic  characteristics  of 
NITINOL-reinforced  composite  plates  (Baz,  Ro  and  Gilheany  [13]).  The 
finite  element  model  developed  will  be  validated  with  experimental  results 
obtained  with  fiberglass  composite  beams. 
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3.  DYNAMIC  CHARACTERISTICS  OF  THE  NITINOL-REINFORCED  BEAMS 


The  dynamic  characteristics  of  NITINOL-reinforced  beams  are  obtained  by 
combining  the  stiffness  matrix  [K^J  with  the  mass  matrix  [M,.]  of  the  beam 
to  form  the  following  element  equations  of  motion 

[M,]  [5]  +  [KJ  [5]  *  (FI.  (19) 

where  [5]  is  the  nodal  acceleration  vector.  The  elements  m^d.j)  of  the 
element  mass  matrix  [Mg]  are  obtained  using  the  consistent  mass  formulation 
(Zienkiewicz  and  Taylor  [14])  as  follows 

L 

=  Pm  J  Ai  Aj  dx.  (20) 

0 

where  Aj  and  Aj  are  the  ith  and  jth  elements  of  the  vector  A  given  by 
equation  (13). 

The  classical  finite  element  approach  is  used  to  form  the  equations  of 
motion  of  the  assembly  of  several  beam  elements  and  the  appropriate 
boundary  conditions  are  then  applied.  The  solution  for  the  eigenvalues  of 
the  resulting  homogeneous  equations  give  the  natural  frequencies  of  the 
composite  beam  as  influenced  by  the  properties  of  the  composite  matrix  and 
the  NITINOL  fibers.  It  is  important  to  note  that  these  properties  are 
influenced  by  the  temperature  distribution  inside  the  beam  which  is 
developed  by  virtue  of  activating  and  de-actlvating  the  NITINOL  fibers. 

4.  EXPERIMENTAL  VALIDATION 

4. 1  Test  facility  and  experimental  beam 

The  characteristics  of  NITINOL-reinforced  beams  are  computed  using  the 
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developed  static  and  dynamic  models.  The  theoretical  predictions  are 
compared  with  experimental  results  obtained  with  a  composite  beam  made  of 
randomly  oriented  glass  fibers  embedded  in  a  low  cure  temperature  polyester 
resin.  The  beam  is  30  cm  long,  2.5  cm  wide  and  0.156  cm  thick  and  is 
mounted  with  fixed-fixed  boundary  conditions.  The  temperature  dependence  of 
the  modulus  of  elasticity  of  the  beam,  shown  in  Figure  (3),  is  obtained 
experimentally  using  the  Dynamic,  Mechanical,  and  Thermal  Analyzer  (DMTA) 
of  Polymer  Laboratories,  Ltd  [15]. 

Four  NITINOL  55  fibers,  that  are  0.55  mm  in  diameter,  are  embedded 
inside  the  beam  through  vulcanized  rubber  sleeves  that  have  outer  diameter 
of  0.95  mm.  Two  sets  of  NITINOL  fibers  were  used.  The  first  set  consisted 
of  trained  fibers  that  have  an  austenite  transformation  temperature  of 
50°C.  However,  the  second  set  is  untrained  and  the  shape  memory  effect  has 
not  been  imparted  to  it.  The  two  sets  are  inserted,  one  at  a  time,  inside 
the  sleeves  and  the  effect  of  the  shape  memory  and  the  associated  phase 
recovery  forces  on  the  performance  of  the  composite  beam  are  monitored  when 
the  beam  is  exposed  to  different  ambient  temperatures.  The  experimental 
set-up,  shown  in  Figure  (4),  is  placed  inside  a  temperature-controlled 
chamber  to  determine  the  natural  frequencies  of  the  fixed-fixed  beam  as  a 
function  of  the  ambient  temperature.  In  the  set-up  the  NITINOL-reinforced 
beam  is  fixed  at  both  ends,  whereas  the  NITINOL  fibers  are  clamped  in  a 
holder  at  one  end  and  connected  to  a  load  cell  at  the  other  end.  The  load 
cell  monitors  the  pre-strain  level  of  the  fibers  when  they  are  in  their 
martensitic  phase,  as  well  as  continuously  measuring  the  recovery  force 
when  the  fibers  undergo  phase  transformation.  In  this  arrangement,  the 


fibers  are  activated  thermally  by  controlling  the  temperature  of  the 
environmental  chamber.  The  measurements  are  carried  out  after  steady-state 
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Figure  (3)  -  Effect  of  operating  temperature  on  modulus  of  elasticity 

of  the  test  fiberglass  composite  beam. 


Figure  (4)  -  Experimental  set-up  used  In  monitoring  the  performance 

of  NITlNOL-relnforced  beams. 


and  thermal  equilibrium  conditions  are  attained.  Under  these  conditions, 
the  composite  matrix  and  the  NITINOL  fibers  are  all  at  the  same  equilibrium 
temperature.  At  each  equilibrium  temperature,  the  composite  beam  is 
subjected  to  random  vibrations  and  the  resulting  response  is  monitored  by 
an  micro-accelerometer  bonded  to  the  beam.  The  response  is  analyzed  in  the 
frequency  domain  to  determine  the  modes  of  vibration  of  the  composite  beam. 

4.2  Natural  frequencies  of  NITINOL-reinforced  beams 

Figure  (5-a)  shows  the  measured  changes  in  the  first  natural  frequency 
of  the  beam  when  it  is  reinforced  with  untrained  NITINOL  fibers  which  are 
pre-strained  at  different  levels.  The  changes  are  normalized  with  respect 
to  the  natural  frequency  Wj,  of  the  un-reinf orced  beam  measured  at  25°C, 
i.e.  50.1  Hz.  The  normalized  characteristics  of  the  un-reinforced  beam  are 
also  plotted  to  serve  as  a  datum  for  defining  both  the  effect  of 
reinforcing  the  beam  with  NITINOL  fibers  and  the  effect  of  the  pre-strain 
level.  It  can  be  seen  that  the  frequency  of  the  un-reinforced  beam  drops 
as  the  ambient  temperature  increases  and  the  beam  losses  its  elastic 
stability  and  start  to  buckle  when  the  temperature  exceeds  40°C  .  The  drop 
in  the  natural  frequency  of  the  un-reinforced  beam  is  attributed  to  the 
softening  of  the  matrix  which  is  clearly  demonstrated  by  the  loss  in  the 
modulus  of  elasticity  of  the  beam  as  shown  in  Figure  (3). 

Reinforcing  the  beam  with  pre-strained  untrained  NITINOL  fibers 
considerably  increases  the  natural  frequency  of  the  beam.  The  extent  of 
the  upward  shift  in  natural  frequency  increases  with  increased  pre-strain 
level.  An  increase  of  about  40%  is  obtained  at  room  temperature  when  the 
pre-strain  level  is  only  0.26%.  However,  as  the  ambient  temperature 
increases  the  frequency  shift  drops  in  a  manner  similar  to  the 
characteristics  of  the  plain  un-reinforced  beam.  Such  a  drop  is  again 
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Effect  of  the  ambient  temperature  and  pre-strain  level  on 
enhancing  the  first  mode  of  vibration  of  fixed-fixed 
composite  beam  reinforced  with  NITINOL  fiber  without  (a) 
and  with  (b)  shape  memory. 


attributed  to  the  softening  effect  of  the  matrix  and  the  fact  that  the 
untrained  NITINOL  fibers  act  as  a  static  pre-tensioning  device  that  produce 
constant  tension  which  is  independent  of  temperature.  Therefore,  the 
frequency  enhancement  is  only  generated  by  the  reinforcement  and  the 
pre-strain  effects,  and  not  by  the  shape  memory  effect.  It  is  important 
to  note  that  considerably  higher  increases  in  the  natural  frequencies  can 
be  obtained  by  further  increases  of  the  pre-strain  level  up  to  its  maximum 
permissible  level  of  67.. 

However,  a  greater  frequency  shift  can  be  achieved  by  imparting  the 
shape  memory  effect  to  the  NITINOL  fibers.  The  fibers  are  trained  over  250 
cycles  using  the  procedure  outlined  by  Johnson  [16].  The  trained  fibers 
are  inserted  into  the  composite  beam  to  replace  the  untrained  set  and  the 
frequency  shifts  become  significant,  particularly  at  high  ambient 
temperatures.  This  is  clearly  demonstrated  in  Figure  (5-b).  For 
temperatures  between  room  temperature  and  40°C,  the  frequency  shifts 
obtained  are  similar  to  those  with  the  untrained  fibers  within  experimental 


accuracy. 

Once  the  ambient 

temperature 

exceeds 

the  50°C,  i.e. 

the 

austeni te 

phase  transformation 

temperature 

of  the 

NITINOL  fibers. 

the 

frequency  shift  characteristics  changes  from  a  gradually  decaying  trend  to 
one  that  is  a  gradually  increasing.  Such  a  sudden  change  is  the  result  of 
the  contribution  of  the  phase  recovery  forces  developed  by  the  shape  memory 
effect  which  is  illustrated  in  Figure  (6). 

The  shape  memory  effect  generates  strain  energy  in  the  NITINOL  fibers 
to  counterbalance  the  softening  effect  of  the  composite  matrix  with 
increasing  temperature.  As  the  amount  of  strain  energy  developed  depends 
on  the  initial  pre-strain  level,  it  can  merely  compensate  for  the  softening 
effect  to  maintain  the  beam  frequency  at  nearly  a  constant  value  which  is 
independent  of  ambient  temperature  as  shown  for  pre-strain  level  of  0.0787.. 
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It  can  also  increase  the  beam  frequency  as  the  ambient  temperature 
increases  as  indicated  for  pre-strain  levels  of  0.22  and  0.26%.  For  a 
pre-strain  level  of  0.26%  and  ambient  temperature  of  90°C,  the  frequency 
increase  reaches  about  70%  as  compared  to  an  18%  increase  when  untrained 
fibers  are  used.  In  this  manner,  the  individual  contributions  of  the 
pre-strain,  matrix  softening  and  shape  memory  effect  on  the  frequency  shift 
are  isolated.  This  facilitates  checking  the  validity  of  the  mathematical 
models  against  the  experimental  results. 

Comparisons  between  the  theoretical  predictions  and  the  measurements 
are  shown  in  Figures  (7-a)  and  (7-b)  for  NITINOL  fibers  without  and  with 
shape  memory  effect,  respectively.  The  figures  include  comparisons  for  the 
first  and  second  modes  of  vibrations.  Close  agreement  between  theory  and 
experiments  is  evident. 


5.  CONCLUSIONS 

The  static  and  dynamic  characteristics  of  NITINOL-reinforced 
composite  beams  have  been  presented.  The  fundamental  issues  governing  the 
behavior  of  this  new  class  of  SMART  composites  have  been  introduced. 
Applications  of  NITINOL  reinforcing  to  control  the  static  and  dynamic 
behavior  of  composite  beams  have  been  successfully  demonstrated. 

Emphasis  has  been  placed  in  the  presentation  on  the  actuation 
capabilities  of  the  NITINOL  fibers.  However,  extensive  efforts  are  in 
progress  to  use  the  NITINOL  fibers  to  extract  modal  and  physical 
displacements  of  structures  with  multi-modes  of  vibration  (Baz,  Poh  and 
Gi Iheany  [  17 ]  ) . 

With  such  built-in  sensing  and  controlling  capabilities, 
NITINOL-reinforced  composites  can  provide  a  means  for  continuously  tuning 
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Figure  (7) 


Comparison  between  the  theoretical  and  experimental 
frequencies  of  composite  beam  reinforced  with  NITINOL 
fibers  without  (a)  and  with  shape  memory  (b). 


the  structural  characteristics  to  adapt  to  changes  in  the  operating 
conditions.  These  features  will  be  particularly  useful  in  many  critical 
structures  that  are  intended  to  operate  autonomously  for  long  durations  in 
isolated  environment  such  as  defense  vehicles,  space  structures  and 
satellites. 
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Wi 

Wsm 

Wat 
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interpolating  function  of  beam  deflection 

ith  element  of  [A] 

cross  sectional  area  of  beam 

[D]  first  and  second  derivatives  of  interpolating  function 
of  beam  deflection 
Young’ s  modulus  of  beam 
vector  of  external  loads  acting  on  beam 
area  moment  of  inertia  of  beam 
stiffness  matrix  of  beam  element 
length  of  beam  element  and  NITINOL  fiber 
external  moment  acting  at  ith  node 
mass  matrix  of  bean  element 
jj  the  element  i.J  of  the  mass  matrix 

mechanical,  net  and  thermal  axial  loads  acting  on  beam 
generalized  coordinate  of  the  nth  vibration  mode  of 
NITINOL  fiber 

generalized  acceleration  of  the  nth  vibration  mode  of 
NITINOL  fiber 
time 

initial  tension  in  a  NITINOL  fiber 

total  tension  in  a  NITINOL  fiber 

shear  force  acting  at  the  ith  node 

transverse  deflection  of  beam  and  NITINOL  fibers 

work  done  by  transverse  loads 

work  done  by  mechanical  axial  loads 

work  done  by  thermal  axial  loads 

strain  energy  of  beam 

strain  energy  of  NITINOL  fiber 

z  cartesian  coordinates  along  beam  neutral  axis  and  cross 
section  respectively 


Greek  letters 

[5]  the  deflection  vector  of  beam  element 

A0  axial  temperature  difference 

angular  deflection  of  ith  node 
mode  shape  of  the  nth  mode 
Wj,  natural  frequency  of  the  nth  mode 
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abstract 

An  Adaptive  Modal  Positive  Position  Feedback 
(AMPPF)  method  is  presented  for  controlling  the 
vibration  and  shape  of  flexible  structures.  The 
proposed  strategy  combines  the  attractive  attributes  of 
the  Independent  Modal  Space  Control  (IMSC)  of 
Heirovltch  and  Positive  Position  Feedback  (PPF)  of  Goh 
and  Caughy.  The  control  ier  is  designed  in  the 
uncoupled  modal  space  using  only  modal  position  signals 
to  damp  the  vibration  of  undamped  modes.  The 
parameters  of  the  AMPPF  controller  are  also  adjusted  in 
an  adaptive  manner  in  order  to  follow  the  performance 
of  an  optimal  reference  model.  In  this  way.  optimal 
damping  and  zero  steady-state  errors  can  be  achieved 
even  in  the  presence  of  uncertain  or  changing 
structural  parameters. 

The  adaptation  laws  governing  the  stable  variation 
of  the  AMPPF  controller  parameters  are  derived  using 
Lyapunov  stability  theorem.  The  effectiveness  of  the 
AMPPF  in  controlling  the  vibration  and  shape  of  a 
variable  mass  cantilevered  beam  is  demonstrated 
cxper imenta 1 ly .  The  performance  obtained  with  the  AMPPF 
algorithm  is  compared  with  those  of  other  ciassicai 
control  algorithms.  The  results  obtained  emphasize  the 
potential  of  the  AMPPF  algorithm  as  an  efficient  means 
for  controlling  flexible  structures  with  uncertainties 
in  real  time. 

1.  INTRODUCTION 

Considerable  emphasis  has  been  recently  placed  on 
the  development  of  wide  variety  of  adaptive  control 
algorithms  to  effectively  control  the  vibration  and 
shape  of  large  flexible  structures.  Such  algorithms  aim 
at  compensating  for  the  problems  arising  from  the 
uncertainty  of  accurateiy  modeling  the  dynamics  of 
these  structures,  truncation  of  their  dynamic  models 
and  the  associated  control  spillover  from  the  residual 
modes,  variation  of  the  structural  parameters  and 
nonstat ionar i ty  of  the  disturbances  acting  on  the 
structures.  Without  the  appropriate  adaptation  to 
these  problems,  fixed  gain  controllers  become  totally 
inadequate  to  meet  the  control  objectives  and  the 
desired  performance  requirements. 

Distinct  among  the  recently  developed  adaptive 


control  algorithms  is  that  of  Bar-Kana  and  Kaufman  (1) 
which  relied  in  its  operation  on  collocated  pairs  of 
sensors/actuators  and  on  feeding  back  position  and 
velocity  outputs.  The  controller  employed  the  direct 
reference  model  approach  to  theoretically  control  the 
vibration  of  beams.  plates  and  frames.  Similar 
approach.  but  with  proportional  and  integral 
controller,  has  been  utilized  by  Ih  et  al  (Zl  to 
experimentally  control  a  5.  6  m  antenna  at  the  JFL. 
K.Ossman  et.al,  [31  have  theoretically  developed  an 
adaptive  model  following  controller,  for  SCOLS.  in 
which  position  and  velocity  outputs  are  forced  via  a 
variable  structure  control  to  track  reference  position 
and  velocity  paths.  K.Ossman  et.al.  developed  also  an 
indirect  adaptive  controller  in  which  the  system 
parameters  are  estimated  by  a  recursive  least-squares 
and  an  LQ  controller  is  designed  to  quickly  damp  out 
the  vibration  (3).  The  effect  of  the  spillover  from 
the  unmodeled  modes  is  found  to  cause  the  control 
Inputs  of  this  algorithm  and  the  system  outputs  to  grow 
without  bound.  Siiverberg  and  Norris  [41  devised  a 
self-tuned  indirect  adaptive  controller  which 
identifies  the  structural  dynamics  and  updates  the 
associated  control  gains  to  uniformly  damp  out  the 
vibration  of  the  structure. 

In  all  the  above  algorithms,  the  indirect  adaptive 
controllers  and  the  direct  reference  modeis  are 
designed  In  the  physical-coupled  space.  For  large 
structures  this  presents  serious  computational 
challenges  particularly  when  Kalman  filtering  Is  needed 
as  in  Ref.  [2).  In  the  present  study,  the  proposed 
Adaptive  Modal  Positive  Position  Feedback  (AMPPF) 
controller  and  its  reference  model  are  designed 
completely  in  the  independent  modal  space  with  the 
open-loop  equations  of  the  system  remaining  uncoupled 
even  after  including  the  modal  controller  Also,  the 
AMPPF  uses  only  modal  position  signals  to  obtain  stable 
and  damped  performance.  Such  performance  is  attained 
by  positively  feeding  back  the  position  signals  through 
tuned  first  order  filters.  The  performance  of  the 
AMPPF  algorithm  is  enhanced  by  augmenting  it  with  a 
"time  sharing"  strategy  that  utilizes  small  number  of 
actuators  to  control  larger  number  of  modes.  The  AMPPF 
algorithm  is  based  on  the  fixed  parameter  Modal 
Positive  Position  Feedback  (MPPF)  algorithm  developed 
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by  Bai  et.al.  which  has  been  successfully  utiliied  to 
control  the  vibration  of  sirnple  beams  151  and  more 
complex  structures  161. 

The  present  study  aims  at  developing  the 
adaptation  laws  of  the  AMPPF  algorithm  and 
experimentally  evaluate  its  effectiveness  in 
controlling  the  vibration  and  shape  of  a  cantilevered 
beam. 

This  paper  is  organized  in  five  sections.  A  brief 
introduction  is  given  in  section  1.  The  concept  of 
the  MPPF  and  the  AMPPF  algorithms  are  presented  In 
sections  2  and  3  respectively.  Section  4  includes  the 
experimental  evaluation  of  the  algorithm  along  with 
comparisons  with  other  algorithms.  Section  S 
summarizes  the  conclusions  of  this  study. 

2.  MODAL  POSITIVE  POSITION  FEEDBACK  (MPPF) 

2.1.  Concept 

The  basic  concept  of  the  MPPF  method  can  be 
clearly  understood  by  considering  the  diagram  shown  in 
Figure  (1).  where  the  controller  is  used  to  control 
the  i^  mode  of  an  undamped  flexible  structure  in  the 
independent  modal  space.  The  controller  feeds  back 
positively  the  modal  displacement  u,  through  a  first 
order  filter  that  has  a  time  constant  t,.  The  filter 
output  Y|  is  added  to  the  desired  reference  modal 
displacement  Up,  and  resulting  signal  is  amplified  by  a 
proportional  controller  gain  K,.  This  gain  is  set  equal 
to  to  be  in  a  form  similar  to  that  of  Goh  and 
Caughey' s  (7],  where  u,  is  the  natural  frequency  of  the 
i  mode.  The  amplified  signal  f,.  i.e.  the  modal 
control  action,  is  then  sent  to  control  the  1  mode  of 
the  structure. 


Figure  (1)  -  Block  diagram  of  the  MPPF  controller. 

Mathematically.  the  interaction  between  the 
structural  mode  and  the  controller  can  be  described  as 
f  ol  lows: 

Structure:  Uj  ♦  u,  =  f,  =  Tj 

Filter:  T,  Y,  ♦  Y,  =  u,  .  (2) 

The  above  structure-filter  system  has  the 

following  closed-loop  transfer  function 

u,/Up,=3r,uf  (t,s*1  )/[T,s^+s^*T,ufs*ufll-T| )  1 .  (3) 

where  s  is  the  Laplace  operator.  Applying  Routh's 
stability  criterion,  the  system  is  asympotical ly  stable 
for  value  of  0  s  y,  s  1  and  Tj  z  0.  Accordingly.  It 
is  possible  for  an  undamped  system  to  attain  asymptotic 
stability  by  feeding  positively  its  position  signal 
through  a  simple  first  order  filter  without  the  need 
for  any  velocity  feedback.  This  constitutes  the  basic 
premise  of  the  present  control  algorithm. 
Conceptually,  the  algorithm  possesses  this  favorable 
stable  performance  because  the  first  order  filter  Is 
in  effect  equivalent  to  an  Integral  controller  with 
negative  feedback.  It  is  also  important  to  note  that  if 


the  modai  position  signal  is  fed  back  negatively, 
instead  of  positively,  through  the  first  order  filter, 
the  system  will  always  be  unstable 

2. 2.  Parameters  of  MPPF 

2.2.1  controller  gain  (j,) 

The  implementation  of  this  modal  control  algorithm 
requires  the  selection  of  two  design  parameters,  i.e  . 
y,  and  t,.  Actually,  only  the  time  constant  t,  of  the 
filter  needs  t^  be  selected  since  y,  must  assume  a 
fixed  value  y,*  0.5  to  eliminate  the  controller 

steady-state  error  as  indicated  by  equations  (1)  and 
(2).  Such  a  value  is  <  1  and  satisfies  the  asymptotic 
stability  condition  previously  discussed.  Therefore, 
this  algorithm  can  be  equally  used  for  accurate  shape 
control  (u^,  *  0  ). 

2.2.2.  filter  time  constant  (Tj) 

The  optimal  value  of  the  time  constant  x,  of  the 
filter  is  determined  by  dividing  the  numerator  and 
denominator  of  equation  (3)  by  T,Wi  to  yield  the 
following  equation 

u,/UR|=Tia^  (s  /a,*l  )/[s^*ajS^*s-(  1-y,  )  a,  1 .  U) 

where  a,  =  1  t,u,  and  s  =  s  /  u.  The  above  system 
has  the  following  characteristic  equation 

s’  ♦  “i  s’  X  s  X  (1-yi)  a,  =  0.  (51 

wjjlch  has  the  root  locus  plot  shown  in  Figure  (2)  for 
y,  =0.5  and  0<O|<co.  The  corresponding  damping  ratio 
of  the  closed-loop  system  obtained  from  the  root  locus 
plot,  is  shown  in  Figure  (3)  as  a  function  of  a,  which 
is  the  only  design  parameter  of  the  system.  Figure  (31 
indicates  that  the  damping  ratio  attains  a  maximum 
value  of  20.07  X  when  a,  =  1  18  This  optimal  value  is 
very  close  to  the  PPF  results  obtai.-.ed  experimentally 
by  Fanson  and  Caughey  181. 


b.i  , 


Figure  (2)  -  Root  locus  of  the  MPPF  method  with  y|=0.S 


Tnnf 

Figure  (3)  -  Effect  of  a,  on  closed-loop  damping  of 
MPPF  method  with  y,=0.5 


In  summary,  using  first  order  filters,  instead  of 
the  second  order  filters  of  Goh  and  Caughey.  has 
simplified  the  design  without  compromising  the  damping 
characteristics  of  the  controller.  More  importantly, 
since  the  analysis  presented  Is  applicable  to  any  mode, 
uniform  damping  for  all  the  modes  results  if  all  their 
filters  are  tuned  to  satisfy  the  optimal  tuning 
condition  (  a,  =  1.18  ).  Accordingly,  a  damping  ratio 
of  20.07  ■/.  can  be  maintained  for  any  mode  1.  that  has 
natural  frequency  u,.  by  selecting  the  time  constant  x, 
of  its  filter  such  that 

T|  =  1/118U|  forl  =  1.2....N.  (6) 

where  N  is  the  number  of  controlled  modes. 

3.  ADAPTIVE-MPPF  with  PAPAMETEH  ADAPTATION  (AMPPF) 

3.1.  Need  for  parameter  adaptation 

In  order  to  effectively  utilize  the  MPPF  method  in 
controlling  the  shape  and  vibration  of  flexible 
structures  in  an  optimal  manner,  it  is  essential  to 
investigate  other  inherent  features  of  the  method. 
Such  features  can  be  easily  revealed  by  considering  the 
root  loci!  of  the  MPPF  characteristic  equation  (5). 
shown  in  Figure  (4).  for  different  values  of  the  gain 
y,  as  a,  is  varied  from  0  to  a.  It  is  evident  that 
increasing  y,  beyond  0.5  while  maintaining  it  below  the 
stability  limit  (i.e.  y|Sl)  results  in  significant 
changes  in  the  shape  of  the  root  locii.  More 
Importantly,  the  branches  of  the  root  locii  converge 
towards  the  real  axis  of  the  s  plane  indicating  that 
higher  closed-loop  damping  ratios  can  be  achieved. 
Figure  (5)  shows  the  combined  effect  of  y,  and  a,  on 
the  attainable  damping  ratios  as  extracted  from  the 
plots  of  Figure  (4).  It  is  clear  that  for  any  value  of 
the  controller  gain  y,  there  Is  an  optimal  value  of  a, 
for  which  the  closed-loop  damping  ratio  assumes  its 
■"aximum.  In  section  2,  it  was  shown  that  If  y^  is  set 

equal  to  0.5  to  ensure  zero  steady-state  error;  a 

maximum  damping  ratio  of  0.207  is  attained  when  a,  is 
1.18.  Increasing  y,  to  0.9  requires  that  a,  be  1.768 
to  reach  a  maximum  damping  ratio  of  l.OO.  Such  a 

significant  increase  in  the  damping  ratio  occurs 
however  at  the  expense  of  increasing  the  steady-state 
error,  (1-U|/Uq,)  =  l-l/( 1/y , -1 ) .  to  800%,  Therefore, 
increasing  the  damping  ratio  above  0.207  cannot  be 
achieved  without  compromising  the  steady-state  error. 
This  will  only  occur  if  the  two  parameters  yj  and  a, 

are  maintained  constant  throughout  the  control  process 
at  values  other  than  0.5  and  1.18  respectively. 
However,  it  is  possible  to  combine  the  high  damping  and 
zero  steady-state  error  characteristics  by  adjusting 
the  parameters  and  a,  in  an  adaptive  manner. 

Specifically,  one  would  start  the  control  process  with 
the  highest  damping  possible  (i.e.  the  largest  value  of 
y,  )  and  as  the  process  approaches  its  completion  the 
gain  y.  Is  adjusted  to  be  as  close  as  possible  to  0.5 
in  order  for  the  structure  to  reach  Its  final  shape  or 
state  with  zero  steady  state  error.  During  the 

process  of  adjusting  the  gain  the  parameter  a,  is 
also  adjusted  accordingly  in  order  to  achieve  the 
maximum  damping. 

Accordingly,  the  AMPPF  method  capitalizes  on  the 
inherent  Interaction  between  the  controller  parameters 
to  optimize  the  transient  performance  of  the  system. 
In  this  manner,  although  the  structural  parameters  are 
time  invariant,  better  transients  can  be  achieved  by 
adaptively  making  the  controller  parameters  (yj  and  a,) 
time  varying  The  desired  time  response  can  be  obtained 
by  selecting  a  desired  reference  behavior  and  change 
the  controller  to  make  the  system  follow  the  reference 
mode  1 . 


Figure  (4)  -  Effect  of  y,  on  shape  of  root  locus 
of  the  MPPF  method. 
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Figure  (5)  -  Effect  of  a,  on  closed-loop  damping  of 
MPPF  method  for  different  gains  y,. 


The  interaction  between  the  reference  model  and 
the  structure-filter  systems  is  shown  in  the  block 
diagram  of  Figure  (6).  The  laws  governing  the 
adaptation  of  the  controller  parameters  yi  and  a,  (or 
T| )  are  derived  in  the  following  section. 

3.2.  Adaptation  laws 

Mathematically,  the  dynamics  of  the  structure  and 
its  model  reference  are  given  by  the  following 
equations: 


A.  Slnjclural  system 

Structure  dynamics;  Up*  Up=y,Wp  (Yp-^Up,). 
Filter  dynamics:  x,  Yp  ♦  Yp  =  Up  . 


I 


i;.  I.  I  .-I 


.■  . . I 


e  *  A,  e  *  [0 


(17) 


I 


» 


I 


I 
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To  calculate  the  adaptation 
Lyapunov  function  is  defined; 

V  =  e’^Pe  ♦  h  (0  .  0) 


law. 


>  0, 


the  following 
(18) 


where 


n 


lal 


^  0^  0,  ■  (19) 

1  =  1 


For  stable  adaptation.  V  has  to  be  strictly  negative 
(91.  l.e. 


V  =  e^Pe*e’^Pe  +  h<0.  (20) 

Substituting  e  from  equation  (15)  gives 

V  =  e^  [  P  *  P  A,)  e  +  2  e’'  P  f  ♦  li  (21  ) 


Figure  (6)  -  Blocl<  diagram  of  the  AMPPF  method  with 
parameter  adaptation. 


To  have  stable  adaptation  (  V  <  0  ).  the  following 
two  conditions  must  be  satisfied; 


B.  Reference  system 

Reference  model;  ('^•*^bi(' 

Reference  filter;  Y  ♦  ^2  ~ 


(3)  2  eV  f  ♦  h  =  0  (22) 

(b)  select  a  stable  reference  model  satisfying 
the  following  Lyapunov  equation; 


The  reference  model,  described  by  equation  (9).  is 
structured  to  be  damped  in  contrast  to  the  undamped 
structural  model  and  its  gain  is  selected  to  ensure 
stability  and  zero  steady-state  error  (i.e.  =  1/2). 

The  interaction  between  the  structure  dynamics  and 
the  reference  model  dynamics  is  displayed  in  the  block 
diagram  of  Figure  (6).  Defining  the  state  vector  Xp 
and  X„  as  Xp=  (u,  u^  Y^l^  and  X,=  (u,  u«  Y„1  .  the 
structure  and  reference  systems  can  be  written  as 


Structure  system  Xp  =  Ap  Xp 
Reference  system  X.,  =  A.,  X., 


Bp  bfii- 
Ba  Up,. 


(11) 

(12) 


0  1  0 

2  n  2 

-u,  0  y,u. 

Bp  = 

0 

2 

7,u, 

1/y,  0  -I/y, 

0 

0  1  0 

a,  32  -  5a, 

and  B.,  = 

0 

-.5a, 

a,  0  -aj 

0 

A^P  +  PA.,  =  -  Q.  (23) 

where  (3  is  a  positive  definite  matrix,  then  equation 
(21)  reduces  to 

Vs—  e^Q  e  .  (24) 

Equation  (22)  gives  the  necessary  adaptation  laws 
for  varying  y,  and  t,  such  that  the  stability  is 
ensured.  i.e.  Lyapunov  stability  theorom  (9J  is 
satisfied.  This  adaptation  law  can  be  determined  by 
considering  the  characteristic  equation  of  the 
reference  system  which  can  be  written  as 

s^-(a2  -  a3)s^*l-a2  a3  -  a,ls*[-0.5ai  a3l  =  0.  (25) 

If  it  is  desired  that  the  system  should  behave 
like  an  optimal  ITAE  (  Integral  of  Time  multiplied  by 
Absoulte  Error  )  third  order  system,  then  it  should 
have  following  characteristic  equation  (10); 

s’  ♦  1.75  wp  s’  *  2.  15  Up  s  +  Up  =  0  .  (26) 

Selecting  Up=  u,  of  the  system,  and  matching  the 
coefficients  of  two  characteristic  equations  (25)  and 
(26).  a,.  32  and  a3  can  be  found. 

The  parameter  error  matrices  0  and  0  can  be  formed 
as  follows; 


The  matrices  A.,  and  B,  of  the  reference  model  are  given 
in  terms  of  the  three  coefficients  a,.  82  and  83. 

Defining  the  error  vector  e  and  its  derivative  e 

e  =  X„  -  Xp  and  e=  )<.,-Xp.  (14) 

and  subtracting  equation  (11)  from  (12)  yields 

e  =  *  (A,-  Ap)Xp  +  (B,-  Bp)Up,  .  A,e  ♦  f.  (15) 

where.  f  =  (  A^-  Ap)  Xp  ♦  (  B,-  Bp)  Up,.  (16) 

Letting  A„-  Ap  =  0  and  B.-  Bp  -  0  .  then 


0 

a,e  u’ 
(33-1/7, ) 


0  0 
32  -(0.  5a2*TiU, ) 

0  -(33-1/7, ) 


(27) 


<!>  =  B,  -  Bp 


0 

-(0. 5a,*y,uf) 
0 


(28) 


Forming  h  using  equation  (19)  yields 


3  3 

h  =  ^  J  ^ 

1=1  J=1 

=  (a,-u^i^*2(a3-l/T,  )^*a2*(0.  5a[+3r,uf ).  (29) 

Differentiating  equation  (29)  with  respect  to  time 
yields: 

h=  4  (aj  -  1/T,)  T|/  rf*  4  (0.  Saj  +  riuf)  (30) 

and  defining  the  error  vector  e 


1 

■  X,™ 

^ip 

’  ei  ' 

e  =  X„  -  Xp  = 

^2p 

= 

ez 

.  , 

.  '^P  . 

. 

then  equation  (22)  becomes 

( e , Pa^ejPi+ejPj ) [( aj  ♦uf ) X, p^aj/jp- (  Saj *Tuf ) ( Yp+Up, )] 

■►(eiPj^ejPs^ejPp)  [(aj-l/T,  )  (X,p-Yp)] 

+2(a3-l/Tj )t,/Tj*2( . Saj+rjuf)  =  0.  (32) 


are  used  to  monitor  the  physical  displacements  of  the 
three  nodes  in  the  transverse  direction.  The  position 
signals  are  sampled  by  a  386-based  micro-processor 
provided  with  an  input/output  board  which  has  a 
conversion  time  of  15  ps  and  a  resolution  of  12  bits 
The  board  analog  outputs  have  a  settling  time  of  30  ps. 

The  micro-processor  uses  the  three  sampled  signals 
to  compute  the  beam  angular  deflections  and  the  linear 
and  angular  velocities  of  the  nodes  The  computed 
state  variables  are  used  to  calculate  the  modal 
coordinates  of  the  flexible  system,  the  mode  that  has 
the  highest  modal  energy,  the  corresponding  optimal 
modal  control  force  f,,  the  physical  control  force 
and  the  necessary  voltage  v  to  be  sent  to  the 
piezo-actuator.  The  Implementation  of  these 
calculations  ,  i.e.  the  AMPPF  algorithm  ,  is  carried 
out  in  real  time  in  3.04  ms. 

Figure  (8)  outlines  a  flow  chart  of  the  AhPPF 
algorithm  indicating  its  main  steps. 
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The  above  equation  yields  the  following  two 
adaptation  laws  : 


Figure  (7)  -  Schematic  drawing  of  the  experimental  beam, 
piezo-actuator  and  sensors. 


T,=-0  5Tf(eiP3*e2P5*e3P5) (  X,p  -  Yp  ),  (33) 

and 

r  =-0.  5(eipa+e2P4*e3Ps)  x 

)X,p*a2X2p-(.  5a,niuf)(Yp*UR,) 

[ - j ^ - ]  (34) 

(j,  (0.  5aj  ) 


which  depend  on  the  elements  of  the  P  matrix.  The 
elements  of  the  P  matrix  depend  In  turn  on  the  Q  matrix 
of  the  Lyapunov  equation  (23)  and  on  the  parameters  a,. 
a2  and  a3  of  the  reference  model  as  Indicated  in  the 
appendix. 


4.  EXPERIMENTAL  PERFORMANCE  OF  OPTIMAL  MPPF  ALGORITHM 
4. 1  EXPERIMENTAL  SET-UP 


A  thin  rectangular  cantilevered  beam  is 
constructed  to  validate  the  developed  algorithm  .  The 
design  parameters  of  the  beam  are  given  in  Table  1. 

The  beam  is  controlled  by  one  piezo-electric  blmorph 
made  from  G1195-  ceramic.  The  actuator  is  available 
commercially  (model  number  R20S  )  from  Piezo-Electric 
Products.  Inc.,  Metuchen,  NJ  08840-4015  .  Table  2  lists 
the  main  design  parameters  of  the  actuator. 

The  experimental  beam  and  the  piezo-actuator  are 
arranged  as  shown  In  Figure  (7).  The  beam  is  divided 
into  three  active  elements.  Bonded  to  the  first 
element,  near  the  fixed  end  of  the  beam,  is  the 
piezo-actuator.  Three  non-contacting  position  sensors 
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Figure  (8)  -  Flow  chart  of  the  AMPPF  method. 
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4.2  Modal  characteristics  of  the  beam-actuator  system 

The  modal  characteristics  of  the  experimental 
beam  are  determined  theoretically  [51  and  validated 
experimentally  using  classical  modal  analysis 
technique  A  comparison  between  the  theoretical  and 
the  experimental  values  of  the  first  five  vibration 
modes  of  the  beam-actuator  system  is  given  in  Table  3. 
The  table  gives  also  the  modal  damping  as  calculated 
from  the  experimental  results  using  the  half  power 
approach  (111. 

4.3  Experimental  results 

In  all  the  experiments  conducted  in  this  study, 
the  beam  is  excited  at  its  second  mode  of  vibration  by 
applying  sinusoidal  excitation  of  20  volts  in  magnitude 
to  the  piezo-actuator.  The  excitations  are  maintained 
for  a  period  of  0. 15  seconds.  The  beam  is  either  left 
to  vibrate  freely  (i  e  Uncontrolled)  or  under  the 
action  of  one  modal  control  algorithm  or  another.  The 
above  excitation  form  is  selected  in  order  to  excite 
the  second  mode  of  vibrations.  Excitation  of  modes 
higher  than  the  second  would  require  faster 
micro-processor  in  order  to  sample  at  least  10  sample 
per  period  to  achieve  meaningful  control.  The 
uncontrolled  performance  is  used  as  datum  for  judging 
the  effectiveness  of  the  different  control  algorithms. 

4.3.1  Vibration  control 

The  time  response  of  the  uncontrolled  beam  is 
shown  in  Figure  {9-a)  indicating  a  very  low  natural 
damping  characteristics.  Figure  (9-b)  shows  the  beam 
time  response  when  it  is  controlled  by  the  AMPPF 
algorithm  The  values  of  the  controller  parameters 
•jTj.  T.,  Tj  and  Tj,  for  the  first  and  second  modes  of 
vibrations,  are  initially  set  at  0.99,  0.836,  0.99  and 
0  0138  respectively.  These  Initial  values  of  the  Ti’s 
and  the  a,'s  are  selected  to  ensure  maximum  initial 
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Figure  (9)  -  Time  response  of  beam 

a.  uncontrolled  and 

b.  controlled  with  AMPPF 


damping  During  the  adaptation  process,  the  resulting 
time  variations  of  the  gains  and  the  filter  time 
constants  ^  are  shown  in  Figures  (10-a)  a.nd  (10-bl 
respectively.  The  figures  indicate  that  all  the 
adaptive  gains  and  time  constants  converge  to  stable 
values  to  ensure  zero  steady-state  errors.  Also,  the 
corresponding  time  histories  of  the  modal  displacements 
of  the  beam  are  shown  in  Figures  (11-a)  and  (11-b)  for 
the  first  and  second  modes  respectively  The  figures 
demonstrate  the  effectiveness  of  the  AMPPF  algorithm  in 
following  the  dynamics  of  the  desired  reference  modal 
model. 

A  better  insight  into  the  effectiveness  of  the 
AMPPF  algorithm  can  be  gained  by  considering  the  Fast 
Fourier  Transform  (FFT)  of  the  beam  response.  Figure 
(12)  shows  the  frequency  content  of  the  response  of  the 
uncontrolled  beam  in  comparison  with  the  controlled 
beam.  These  characteristics  are  obtained  by  sampling 
the  beam  tip  position  signal  by  a  spectrum  analyzer  and 
performing  on  it  an  FFT  analysis.  The  figure 
effectiveness  of  the  AMPPF  method  in  simultaneously 
suppressing  the  vibration  of  the  first  two  modes  of 
vibration  using  a  single  piezo-electric  actuator 

Comparisons  between  the  theoretical  and 
experimental  time  responses  of  the  uncontrolled  and 
controlled  beam  are  shown  in  Figures  (13-a)  and  (I3-b) 
respectively.  The  displayed  results  show  close 
agreements  between  theory  and  experiments. 
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Figure  (10)  -  Time  history  of  the  control  parameters 
a.  gain  yj  2  H’’’®  constant  ^ 

4.3.2  Vibration  control  of  a  variable  mass  system 

A  better  Insight  of  the  effectiveness  of  the 
adaptive  control  algorithm  can  be  gained  by  adaptively 
controlling  a  variable  mass  system.  The  mass  variation 
of  the  flexible  system  Is  achieved  by  adding  a  mass  to 
the  free  end  of  the  beam.  This  additional  mass  changes 
the  first  mode  of  vibration  by  20 

Figure  (14)  shows  the  time  response  of  the 
controlled  beam  when  an  additional  mass  weighing  55  gm 
attached  to  the  free  end  of  the  beam.  Figures  (15-a) 
and  (15-b)  show  the  corresponding  time  histories  of  the 
adaptive  gains  Tj.y  lod  lime  constant  2 
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variable  mass  system. 

It  Is  evident  that  the  AMPPF  algorithm  has 
successfuly  adapted  to  the  mass  changes  and  has 
effectively  controlled  the  beam  vibrations. 
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Figure  (13)  -  Comparison  between  theoretical  and 

experimental  time  response  of  the  beam 
a.  uncontrolled  and  b.  controlled 
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Figure  (12)  -  Frequency  response  of  the  free  end  of 
the  controlled  and  uncontrolled  beam 

4.3.3  Shape  control 

The  use  of  the  AMPPF  algorithm  In  controlling  the 
shape  of  the  beam  system  Is  demonstrated  In  Figure  (16) 
when  the  beam  is  deflected  0.1  mm  off  its  initial 
zero-load  position.  The  resulting  time  response  of  the 
beam  and  the  corresponding  control  voltage  are  shown  in 
Figures  (16-a)  and  (16-b)  respectively.  The  two 
figures  clearly  demonstrate  that  the  applied  control 
voltage  reaches  a  steady-state  value  after  the  required 
shapes  are  attained.  The  model  following  capability  of 
the  AMPPF  algorithm  Is  Indicated  In  Figure  (16-c)  where 
the  modal  displacement  at  the  first  mode  of  vibration 
Is  shown  to  track  the  desired  reference  model  output 
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Figure  (14)  -  Time  response  of  the  variable  mass  beam. 


4.3.4  System  comparison 

Figure  (17)  shows  the  time  response  of  the  beam 
when  the  beam  is  controlled  by  the  non-adaptive  MPPF 
method.  A  comparison  with  the  beam  response  with  the 
AMPPF.  shown  in  Figure  (9-b).  indicates  that  the 
adaptation  results  in  improving  the  system  performance 
considerably.  Quantitative  comparison  between  the  two 
methods  can  be  obtained  by  considering  the  displacement 
index  (J^  defined  as  the  summation  of  the  squared 
position  error  of  the  beam  tip  over  time: 


t 


t=o 


and  the  control  votage  index  given  by 

t 

Ur  =  ^  v2(t)  . 

1=0 


(36} 
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Figure  (15)  -  Time  history  of  the  control  parameters 
for  the  variable  mass  beam 
a.  gain  ^  b.  time  constant  Tj  2 
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Figure  (16)  -  Time  response  of  the  beam  during  shape 

control  process 

a.  lip  displacement,  b.  control  voltage 
and  c.  modal  displacement 
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Figure  (17)  -  Time  response  of  the  beam  when 
controlled  with  MPPF  method 


where  t  is  the  time  of  each  experimental  run.  Table  4 
lists  the  displacement  and  control  voltage  indices  for 
the  different  algorithms. 

The  experimental  results  obtained  indicate  the 
effectiveness  of  the  AMPPF  algorithm  in  suppressing 
structural  vibrations  and  controlling  the  shape  of 
flexible  structures  in  the  presence  of  modeling  errors 
and  structural  system  changes.  Comparisons,  carried  out 
between  the  new  algorithm  and  the  non-adaptive  MPPF 
method,  emphasizes  its  favorable  vibration  damping  and 
shape  control  characteristics. 

5.  CONCLUSIONS 

This  study  has  presented  an  adaptive  modal  control 
algorithm  which  is  based  on  the  Positive  Position 
Feedback  method.  The  algorithm  utilizes  only  modal 
position  signals,  fed  through  first  order  filters,  to 
damp  out  the  vibration  of  undamped  flexible  systems. 
The  theory  behind  the  algorithm  is  presented.  The 
adaptation  laws  necessary  for  tuning  the  filters  and 
the  controller  gains  are  obtained,  for  all  the 
controlled  modes,  using  Lyapunov’s  stability  theory. 
The  algorithm  is  validated  experimentally  using  a 
single  piezo-electric  actuator  to  control  the  vibration 
and  shape  of  a  flexible  cantilevered  beam.  The 
results  obtained  indicate  the  effectiveness  of  the 
adaptation  in  improving  the  performance  of  the 
controller  particularly  when  it  is  used  to  control 
beams  with  varying  structural  parameters. 

The  study  demonstrates  clearly  the  simplicity  and 
potential  of  the  method  as  an  effective  method  for 
controlling  large  number  of  vibration  modes  with  a 
smaller  number  of  actuators  These  features  have 
important  practical  implications  that  make  the 
algorithm  viable  for  controlling  large  structures  in 
real  time. 
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Ta0l«  1  -  Main  design  parameters  of  test  beam 


Leng-.h  i  Width  |  ThicXnsss 

j  (cml  [  !eml  1  (cm) 

Young's  modulus 
(GN/m^)  1 

density 

;  2b  3  I  0  075 

2.96  j 

1. 31 

ELEMENTS  OF  A„  and  P  MATRICES 
A.  1.  Elements  of  A^  matrix 

The  parameters  a,,  aj  and  aj  which  define  the  elements 
of  the  A„  matrix  are  obtained  by  matching  the 
coefficients  of  the  characteristic  equation  (25)  of  the 
reference  system  with  the  coefficients  of  the  optimal 
ITAE  characteristic  equation  (26).  This  gives  a;  as 
the  negative  real  solution  of  the  following  equation 


aj  a-  2.15  af  ♦  3.5  u\  a,*  4  u*  =  0  (A-1  ) 

where  is  the  natural  frequency  of  the  optimal 

reference  system. 

The  parameters  aj  and  a3  can  then  be  determined  from: 

a2  =  -1.75  Ug  *  aj  (A-2) 

and  aj  =  -2  /  a,  {A-3) 

A.  2.  Elements  of  P  matrix 


The  matrix  P  is  given  by 


•  Pt  Pz  Ps 

Pz  P4  Ps 

■  Ps  Ps  Pa 


(A-4) 


where  the  elements  p,  through  p^,  are  obtained  by 
solving  the  Lyapunov  equation  (23)  which  becomes: 
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Table  Z  -  Main  design  parameters  of  the  actuator 


1 

1  Length 
j  (cm  1 

width 

(cm) 

thicK 

( c.m  1 

1  charge  coeff. 
j  ( m/v ) 

max. volt 

{v/«i 1 ) 

Young's  !nod. 
(C.M/m^) 

density 

(gm/cm^) 

1  4.35 

1. 375 

0.  1 

j  lOOxiO"'^ 

25 

63 

7,3 

where  Q|j’s  are  the  elements  of  the  positive  definite 
matrix  Q  which  are  assumed  by  the  control  system 
designer.  Equation  (A-U,  (A-21,  (A-3)  and  (A-5)  can 
then  be  used  to  compute  the  elements  Pi’s  of  the  P 
matrix. 

NOMENCLATURE 


Table  3  -  Modal  characteristics  of  the  beam  system 


Mode  number 

’  ^ 

Theor  mode 

-  HZ 

2  25 

13.57 

|E/per  mode 

-  Hz 

2.  20 

13.75 

j  Error 

2.  22 

1.  10 

I - 

Modal  damping 

,  031 

.018 

Table  4  -  Displacement  and  control  voltage  Indices 


j  Algorithms 

Displacement  Indexlmm^i 

Voltage  Index  (volt^) 

I  AMPPP 

0  36 

27,  S3 

i  MPPF 

0  73 

36  97 

a,  2  3  parameters  of  reference  model 

A,’p’  reference  and  structure-filter  system 

matrices 

B„  p  reference  and  structure-filter  system 

input  matrices 

e,  e  error  vector  and  its  derivative 

e,  ith  error  component  (i=l-3) 

f,  modal  control  force  of  the  1th  mode 

f  vector  defined  by  equation  (16) 

F,.  physical  force  vector 

h.  li  positive  definite  function  defined  by 

equation  (19)  and  its  derivative. 

K,  controller  gain 

P  mat',  lx  satisfying  Lyapunov  equation 

p,  ith  '  lement  of  matrix  P 

Q  positive  definite  matrix  defined  in 

Lvapunov  equation  (23) 


s,  s  Laplace  and  normalized  Laplace  operators 

t  t  ime 

t"  duration  of  experimental  run 

Uj  P  displacement  and  control  effort  indices 

u,,U|,U|  modal  displacement,  velocity  and 
acceleration  of  ilh  mode 
modal  displacement,  velocity  and  accelerat¬ 
ion  of  Ith  mode  of  reference  model. 

Up. Up. Up  modal  displacement,  velocity  and  accelerat¬ 
ion  of  1th  mode  of  structure. 

Up,  reference  modal  di.splacement  of  ith  mode 

V  voltage  applied  across  actuator 

V.  V  Lyapunov  function  and  Its  derivative 

X,  2  position  and  velocity  state  variables. 

X„_p  state  vectors  of  reference  and 

structure-filter  systems 

Y, .  Y,  output  of  ith  mode  filter  and  its  derivative 

y^ip  transverse  displacement  of  beam  tip. 

Greek  Symbols 

a,  dimensionless  time  constant  of  filters 

T,  dimensionless  gain  of  the  controller 

time  constant  of  the  ith  mode 
parameter  error  of  A  matrix 
parameter  error  of  B  matrix 
Up  frequency  of  reference  model 

u,  frequency  of  the  ith  mode 
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OPTIMAL  VIBRATION  CONTROL 
WITH 

MODAL  POSITIVE  POSITION  FEEDBACK 
A. Baz  and  S.Poh 

Mechanical  Engineering  Department 
The  Catholic  University  of  America 
Washington,  DC  20064 

SUMMARY 

The  vibrations  of  flexible  structures  are  controlled  by  an  Optimal 
Modal  Positive  Position  Feedback  (OMPPF)  algorithm  whose  control  forces  are 
generated  by  only  using  modal  position  signals  to  provide  damping  action  to 
undamped  structural  modes.  The  sub-optimal  parameters  of  the  OMPPF 
controller  are  obtained  by  casting  the  synthesis  problem  as  an  optimal 
control  problem  with  incomplete  state  feedback.  The  performance  of  the 
OMPPF  algorithm  is  enhanced  by  augmenting  it  with  a  "time  sharing"  strategy 
to  share  a  small  number  of  actuators  between  larger  number  of  vibration 
modes.  The  effectiveness  of  the  algorithm  in  damping  out  the  vibration  of 
flexible  structures  is  validated  experimentally  using  a  cantilevered  beam 
whose  multi-modes  of  vibration  are  controlled  by  a  single  piezo-electric 
actuator.  Theoretical  performance  predictions  are  found  to  be  in  close 
agreement  with  experimental  results. 

KEY  WORDS:  Active  vibration  control 

Modal  Positive  Position  Feedback 
Incomplete  state  feedback 
time  sharing  of  actuators 


1.  INTRODUCTION 


Considerable  emphasis  has  been  placed,  during  the  past  few  years,  on 

actively  controlling  the  vibration  of  a  wide  variety  of  flexible  structures. 

Several  control  algorithms  have  been  considered  ranging  from  the  simple 

1  2 

velocity  feedback  control  law  ’  to  the  more  imaginative  methods  such  as 

3  4 

the  Independent  Modal  Space  Control  (IMSC)  of  Meirovitch  ’  and  the 

Positive  Position  Feedback  (PPF)  of  Goh  and  Caughey^.  1988,  Baz  and 
6  7 

co-workers  ’  modified  the  IMSC  to  account  for  the  control  spillover  and 
devised  a  time  sharing  strategy  to  share  small  number  of  actuators  between 
larger  number  of  modes.  The  Modified  IMSC  method  (MIMSC)  has  been  shown  to 
have  favorable  vibration  damping  characteristics  as  compared  to  the  IMSC  and 

o 

the  Pseudo-Inverse  (PI)  methods  .  However  the  IMSC,  PI,  MIMSC  and  other 

modal  control  methods,  rely  in  their  operation  on  feeding  back  both  the 

modal  position  and  velocity  signals  of  the  controlled  modes  to  achieve  the 

required  vibration  damping  Extraction  of  these  signals  from  physical 

measurements  is  both  time  consuming  and  computationally  intensive  especially 
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when  dealing  with  large  structures.  In  1989,  Baz  et.al.  ’  developed  the 
Modal  Positive  Position  Feedback  (MPPF)  method  to  combine  the  attractive 
attributes  of  the  IMSC,  MIMSC  and  the  PPF  methods.  In  the  MPPF  method,  the 
controller  is  designed  completely  in  the  independent  modal  space  with  the 
open-loop  equations  of  the  system  remaining  uncoupled  even  after  including 
the  modal  controller.  Also,  the  MPPF  uses  only  the  modal  position  signals 
to  obtain  stable  and  damped  performance.  Such  performance  is  attained  by 
positively  feeding  back  the  position  signals  through  tuned  first  order 
filters.  Closed-form  expressions  are  given  for  determining  the  controller 
gains  and  the  time  constants  of  the  filters  in  order  to  ensure  stability, 
zero  steady-state  errors  and  maximum  closed-loop  damping  ratio.  The 


attractive  features  of  the  MPPF  method  have  been  successfully  demonstrated 
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by  controlling  the  vibration  of  simple  beams  and  more  complex  structures 
However,  no  attempt  has  been  made  to  develop  optimal  control  strategies  to 
enable  the  selection  of  the  controller  parameters  in  such  a  way  that  weighs 
the  relative  merits  of  the  vibrational  energy  vis-a-vis  the  control  effort. 

It  is  therefore  the  goal  of  this  study  to  develop  such  strategies  and 
experimentally  evaluate  the  effectiveness  of  the  Optimal  Modal  Positive 
Position  Feedback  (OMPPF)  algorithm  in  controlling  the  vibration  of  a 
cantilevered  beam. 

This  paper  is  organized  in  five  sections.  A  brief  introduction  is 
given  in  section  1.  The  concept  of  the  OMPPF  algorithm  and  the  selection 
of  its  optimal  parameters  are  presented  in  sections  2  and  3  respectively. 
Section  4  includes  the  experimental  evaluation  of  the  algorithm  along  with 
comparisons  with  other  algorithms.  Section  5  summarizes  the  conclusions  of 
this  study. 


2.  THE  CONCEPT  OF  THE  "MPPF"  ALGORITHM 

The  proposed  method  can  be  clearly  understood  by  considering  the  block 
diagram  shown  in  Figure  (1).  In  the  figure,  the  controller  is  used  to 
control,  in  the  independent  modal  space,  the  ith  mode  of  an  undamped 
flexible  structure. 

The  controller  feeds  back  positively  the  modal  displacement  qj  through 
a  first  order  filter  that  has  a  time  constant  Tj.  The  filter  output  Yj  is 
amplified  by  a  proportional  controller  gain  Kj  which  is  set  equal  to  yj  Uj, 
where  Wj  is  the  natural  frequency  of  the  ith  mode.  The  amplified  signal  fj, 
i.e.  the  modal  control  action,  is  then  sent  to  control  the  ith  mode  of  the 


Figure  (1)  -  Block  diagram  of  the  MPPF  controller. 


structure. 


Mathematically,  the  interaction  between  the  structural  mode  and  the 
controller  can  be  described  as  follows  : 


The 

structure 

qi  +  qi  =  fi  =  Ti  <^1  Yj 

(1) 

The 

filter 

Ti  Yi  +  Yi=  qi 

(2) 

where  the  dots  denote  differentiation  with  respect  to  time  t. 

Defining  a  normalized  time  t*  =  Wj  t,  then  equations  (1)  and  (2)  reduce 
to: 


f* 


The 

structure 

+ 

II 

(3) 

The 

filter 

/ 

Vi 

+  a,Yi*  at  qi 

(4) 

where  the  primes  denote  differentiation  with  respect  to  the  normalized  time 
t*  and  a^  defines  a  dimensionless  time  constant  of  the  filter  which  is  given 
by 

ttj  =  1  /  (  Tj  Uj  )  (5) 

Equations  (3)  and  (4)  can  be  combined  in  the  following  equation: 

/  /  /  /  It 

qi  +  q^  =  u  =  -  tti  (1  -  y,)  qi"  q^  (6) 

The  selection  of  the  optimal  gain  and  dimensionless  time  constant  aj 
of  the  filter  can  be  achieved  by  considering  the  following  state-space 
representation  of  equation  (6): 


X 


A  X  +  B  u 


(7) 


I 
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and 

such  that 

where 

X 

with 

Considering  the  output  matrix  C,  it  is  evident  that  the  above  state 
representation  conforms  with  the  premise  that  only  the  modal  displacement  qj 
is  accessible  for  measurements.  Also,  the  right  hand  side  of  equation  (6) 
and  equation  (9)  define  a  fictitious  control  action  u  which  is  based  on  an 
incomplete  state  feedback  as  the  second  element  of  the  gain  matrix  D  is  zero 
as  indicated  by  equation  (10).  Once  the  gains  dj  and  62  of  the  fictitious 
controller  are  optimally  determined,  as  will  described  in  section  3,  the 
gain  and  the  dimensionless  time  constant  ttj  of  the  original  modal 
controller  can  be  realized  and  the  algorithm  can  be  implemented  according  to 
equations  (1)  and  (2). 


3.  SELECTION  OF  THE  OPTIMAL  CONTROL  PARAMETERS 


I 


» 


3.1.  Formulation  of  the  optimal  control  problem 

The  optimal  gains  dj  and  dg  are  determined  by  considering 
first  the  following  optimal  control  problem: 


-  6  - 
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In  the  above  problem,  it  can  be  easily  shovm  ’  that  the  performance 
index  J  can  be  expressed  as  a  quadratic  form  of  the  initial  state  as 
follows 

J  =  xj  P(0)  Xg  (13) 

where  P  is  symmetric  and  positive-definite  matrix  which  is  a  solution  of  the 
following  Lyapunov  matrix  equation: 

P  =  (A  +  B  D)^  P  +  P  (A  +  B  D)  +  Q  C  +  R  D  (14) 

such  that  P(T)  =  0.  For  time  invarying  and  stable  systems,  P  approaches  a 
constant  value  determined  by  setting  P  =  0.  Equation  (14)  is  obtained  by 
differentiating  the  performance  index  with  respect  to  time  and  using 
equations  (7)  through  (9). 

It  is  important  to  note  that  the  above  optimal  control  problem,  given 
by  equation  (12),  can  not  be  directly  solved  using  the  solution  (D  =  -  R"* 
B^  P)  of  the  classical  linear  quadratic  regulator  problem^ ^  as  the  control 


action  u  is  based  on  an  incomplete  state  feedback. 


For  this  class  of 


problems,  the  solution  is  found  to  depend  on  the  initial  state  of  the 
12-14 

system  This  results  in  time  varying  feedback  which  is  not  efficient 

to  implement  particularly  for  large  systems.  Furthermore,  it  is  not  at  all 
certain  that  a  system  optimized  for  one  initial  disturbaince  will  perform 
satisfactorily  for  another  disturbance.  In  order  to  ensure  the 

acceptability  of  the  system  performance  for  all  initial  conditions,  the 
performance  of  the  worst  case  is  optimized.  This  is  achieved  by  considering 
the  following  performance  index  M  instead  of  the  original  performance  index 
J: 

M  =  max  J  /  (  Xo  X„  )  (15) 


This  performance  index  normalizes  the  original  performance  index  J  with 
respect  to  the  initial  state  of  the  system  X^,.  Substituting 
equation  (13)  into  equation  (15)  gives: 


M  =  max  (  x!  P(0)  X„  )  /  (  x!  X„  ) 


(16) 


It  can  be  shown,  as  outlined  in  appendix  A,  that  M  is  equal  to  the  maximiM 
eigen  value  of  the  P  matrix 


M  =  (P) 


(17) 


Therefore,  the  normalized  performance  index  M  becomes  independent  of  the 
initial  states  of  the  system  and  the  optimal  control  problem,  given  by 
equation  (12),  can  be  rewritten  as 


Find  the  controller  gains  dj  and  d2 
to  minimize  the  performance  index  M 
H  =  ^nax  (P) 

such  that: 

(A  +  B  D)"^  P  +  P  (A  +  B  D)  +  Q  C  +  R  D 


(18) 
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where  A,  B,  C  and  D  are  as  given  in  equation  (10). 


In  other  words,  the  above  optimal  control  problem  reduces  to  finding 
the  sub-optimal  fictitious  controller  gains  dj  and  d2  which  minimize  the 
maximum  eigen  value  of  the  matrix  P  that  satisfies  the  steady-state  Lyapunov 
equation.  The  solution  of  such  min-max  problem  yields  dj  and  d2  which  can 
be  used  along  with  equation  (11)  to  compute  the  optimal  gain  and  filter 
time  constant  aj  of  the  MPPF  algorithm. 

3.2.  Solution  of  the  optimal  control  problem 

The  solution  of  the  optimal  control  problem,  given  by  equation  (18),  is 

carried  out  according  to  the  flow  chart  shown  in  Figure  (2).  For  a  given 

weighting  parameter  R  and  an  initial  guess  of  the  control  gains  dj  and  d2, 

the  elements  of  the  matrix  P  are  computed  as  outlined  in  appendix  B.  The 

positive  definiteness  of  the  matrix  P  is  checked  and  its  eigen  values  Xj  are 

computed  using  the  Jacobi  method^^.  The  eigen  values  computed  are  ordered 

to  determine  the  maximum  value  .  This  value  is  minimized  by  an 

J 

optimization  routine  based  on  Powell’s  conjugate  direction  method'^  .  The 

subroutine  modifies  the  initial  guess  by  finding  an  improved  combination  of 

dj  and  d2.  The  improved  gains  are  used  again  to  compute  the  elements  of  P 

and  find  the  corresponding  X^^^  .  If  the  resulting  X^^^^  is  smaller 

j+i 

than  the  initial  value  X^^^  ,  the  process  is  repeated  again.  When  no 

J 

further  improvement  can  be  attained,  the  optimum  is  reached. 


3.3.  Optimal  parameters 


Figure  (3)  shows  the  contours  of  constant  maximum  eigen  values 
of  the  matrix  P  drawn,  as  a  function  of  dj  and  d2,  in  the  yj  and  aj 
plane  when  the  weighting  parameter  R  =  0.01.  The  figure  shows  also  the 
path  followed  from  an  initial  guess  (yj=0.  9  and  a^=2.5)  until  the  optimum  is 
attained.  Table  1  summarizes  the  optimal  values  of  and  aj  for  different 
values  of  the  weighting  parameter  R. 

Table  1  indicates  that  the  optimal  value  of  the  control  gain 
remains  nearly  at  0.5  irrespective  of  the  value  of  the  weighting  parameter 
R.  This  optimal  value  ensures  both  stability  and  zero  steady-state  error  of 
the  controller  as  can  be  seen  from  equations  (1)  and  (2).  However,  the 

optimal  values  of  the  dimensionless  filter  time  constant  aj  are  found  to 
decrease  as  the  weighting  parameter  R  is  increased.  For  a  given  mode  of 
vibration  Wp  this  means  that  the  filter  time  constant  Tj  increases  with 
increasing  R  as  indicated  by  equation  (5).  Therefore,  if  attenuation  of  the 
vibration  is  weighted  to  be  more  important  than  the  control  energy,  i.e.  R 
is  small,  then  the  filter  should  be  very  fast  and  Xj  should  be  small. 
Conversely,  when  the  control  energy  is  considered  to  be  more  important  then 
the  damping  of  the  vibration,  i.e.  R  is  large,  then  the  filter  is  tuned  to 
be  slow  and  its  time  constant  Xj  is  accordingly  large. 

It  is  very  important  to  note  that  the  presented  analysis  being 
applicable  to  any  mode,  results  in  the  same  optimal  tuning  conditions  for 
all  the  modes  of  vibration.  Accordingly,  optimal  vibration  control  can  be 
achieved  for  any  mode  i,  that  has  natural  frequency  Uj,  by  selecting  the 
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optiaial  control  gain  jtj  to  be  equal  to  0.5  and  the  optimal  time  constant  Tj 
of  its  filter  such  that 

Tj  =  ttj  /  Uj  for  i=l,..,  N  (19) 

where  N  is  the  number  of  controlled  modes  and  is  the  optimal 

dimensionless  filter  time  constant  given  in  Table  1. 

Figure  (4)  shows  a  flow  chart  of  the  proposed  OMPPF  algorithm.  The 

algorithm  performance  is  enhanced  by  augmenting  it  with  a  "time  sharing" 

strategy  to  share  a  small  number  of  actuators  to  control  larger  number  of 

9-11 

vibration  modes.  Such  "time  sharing"  strategy  has  been  shown  to  be  very 

effective  in  controlling  large  structures  with  a  small  number  of  actuators. 

3  4 

This  is  unlike  the  classical  IMSC  method  of  Meirovitch  ’  where  the  number 
of  actuators  needed  must  be  equal  to  the  number  of  modes  to  be  controlled. 
The  effectiveness  of  the  OMPPF  algorithm  in  damping  the  vibration  of 
flexible  systems  is  validated  experimentally  in  the  what  follows. 


4.  EXPERIMENTAL  PERFORMANCE  OF  OPTIMAL  MPPF  ALGORITHM 

4.1  EXPERIMENTAL  SET-UP 

A  thin  rectangular  cantilevered  beam  is  constructed  to  validate  the 
developed  algorithm  .  The  design  parameters  of  the  beam  are  given  in  Table 
2.  The  beam  is  controlled  by  one  piezo-electric  bimorph  made  from  G1195- 
ceramic.  The  actuator  is  available  commercially  (model  number  R205  )  from 
Piezo-Electric  Products,  Inc. ,  Metuchen,  NJ  08840-4015  .  Table  3  lists  the 
main  design  parameters  of  the  actuator. 


-  13  - 


Figure  (4)  -  Flow  chart  of  the  OMPPF  algorithm. 


The  experimental  beam  and  the  piezo-actuator  are  arranged  as  shown  in 
Figure  (5).  The  beam  is  divided  into  three  finite  elements  .  Bonded  to  the 
first  element,  near  the  fixed  end  of  the  beam,  is  the  piezo-actuator.  Three 
non-contacting  position  sensors  are  used  to  monitor  the  physical 

displacements  of  the  three  nodes  in  the  transverse  direction.  The  position 
signals  are  sampled  by  a  386-based  micro-processor  provided  with  an 
input/output  board  which  has  a  conversion  time  of  15  ps  and  a  resolution  of 
12  bits.  The  board  analog  outputs  have  a  settling  time  of  30  ps. 

The  micro-processor  uses  the  three  sampled  signals  to  compute  the  beam 
angular  deflections  and  the  linear  and  angular  velocities  of  the  nodes  .  The 
computed  state  variables  are  used  to  calculate  the  modal  coordinates  of  the 
flexible  system,  the  mode  that  has  the  highest  modal  energy,  the 

corresponding  optimal  modal  control  force  fj,  the  physical  control  force 
Fq  and  the  necessary  voltage  v  to  be  sent  to  the  piezo-actuator.  The 
implementation  of  these  calculations  ,  i.e.  the  OMPPF  algorithm  ,  is  carried 
out  in  real  time  in  3.04  ms. 

4. 2  Modal  characteristics  of  the  beam-actuator  system 

The  modal  characteristics  of  the  experimental  beam  are  determined 

9 

theoretically  and  validated  experimentally  using  classical  modal  analysis 

technique.  A  comparison  between  the  theoretical  and  the  experimental 

values  of  the  first  five  vibration  modes  of  the  beam-actuator  system  is 

given  in  Table  4.  The  table  gives  also  the  modal  damping  as  calculated  from 

17 

the  experimental  results  using  the  half  power  approach 

4.3  Experimental  results 


In  all  the  experiments  conducted  in  this  study,  the  beam  is 
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excited  at  its  second  mode  of  vibration  by  applying  sinusoidal  excitations 
of  20  volts  in  magnitude  to  the  piezo-actuator.  The  excitations  are 
maintained  for  a  period  of  0.3  seconds.  The  beam  is  either  left  to  vibrate 
freely  (i.  e.  uncontrolled)  or  under  the  action  of  one  modal  control  algorithm 
or  another.  The  uncontrolled  tip  displacement,  shovm  in  Figure  (6-a),  is 
used  as  a  datum  for  judging  the  effectiveness  of  the  control  algorithm. 

Figures  (6-b)  and  (6-c)  show  the  corresponding  time  response  of  the 
beam  when  it  is  optimally  controlled  by  the  OMPPF  algorithm  with  the 
weighting  parameter  R  set  at  0.01  and  the  piezo-actuator  is  dedicated  either 
to  the  first  mode  alone  or  time  shared  between  all  the  modes  respectively. 
It  is  evident  that  time  sharing  the  actuator  is  more  effective  in  damping 
out  the  beam  vibration.  Such  effectiveness  comes  about  because  the 
actuator  is  dedicated  to  control  the  mode  that  has  the  highest  instantaneous 
modal  energy.  This  is  not  necessarily  the  lowest  mode  of  vibration  as  it 
depends  on  the  nature  of  the  external  disturbance. 

Figure  (7-a)  and  (7-b)  show  the  corresponding  optimal  control  voltage 
sent  to  the  piezo-actuator  when  it  is  used  to  control  the  lowest  mode  or 
when  it  is  time  shared  between  the  modes  respectively.  The  figures  indicate 
that  higher  control  voltage  is  needed  when  the  actuator  is  time-shared 
between  modes.  This  is  attributed  to  the  need  for  fast  response  filters  to 
control  higher  order  modes,  as  implied  by  equation  (5),  which  results  in 
turn  in  higher  control  actions  as  indicated  by  equation  (6). 

The  effect  of  varying  the  weighting  parameter  R  on  the  performance  of 
the  OMPPF  algorithm,  with  its  time  sharing  capability,  is  shown  in  Figure 
(8).  Increasing  R  results  in  slower  vibration  damping  as  more  emphasis  is 


placed  on  the  importance  of  the  control  effort  over  the  vibration  energy. 
This  emphasis  is  manifested  by  the  decrease  in  the  control  voltages  as  R  is 
increased  as  indicated  in  Figure  (9). 

A  better  insight  into  the  effectiveness  of  the  OMPPF  algorithm  can  be 
gained  by  considering  the  Fast  Fourier  Transform  (FFT)  of  the  beam 
response.  Figure  (10)  shows  the  frequency  content  of  the  response  of  the 
uncontrolled  beam  in  comparison  with  the  optimally  controlled  beam.  These 
characteristics  are  obtained  by  sampling  the  beam  tip  position  signal  by  a 
spectrum  analyzer  and  performing  on  it  an  FFT  analysis.  The  Figure 
emphasizes  the  effectiveness  of  the  new  algorithm  particularly  when  R  is 
increased  and  when  it  is  augmented  with  the  time  sharing  capability. 

Comparisons  between  the  theoretical  and  experimental  time  responses  of 
the  uncontrolled  and  optimally  controlled  beam  (with  R=0.01)  are  shown  in 
Figures  (11-a)  and  (11-b)  respectively.  The  displayed  results  show  close 
agreements  between  theory  and  experiments. 


5.  CONCLUSIONS 

This  study  has  presented  an  optimal  modal  control  algorithm  which  is 
based  on  the  Positive  Position  Feedback  method.  The  algorithm  utilizes  only 
modal  position  signals,  fed  through  first  order  filters,  to  damp  out  the 
vibration  of  undamped  flexible  systems.  The  theory  behind  the  algorithm  is 
presented.  Optimal  tuning  of  the  filters  and  the  controller  gains  are 
obtained,  for  all  the  controlled  modes,  using  optimal  control  theory  for 
systems  with  incomplete  state  feedback.  The  algorithm  is  validated 
experimentally  using  a  single  piezo-electric  actuator  to  control  the 
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Figure  (10)-  Frequency  response  of  optimally  controlled  beam 

a.  at  different  R 

b.  with  and  without  time  sharing 


Displacement  (mm)  Displacement  (mm) 


vibration  of  a  flexible  cantilevered  beam. 


The  results  obtained  indicate 


its  effectiveness  in  suppressing  structural  vibration  particularly  when  it 
is  provided  with  time  sharing  capabilities. 

The  study  demonstrates  clearly  the  simplicity  and  potential  of  the 
method  as  an  effective  method  for  controlling  large  number  of  vibration 
modes  with  a  smaller  number  of  actuators  .  These  features  have  important 
practical  implications  that  make  the  algorithm  viable  for  controlling  large 
structures  in  real  time. 
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APPENDIX-A 


The  Performance  Index  (M) 

The  performance  index  M  is  given  by 

M  =  nax  J  /  (  Xq  X„  )  (A-1) 

*0 

This  performance  index  normalizes  the  original  performance  index  J  with 
respect  to  the  initial  state  of  the  system  X^.  Substituting 

equation  (13)  into  equation  (A-1)  gives: 

M  =  max  (  X^  P(0)  X„  )  /  (  X^  X^  )  (A-2) 

*0 

As  P  is  symmetric,  it  has  real  eigen  values  Xj  and  the  corresponding  eigen 
vectors  may  be  formed  as  an  orthonormal  set.  If  V  is  the  eigen  vector 
matrix  ,  then 

V  =  I,  (A-3) 

and  V’r  P(0)  V  =  A  (A-4) 

where  A  is  a  diagnoal  matrix  with  elements  Xj.  Assume  that; 

=  y  z  (A-5) 

then, 

J  /  (  X^  X„  )  =  (  xJ  P(0)  Xo  )  /  (  Xo  Xo  ) 

=  (  z^  A  z  )  /  (  z  ) 

=  (  X,  zf  )  /  {  j]  zf  )  (A-6) 

but  as,  X„i„  [  zf  £  [  X,  z?  <  X„,^  [  zj  (A-7) 

then,  from  equations  (A-1),  (A-6)  and  (A-7)  we  have 

M  =  X„,,  (P)  (A-8) 
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APPENDIX-B 


ELEMENTS  OF  MATRIX  P 


The  matrix  P  is  given  by 


-  Pi 

P2 

P3  ■ 

p  = 

P2 

P4 

Ps 

^  P3 

Ps 

Ps 

(B-i: 


where  the  elements  pj  through  pg,  which  satisfy  the  steady-state  Lyapunov 
equation  (14),  are  given  by 


and 


Pi  “  P3  ~  *^i  Ps* 

(B-2) 

P2  =  d2  /  [  2di 

(l-di/d2)  1, 

(B-3) 

P3  =  -  (  1  /  di 

-  di  R  )  /  2, 

(B-4) 

P4  ■  Ps  ■  *^2  Ps 

-  P3* 

(B-5) 

Ps  ~  P2’ 

(B-6) 

Pg  =  -  (  d2  R  / 

2  -  Ps  /  d2  ) 

(B-7) 

with  R  denoting  a  scalar  weighting  factor. 


For  a  given  d^  and  dj,  Pg  and  P3  can  be  calculated  first,  then  P5 
followed  by  Pj,  pg  then  p^. 
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NOMENCLATURE 


A 

B 

C 

D 

di,  d2 
fi 


I 

J 

Ki 

M 

P 

Pi 

Q 

cii. 

R 

T 

t 

t» 

u 

V 


V 

X 

X„ 


Qi.q 


structure-filter  system  matrix  (3x3) 

input  matrix  (3x1) 

output  matrix  (1x3) 

fictitious  controller  gain  matrix 

fictitious  controller  gains 

controlled  modal  force  of  the  ith  mode 

Physical  force  vector 

identity  matrix  (3x3) 

performance  index 

controller  gain 

normalized  performance  index 

matrix  satisfying  Lyapunov  equation  (3x1) 

ith  element  of  matrix  P 

weighting  matrix 

modal  displacement,  velocity  and  acceleration  of  ith  mode 
weighting  parameter 
final  time 
time 

dimensionless  time 
control  effort 

voltage  applied  across  actuator 

ith  eigen  vector  of  matrix  P 

eigen  vector  matrix  of  matrix  P 

state  vector  of  structure-filter  system  (3x1) 

initial  state  vector 

system  output 

the  output  of  the  ith  filter 
any  vector  (3x1 ) 


Greek  Symbols 


“i 


^max’ \nln 

A 


T 

(ji 


i 

I 


dimensionless  time  constant  of  filters 
dimensionless  gain  of  the  controller 
maximum  and  minimum  eigenvalues  of  P  matrix 
matrix  of  eigen  values  of  P  (3x3) 
time  constant  of  the  ith  mode 
frequency  of  the  ith  mode 


Superscripts 


1  dots  on  letter  denote  differentiation  with  respect  to  time 

2  primes  on  letter  denote  differentiation  with  respect  to 

dimensionless  time 
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Table  1  -  optimal  values  of  and  ttj  for  different  R 


R 

0.01 

1.00 

10.00 

100.00 

“i 

1.  139 

1.043 

0.  598 

0. 1999 

0.  505 

0.  503 

0. 5001 

0. 4988 

Table  2  -  Main  design  parameters  of  test  beam 


Length 

(cm) 

Width 

(cm) 

Thickness 

(cm) 

Young’ s  modulus 
(GN/m^) 

density 
(gm/cm  ) 

25. IS 

3.  75 

0.075 

2.  96 

1. 31 

Table  3  -  Main  design  parameters  of  the  actuator 


width 

thick. 

charge  coeff. 

max. volt 

Young’ s  mod. 

density 

(cm) 

(cm) 

(m/v) 

(v/mil ) 

(GN/m^) 

(gm/cm^) 

4.85 

1. 375 

0. 1 

190x10'^^ 

25 

63 

7.8 
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Table  4  -  Modal  characteristics  of  the  beam  system 


Mode  number 

2 

3 

5 

Theor.  mode 

-  Hz 

1.62 

11.0 

64.7 

87.  4 

157.  7 

Exper.  mode 

-  Hz 

1.55 

10.  1 

49.  5 

80.  5 

167.  5 

Error 

-  •/. 

■B 

m 

-6.2 

Modal  damping 

.038 

.022 

.016 

.  015 

.  010 
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vibrutloa  of  the  beam  when,  all  the  four  NITINOL  fibers  are  acti-v^ated  or 
unactlvated.  It  is  evident  that  the  activation  of  the  NITINOL  fibers 
shifts  the  first  three  modes  of  vibration  to  higher  frequency  bands  and 
that  the  enhanced  damping  characteristics  of  the  heated  beam  result  in 
considerable  vibration  attenuation  for  all  the  modes. 

A  summary  of  the  amplitude  attenuation  obtained  at  different  flow 
speeds  is  given  in  Figure  C33).  Attenuations  of  about  4D?4  are  observed 
for  the  three  modes  of  vibration  over  the  speed  range  considered. 


w 


The  effect  of  acti-srating  different  number  of  fibers  on  the 
attenci^ion  of  the  flow-induced  vibrations  is  shown  in'  Figvire  (34) 
at  flow  speed  of  8.32  m/s.  Increasing  the  number  of  activated  fibers 
results  in  a  proportionate  reduction  in  the  amplitude  of  vibration. 

A  summary  of  the  effect  of  number  of  activated  NITINOL  fibers  on 
the  amplitude  of  vibration  normalized  with  respect  to  the  amplitude  of 
vibration  of  the  uncontrolled  beam,  is  shown  in  Figure  (35). 

The  results  obtained  demonstrate  the  effectiveness  of 
NITINOL-reinforced  composites  In  suppressing  flow-induced  vibrations 
over  a  wide  range  of  flow  speeds. 

5.  SUMMARY 


The  static,  dynamic  and  thermal  characteristics  .  of 
NITINOL-reinforced  composite  beams  have  been  presented.  The 
fundamental  issues  governing  the  behavior  of  this  new  class  of  SMART 
Composites  have  been  Introduced.  Applications  of  NITINOL  r«!:inforclng 
fibers  la  the  control  of  buckling  and  flow-induced  vibrations  are 
Successfully  demonstrated. 


